Skip to content
Snippets Groups Projects
Commit 855e15c4 authored by Jean-François Remacle's avatar Jean-François Remacle
Browse files

I did that commit because without it, the solver crashes : I initialize

the partition to 1.

In Blasius 2D , we control the number of elements.
parent f3a9cf59
No related branches found
No related tags found
No related merge requests found
......@@ -692,7 +692,7 @@ void meshMetric::scaleMetric( int nbElementsTarget,
m.eig(V,S);
S(0) = std::min(std::max(S(0),1/(hmax*hmax)),1/(hmin*hmin));
S(1) = std::min(std::max(S(1),1/(hmax*hmax)),1/(hmin*hmin));
S(2) = std::min(std::max(S(2),1/(hmax*hmax)),1/(hmin*hmin));
if (_dim == 3)S(2) = std::min(std::max(S(2),1/(hmax*hmax)),1/(hmin*hmin));
SVector3 t1 (V(0,0),V(1,0),V(2,0));
SVector3 t2 (V(0,1),V(1,1),V(2,1));
SVector3 t3 (V(0,2),V(1,2),V(2,2));
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment