From 8b677e0e882bceaa60eca07fa8a92b786de7f309 Mon Sep 17 00:00:00 2001 From: Arne Morten Kvarving Date: Tue, 15 Nov 2016 17:05:39 +0100 Subject: [PATCH] fixed: add element corners to boundary integrals in ASMu2D(mx) --- src/ASM/LR/ASMu2D.C | 2 ++ src/ASM/LR/ASMu2Dmx.C | 3 +++ 2 files changed, 5 insertions(+) diff --git a/src/ASM/LR/ASMu2D.C b/src/ASM/LR/ASMu2D.C index ed545820..35387607 100644 --- a/src/ASM/LR/ASMu2D.C +++ b/src/ASM/LR/ASMu2D.C @@ -1289,6 +1289,8 @@ bool ASMu2D::integrate (Integrand& integrand, int lIndex, // Get integration gauss points over this element this->getGaussPointParameters(gpar[t2-1],t2-1,nGP,iel,xg); + if (integrand.getIntegrandType() & Integrand::ELEMENT_CORNERS) + this->getElementCorners(iel,fe.XC); // --- Integration loop over all Gauss points along the edge ------------- diff --git a/src/ASM/LR/ASMu2Dmx.C b/src/ASM/LR/ASMu2Dmx.C index ec28eae4..72a6662d 100644 --- a/src/ASM/LR/ASMu2Dmx.C +++ b/src/ASM/LR/ASMu2Dmx.C @@ -445,6 +445,9 @@ bool ASMu2Dmx::integrate (Integrand& integrand, int lIndex, return false; } + if (integrand.getIntegrandType() & Integrand::ELEMENT_CORNERS) + this->getElementCorners(iel,fe.XC); + // Get integration gauss points over this element this->getGaussPointParameters(gpar[t2-1],t2-1,nGP,geoEl,xg);