No commit message

No commit message

benchmarks/3d/CubeAniso.geo

+Mesh.Algorithm = 7; //Bamg
+Mesh.Algorithm3D= 7; //mmg3d
+lc = 0.1;
+
+Point(1) = {0.0,0.0,0.0,lc};         
+Point(2) = {1,0.0,0.0,lc};         
+Point(3) = {1,1,0.0,lc};         
+Point(4) = {0,1,0.0,lc};         
+Line(1) = {4,3};         
+Line(2) = {3,2};         
+Line(3) = {2,1};         
+Line(4) = {1,4};         
+Line Loop(5) = {2,3,4,1};         
+Plane Surface(6) = {5};         
+Extrude Surface { 6, {0,0.0,1} }; 
+
+/*
+Field[1] = Attractor;
+Field[1].FacesList = {23};
+Field[1].NNodesByEdge = 100;
+
+Field[2] = Threshold;
+Field[2].IField = 1;
+Field[2].DistMax = 0.3;
+Field[2].DistMin = 0.05;
+Field[2].LcMax = 0.1;
+Field[2].LcMin = 0.01;
+Background Field = 2;
+
+Field[3] = BoundaryLayer;
+Field[3].IField = 1;
+Field[3].hfar = 0.1;
+Field[3].hwall_n = 0.005;
+Field[3].hwall_t = 0.1;
+Background Field = 3;
+*/
+
+/*
+Field[1] = BoundaryLayer;
+Field[1].FaceList = {23};
+Field[1].hfar = 0.01;
+Field[1].hwall_n = 0.005;
+Field[1].hwall_t = 0.1;
+Field[1].ratio = 1.0;
+Field[1].thickness = 0.5;
+Background Field = 1;
+*/
+
+Field[2] = MathEvalAniso;
+Field[2].m11 = "1./(0.1)^2";
+Field[2].m12 = "0";
+Field[2].m13 = "0";
+Field[2].m22 = "y+1/(0.02)^2";
+Field[2].m23 = "0";
+Field[2].m33 = "1/(0.1)^2";
+Background Field = 2;

benchmarks/3d/CubeAniso2.geo

+lc = 0.1;
+
+Mesh.RemeshParametrization=1; //(0) harmonic (1) conformal 
+Mesh.RemeshAlgorithm=1; //(0) nosplit (1) automatic (2) split only with metis ///Default: 1 
+
+Merge "CubeAniso.msh";
+CreateTopology;
+
+Compound Surface(106) = {6};
+Compound Surface(115) = {15};
+Compound Surface(119) = {19};
+Compound Surface(123) = {23};
+Compound Surface(127) = {27};
+Compound Surface(128) = {28};
+
+Surface Loop(100) = {106,115,119,123,127,128};
+Volume(1) = {100};
+
+Field[2] = MathEvalAniso;
+Field[2].m11 = "1./(0.1)^2";
+Field[2].m12 = "0";
+Field[2].m13 = "0";
+Field[2].m22 = "y+1/(0.02)^2";
+Field[2].m23 = "0";
+Field[2].m33 = "1/(0.1)^2";
+Background Field = 2;