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

ASMu3D::constrainFace

Browse Source 
use midpoint()
This commit is contained in:
Arne Morten Kvarving
2023-09-01 11:56:20 +02:00
parent 3ad0b872c8
commit 82ccdc49e3
1 changed files with 2 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
src/ASM/LR/ASMu3D.C
View File

@@ -531,14 +531,9 @@ void ASMu3D::constrainFace (int dir, bool open, int dof, int code, char basis)
de.MNPC.reserve(faceElements.size());
for (LR::Element* el : faceElements)
{
// for mixed FEM models, let MLGE point to the *geometry* index
// for mixed FEM models, let MLGE point to the integration basis
if (de.lr != this->lrspline.get())
{
double umid = (el->umax() + el->umin())/2.0;
double vmid = (el->vmax() + el->vmin())/2.0;
double wmid = (el->wmax() + el->wmin())/2.0;
de.MLGE.push_back(lrspline->getElementContaining(umid,vmid,wmid));
}
de.MLGE.push_back(lrspline->getElementContaining(el->midpoint()));
else
de.MLGE.push_back(el->getId());