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

fixed: add element corners to boundary integrals in ASMu2D(mx)

Browse Source
This commit is contained in:
Arne Morten Kvarving 2016-11-15 17:05:39 +01:00 committed by Kjetil Andre Johannessen
parent b2365bbeab
commit 8b677e0e88
2 changed files with 5 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/ASM/LR/ASMu2D.C
View File

@ -1289,6 +1289,8 @@ bool ASMu2D::integrate (Integrand& integrand, int lIndex,
// Get integration gauss points over this element // Get integration gauss points over this element
this->getGaussPointParameters(gpar[t2-1],t2-1,nGP,iel,xg); 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 ------------- // --- Integration loop over all Gauss points along the edge -------------

3
src/ASM/LR/ASMu2Dmx.C
View File

@ -445,6 +445,9 @@ bool ASMu2Dmx::integrate (Integrand& integrand, int lIndex,
return false; return false;
} }
if (integrand.getIntegrandType() & Integrand::ELEMENT_CORNERS)
this->getElementCorners(iel,fe.XC);
// Get integration gauss points over this element // Get integration gauss points over this element
this->getGaussPointParameters(gpar[t2-1],t2-1,nGP,geoEl,xg); this->getGaussPointParameters(gpar[t2-1],t2-1,nGP,geoEl,xg);