Commit ffe49eda authored by Timo Koch's avatar Timo Koch
Browse files

[1d3d][avg] Skip point source if the circle stencil is empty (zero coupling term)

In this case the integration element is zero anyways and there is no coupling  but the pressure may
become NaN and makes the simulation crash due to NaN residuals.
parent 41345cde
......@@ -253,6 +253,11 @@ public:
}
}
// if the circle stencil is empty (that is the circle is entirely outside the domain)
// we do not add a (zero) point source
if (circleStencil.empty())
continue;
// export low dim circle stencil
if constexpr (isBox<bulkIdx>())
{
......
Supports Markdown
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