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

ASMu3D::evalPoint: use geometry basis

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

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

@ -1389,10 +1389,11 @@ bool ASMu3D::diracPoint (Integrand& integrand, GlobalIntegral& glInt,
int ASMu3D::evalPoint (const double* xi, double* param, Vec3& X) const
{
if (!lrspline) return -3;
const LR::LRSplineVolume* geo = this->getBasis(ASM::GEOMETRY_BASIS);
if (!geo) return -3;
for (int i = 0; i < 3; i++)
param[i] = (1.0-xi[i])*lrspline->startparam(i) + xi[i]*lrspline->endparam(i);
param[i] = (1.0-xi[i])*geo->startparam(i) + xi[i]*geo->endparam(i);
int iel = 0;
return this->evalPoint(iel,param,X);
@ -1401,9 +1402,10 @@ int ASMu3D::evalPoint (const double* xi, double* param, Vec3& X) const
int ASMu3D::evalPoint (int, const double* param, Vec3& X) const
{
const LR::LRSplineVolume* geo = this->getBasis(ASM::GEOMETRY_BASIS);
Go::Point X0;
lrspline->point(X0,param[0],param[1],param[2]);
for (int i = 0; i < 3 && i < lrspline->dimension(); i++)
geo->point(X0,param[0],param[1],param[2]);
for (int i = 0; i < 3 && i < geo->dimension(); i++)
X[i] = X0[i];
return 0;