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

ASMu3D::faceL2projection: add support for a separate geometry basis

Browse Source
This commit is contained in:
Arne Morten Kvarving 2023-09-01 10:06:41 +02:00
parent 0576f121b0
commit 494a8df6b8
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/ASM/LR/ASMu3Drecovery.C
View File

@ -500,6 +500,8 @@ bool ASMu3D::faceL2projection (const DirichletFace& face,
StdVector B(n*m);
A.resize(n,n);
const LR::LRSplineVolume* geo = this->getBasis(ASM::GEOMETRY_BASIS);
// find the normal direction for the face
int faceDir;
switch (face.edg)
@ -533,6 +535,7 @@ bool ASMu3D::faceL2projection (const DirichletFace& face,
for (size_t ie = 0; ie < face.MLGE.size(); ie++) // for all face elements
{
int iel = 1 + face.MLGE[ie];
int ielG = geo == lrspline.get() ? iel : geo->getElementContaining(lrspline->getElement(iel-1)->midpoint()) + 1;
std::array<Vector,3> gpar;
for (int d = 0; d < 3; d++)
@ -579,7 +582,7 @@ bool ASMu3D::faceL2projection (const DirichletFace& face,
if (gpar[2].size() > 1) w = param[2] = gpar[2](k3+1);
// Evaluate basis function derivatives at integration points
this->evaluateBasis(iel-1, u, v, w, N, dNdu, myGeoBasis);
this->evaluateBasis(ielG-1, u, v, w, N, dNdu, ASM::GEOMETRY_BASIS);
// Compute basis function derivatives
double dJxW = dA*wg[i]*wg[j]*utl::Jacobian(Jac,X,dNdX,Xnod,dNdu,t1,t2);
@ -590,7 +593,7 @@ bool ASMu3D::faceL2projection (const DirichletFace& face,
X.t = time;
// For mixed basis, we need to compute functions separate from geometry
if (face.lr != lrspline.get())
if (face.lr != geo)
{
// different lrspline instances enumerate elements differently
Go::BasisDerivs spline;