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

ASMs3D::evalPoint: use geometry basis

Browse Source
This commit is contained in:
Arne Morten Kvarving 2023-09-07 10:59:14 +02:00
parent 7ebdd7acab
commit 91959fdf27
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/ASM/ASMs3D.C
View File

@ -2982,12 +2982,13 @@ bool ASMs3D::integrateEdge (Integrand& integrand, int lEdge,
int ASMs3D::evalPoint (const double* xi, double* param, Vec3& X) const
{
if (!svol) return -3;
const Go::SplineVolume* geo = this->getBasis(ASM::GEOMETRY_BASIS);
if (!geo) return -3;
for (int i = 0; i < 3; i++)
param[i] = (1.0-xi[i])*svol->startparam(i) + xi[i]*svol->endparam(i);
param[i] = (1.0-xi[i])*geo->startparam(i) + xi[i]*geo->endparam(i);
SplineUtils::point(X,param[0],param[1],param[2],svol);
SplineUtils::point(X,param[0],param[1],param[2],geo);
// Check if this point matches any of the control points (nodes)
return this->searchCtrlPt(svol->coefs_begin(),svol->coefs_end(),