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

ASMs2D::evalPoint: use geometry basis

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

9
src/ASM/ASMs2D.C
View File

@ -2546,11 +2546,12 @@ bool ASMs2D::integrate (Integrand& integrand, int lIndex,
int ASMs2D::evalPoint (const double* xi, double* param, Vec3& X) const
{
if (!surf) return -2;
const Go::SplineSurface* geo = this->getBasis(ASM::GEOMETRY_BASIS);
if (!geo) return -2;
param[0] = (1.0-xi[0])*surf->startparam_u() + xi[0]*surf->endparam_u();
param[1] = (1.0-xi[1])*surf->startparam_v() + xi[1]*surf->endparam_v();
SplineUtils::point(X,param[0],param[1],surf);
param[0] = (1.0-xi[0])*geo->startparam_u() + xi[0]*geo->endparam_u();
param[1] = (1.0-xi[1])*geo->startparam_v() + xi[1]*geo->endparam_v();
SplineUtils::point(X,param[0],param[1],geo);
// Check if this point matches any of the control points (nodes)
return this->searchCtrlPt(surf->coefs_begin(),surf->coefs_end(),