diff --git a/tutorial/t13.geo b/tutorial/t13.geo index d03731080ef4429d89debc58fb177f02122368a3..4e2ca14a4f1072afcd873547ad4803ba3d92a3fa 100644 --- a/tutorial/t13.geo +++ b/tutorial/t13.geo @@ -41,20 +41,17 @@ Volume(1) = {1}; Physical Surface(1) = {s : s + #ss[]-1}; Physical Volume(1) = 1; -uniform = 1; -If(uniform) - // uniform mesh size... - Mesh.CharacteristicLengthMin = 2.5; - Mesh.CharacteristicLengthMax = 2.5; -EndIf -If(!uniform) - // ... or apply a funny mesh size field, just because we can :-) - Field[1] = MathEval; +// element size imposed by a size field +Field[1] = MathEval; +Field[1].F = "2.5"; +Background Field = 1; + +DefineConstant[ funny = {0, Choices{0,1}, Name "Apply funny mesh size field?"} ]; +If(funny) Field[1].F = "2*Sin((x+y)/5) + 3"; - Background Field = 1; EndIf -Mesh.RemeshAlgorithm = 1; // (0) no split (1) automatic (2) automatic only with metis -Mesh.RemeshParametrization = 7; // (0) harmonic (1) conformal spectral (7) conformal finite element +Mesh.RemeshAlgorithm = 1; // automatic +Mesh.RemeshParametrization = 7; // conformal finite element Geometry.HideCompounds = 0; // don't hide the compound entities Mesh.Algorithm = 6; // Frontal