Skip to content
Snippets Groups Projects
Commit 40c57619 authored by Ujwal Kishore Jinaga's avatar Ujwal Kishore Jinaga :clown:
Browse files

corrected

parent 40e585e1
No related branches found
No related tags found
No related merge requests found
...@@ -583,8 +583,8 @@ void mlawNonLinearTVENonLinearTVP::predictorCorrector_TVP_nonAssociatedFlow_nonL ...@@ -583,8 +583,8 @@ void mlawNonLinearTVENonLinearTVP::predictorCorrector_TVP_nonAssociatedFlow_nonL
DdevPhidGamma_RHS(i,j) += ( 0.5*2*_b*Ke*1/q1->_Jacobian*(ptilde+Gamma*dPhiPdGamma)*_I(i,j) - 3*pow(kk,1)*Hb/Cxdev*devPhi(i,j) ); // pow(kk,2) DEBUG DdevPhidGamma_RHS(i,j) += ( 0.5*2*_b*Ke*1/q1->_Jacobian*(ptilde+Gamma*dPhiPdGamma)*_I(i,j) - 3*pow(kk,1)*Hb/Cxdev*devPhi(i,j) ); // pow(kk,2) DEBUG
for (int k=0; k<3; k++) for (int k=0; k<3; k++)
for (int l=0; l<3; l++){ for (int l=0; l<3; l++){
DdevPhidDgamma_RHS(i,j) += 0.5* 2.*_b/3.*Gamma*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l))*_I(k,l)*Hp*1/q1->_Jacobian; DdevPhidDgamma_RHS(i,j) += 0.5* 2.*_b/3.*Gamma*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l)*1/q1->_Jacobian)*_I(k,l)*Hp;
DdevPhidGamma_RHS(i,j) += 0.5* ( -3.*Ge_Tensor(i,j,k,l)*1/q1->_Jacobian*devPhi(k,l) + 2.*_b/3.*Gamma*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l))*_I(k,l)*dPhiPdGamma ); DdevPhidGamma_RHS(i,j) += 0.5* ( -3.*Ge_Tensor(i,j,k,l)*1/q1->_Jacobian*devPhi(k,l) + 2.*_b/3.*Gamma*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l)*1/q1->_Jacobian)*_I(k,l)*dPhiPdGamma );
DdevPhidGamma_RHS(i,j) += 0.5*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l)*1/q1->_Jacobian)*N(k,l); // DEBUG DdevPhidGamma_RHS(i,j) += 0.5*(-Ge_TrEeTensor(i,j,k,l)*1/q1->_Jacobian + Dho(i,j,k,l)*1/q1->_Jacobian)*N(k,l); // DEBUG
} }
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment