Skip to content
Snippets Groups Projects
Commit 23b191d0 authored by Robert K's avatar Robert K
Browse files

ScalingLimiter: relax physical check.

parent 2b8791f0
No related branches found
No related tags found
1 merge request!3Scaling limiter
Pipeline #
......@@ -530,27 +530,23 @@ namespace Fem
RangeType minVal( enVal );
RangeType maxVal( enVal );
// evaluate uEn on all quadrature points on the intersections
for (const auto& intersection : intersections(gridPart_, en) )
{
FaceQuadratureType faceQuadInner(gridPart_,intersection, faceQuadOrd_, FaceQuadratureType::INSIDE);
if( !checkPhysicalQuad( faceQuadInner, uEn, minVal, maxVal ) )
{
limiter = true;
}
}
CornerPointSetType cornerquad( en );
if( ! checkPhysicalQuad( cornerquad, uEn, minVal, maxVal ) )
{
limiter = true;
}
// evaluate uEn on all interior quadrature points
VolumeQuadratureType quad( en, spc_.order( en ) );
if( ! checkPhysicalQuad( quad, uEn, minVal, maxVal ) )
{
limiter = true;
// compute full set of min and max values
// evaluate uEn on all quadrature points on the intersections
for (const auto& intersection : intersections(gridPart_, en) )
{
FaceQuadratureType faceQuadInner(gridPart_,intersection, faceQuadOrd_, FaceQuadratureType::INSIDE);
checkPhysicalQuad( faceQuadInner, uEn, minVal, maxVal );
}
// evaluate uEn on all interior quadrature points
VolumeQuadratureType quad( en, spc_.order( en ) );
checkPhysicalQuad( quad, uEn, minVal, maxVal );
}
// scale function
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment