Commit 4e1d4427 authored by Benoît LÉ's avatar Benoît LÉ

[xFEM] Corrected bug in xGradOperatorAxisym and xGradLocalOperatorAxisym

Use geo_integ instead of geo_appro.
parent 0c0aaefa
......@@ -64,7 +64,7 @@ class xGradOperatorAxisym
vals_tmp.reserve(vals.size() + f->size());
xField<>::getFF(f->begin(), f->end(), vals_tmp, geo_appro, geo_integ, xtool::xIdentity<xtensor::xVector<>>());
auto uvw = geo_appro->getXYZ();
auto uvw = geo_integ->getXYZ();
double r = uvw(0);
for (int i = 0; i < vals.size() / 2; ++i)
{
......@@ -103,7 +103,7 @@ class xGradLocalOperatorAxisym
vals_tmp.reserve(vals.size() + f->size());
xField<>::getFF(f->begin(), f->end(), vals, geo_appro, geo_integ, xtool::xIdentity<xtensor::xVector<>>());
auto uvw = geo_appro->getXYZ();
auto uvw = geo_integ->getXYZ();
double r = uvw(0);
for (int i = 0; i < vals.size() / 2; ++i)
{
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment