OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fixed: set fe.iel in ASMsxDLag::evalSolution

Browse Source
This commit is contained in:
Arne Morten Kvarving 2020-11-10 11:50:18 +01:00
parent 5a671fa594
commit 6a2959d9ca
2 changed files with 4 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/ASM/ASMs2DLag.C
View File

@ -738,6 +738,8 @@ bool ASMs2DLag::evalSolution (Matrix& sField, const IntegrandBase& integrand,
if (!Lagrange::computeBasis(fe.N,dNdu,p1,fe.xi,p2,fe.eta)) if (!Lagrange::computeBasis(fe.N,dNdu,p1,fe.xi,p2,fe.eta))
return false; return false;
fe.iel = MLGE[iel-1];
// Compute the Jacobian inverse // Compute the Jacobian inverse
fe.detJxW = utl::Jacobian(Jac,fe.dNdX,Xnod,dNdu); fe.detJxW = utl::Jacobian(Jac,fe.dNdX,Xnod,dNdu);

2
src/ASM/ASMs3DLag.C
View File

@ -981,6 +981,8 @@ bool ASMs3DLag::evalSolution (Matrix& sField, const IntegrandBase& integrand,
if (!Lagrange::computeBasis(fe.N,dNdu,p1,fe.xi,p2,fe.eta,p3,fe.zeta)) if (!Lagrange::computeBasis(fe.N,dNdu,p1,fe.xi,p2,fe.eta,p3,fe.zeta))
return false; return false;
fe.iel = MLGE[iel-1];
// Compute the Jacobian inverse // Compute the Jacobian inverse
fe.detJxW = utl::Jacobian(Jac,fe.dNdX,Xnod,dNdu); fe.detJxW = utl::Jacobian(Jac,fe.dNdX,Xnod,dNdu);