From 6a2959d9ca90f9252b01fa7c5f373215ede07fc6 Mon Sep 17 00:00:00 2001 From: Arne Morten Kvarving Date: Tue, 10 Nov 2020 11:50:18 +0100 Subject: [PATCH] fixed: set fe.iel in ASMsxDLag::evalSolution --- src/ASM/ASMs2DLag.C | 2 ++ src/ASM/ASMs3DLag.C | 2 ++ 2 files changed, 4 insertions(+) diff --git a/src/ASM/ASMs2DLag.C b/src/ASM/ASMs2DLag.C index 39d18d0a8..6651e500f 100644 --- a/src/ASM/ASMs2DLag.C +++ b/src/ASM/ASMs2DLag.C @@ -738,6 +738,8 @@ bool ASMs2DLag::evalSolution (Matrix& sField, const IntegrandBase& integrand, if (!Lagrange::computeBasis(fe.N,dNdu,p1,fe.xi,p2,fe.eta)) return false; + fe.iel = MLGE[iel-1]; + // Compute the Jacobian inverse fe.detJxW = utl::Jacobian(Jac,fe.dNdX,Xnod,dNdu); diff --git a/src/ASM/ASMs3DLag.C b/src/ASM/ASMs3DLag.C index c3afa624d..e6d8d44e3 100644 --- a/src/ASM/ASMs3DLag.C +++ b/src/ASM/ASMs3DLag.C @@ -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)) return false; + fe.iel = MLGE[iel-1]; + // Compute the Jacobian inverse fe.detJxW = utl::Jacobian(Jac,fe.dNdX,Xnod,dNdu);