From e3caca11eb3816f94c87fa7e3c2f2fd2fd323a35 Mon Sep 17 00:00:00 2001 From: Arne Morten Kvarving Date: Thu, 7 Sep 2023 10:59:14 +0200 Subject: [PATCH] ASMu3D::evalPoint: use geometry basis --- src/ASM/LR/ASMu3D.C | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) diff --git a/src/ASM/LR/ASMu3D.C b/src/ASM/LR/ASMu3D.C index 9af6b784..60407553 100644 --- a/src/ASM/LR/ASMu3D.C +++ b/src/ASM/LR/ASMu3D.C @@ -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;