Commit cfe6ecb8 authored by Benoît LÉ's avatar Benoît LÉ

[xFEM] Corrected a few bugs in xEvalGradFieldAxisym

- Use integ instead of appro
- Deal with the case r = 0
parent d5c4a868
......@@ -409,12 +409,19 @@ class xEvalGradFieldAxisym : public xEval<typename UnaryOperator::result_type>
xtensor::xVector<> val;
f.getVal(appro, integ, val);
auto uvw = appro->getXYZ();
auto uvw = integ->getXYZ();
double r = uvw(0);
f.getGrad(appro, integ, v);
result = funct(v);
result(2, 2) = val(0) / r;
if (r > 0.)
{
result(2, 2) = val(0) / r;
}
else
{
result(2, 2) = result(0, 0);
}
}
private:
......
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