remove test-cases too large for nightly testing

benchmarks/2d_large/poisson_very_fine_mesh.geo

-Point (1) = {0, 0, 0, 1};
-Point (2) = {500, 0, 0, 1};
-Point (3) = {500, 25, 0, 1};
-Point (4) = {0, 25, 0, 1};
-Point (5) = {0, -2000, 0, 1};
-Point (6) = {2000, -2000, 0, 1};
-Point (7) = {2000, 0, 0, 1};
-Point (8) = {1500, 0, 0, 1};
-Point (9) = {1500, 25, 0, 1};
-Point (10) = {2000, 25, 0, 1};
-Point (11) = {2000, 50, 0, 1};
-Point (12) = {0, 50, 0, 1};
-Point (13) = {525, 50, 0, 1};
-Point (14) = {525, 25, 0, 1};
-Point (15) = {1475, 25, 0, 1};
-Point (16) = {1475, 50, 0, 1};
-Point (17) = {0, 2000, 0, 1};
-Point (18) = {2000, 2000, 0, 1};
-Line (1) = {1, 2};
-Line (2) = {2, 3};
-Line (3) = {3, 4};
-Line (4) = {4, 1};
-Line (5) = {7, 8};
-Line (6) = {9, 8};
-Line (7) = {9, 10};
-Line (8) = {10, 7};
-Line (9) = {10, 11};
-Line (10) = {11, 16};
-Line (11) = {16, 15};
-Line (12) = {15, 14};
-Line (13) = {14, 13};
-Line (14) = {13, 12};
-Line (15) = {12, 4};
-Line (16) = {2, 8};
-Line (17) = {5, 1};
-Line (18) = {5, 6};
-Line (19) = {6, 7};
-Line (20) = {17, 18};
-Line (21) = {18, 11};
-Line (22) = {17, 12};
-Line Loop (24) = {20, 21, 10, 11, 12, 13, 14, -22};
-Plane Surface (24) = {24};
-Line Loop (26) = {18, 19, 5, -16, -1, -17};
-Plane Surface (26) = {26};
-Line Loop (28) = {1, 2, 3, 4};
-Plane Surface (28) = {28};
-Line Loop (30) = {6, -5, -8, -7};
-Plane Surface (30) = {30};
-Line Loop (32) = {6, -16, 2, 3, -15, -14, -13, -12, -11, -10, -9, -7};
-Plane Surface (32) = {32};

benchmarks/3d_large/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 {0,0.0,1}{ Surface {6}; }
-
-/*
-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.005)^2";
-Field[2].m23 = "0";
-Field[2].m33 = "1/(0.1)^2";
-Background Field = 2;

benchmarks/3d_large/Falcon/CreateFullInitialMesh.geo

-Merge "InitialMeshFalcon.msh";
-Plugin(NewView).Run;
-Plugin(Transform).A33 = -1;
-Plugin(Transform).Tz = 2 * 0.272342;
-Plugin(Transform).SwapOrientation = 1;
-Plugin(Transform).Run;
-Physical Surface(1) = 1; // keep same num
-For i In {2 : 16}
-Physical Surface(100 + i) = i;
-EndFor
-Save "tmp.msh";
-
-Delete All;
-
-Mesh.SwitchElementTags = 1;
-Merge "tmp.msh";
-
-Mesh.SwitchElementTags = 0;
-Merge "InitialMeshFalcon.msh";
-
-Geometry.Tolerance = 1.e-5;
-Coherence Mesh;
-
-Mesh.SaveAll = 1;
-Save "FullInitialMeshFalcon.msh";

benchmarks/3d_large/Falcon/CreateVolumeAroundFalcon_FrontalISO.geo

-Mesh.Algorithm3D=4; //1-Delaunay, 4-Frontal, 5-Frontal Delaunay 6-Frontal hex, 7-MMG3D, 9-RTree
-
-lc = 2;
-Mesh.CharacteristicLengthMin = lc;
-Mesh.CharacteristicLengthMax = lc;
-
-//Merge "InitialMeshFalcon.stl";
-Merge "SurfaceMeshUniform.stl";
-//Merge "SurfaceMeshAnisoCurvature.stl";
-
-CreateTopology;
-
-z0 = 0.272342;
-zmax = 11;
-xmin = -20;
-xmax = 3;
-ymin = -3.0;
-ymax = 6.0;
-
-Point(20001) = {xmin, ymin, z0};
-Point(20002) = {xmax, ymin, z0};
-Point(20003) = {xmax, ymax, z0};
-Point(20004) = {xmin, ymax, z0};
-Point(20005) = {xmin, ymin, zmax};
-Point(20006) = {xmax, ymin, zmax};
-Point(20007) = {xmax, ymax, zmax};
-Point(20008) = {xmin, ymax, zmax};
-
-Line(10) = {20001, 20002};
-Line(11) = {20002, 20003};
-Line(12) = {20003, 20004};
-Line(13) = {20004, 20001};
-Line(14) = {20005, 20006};
-Line(15) = {20006, 20007};
-Line(16) = {20007, 20008};
-Line(17) = {20008, 20005};
-Line(18) = {20001, 20005};
-Line(19) = {20002, 20006};
-Line(20) = {20003, 20007};
-Line(21) = {20004, 20008};
-
-Line Loop(22) = {10,19,-14,-18};
-Line Loop(23) = {11,20,-15,-19};
-Line Loop(24) = {12,21,-16,-20};
-Line Loop(25) = {13,18,-17,-21};
-Line Loop(26) = {14,15,16,17};
-Line Loop(27) = {1,10,11,12,13};
-
-Plane Surface(2) = {22};
-Plane Surface(3) = {23};
-Plane Surface(4) = {24};
-Plane Surface(5) = {25};
-Plane Surface(6) = {26};
-Plane Surface(7) = {27};
-
-Surface Loop(8) = {1,2,3,4,5,6,7};
-Volume(1) = {8};
-
-Physical Surface(1) = {1};
-Physical Surface(2) = {5, 4, 3, 2, 6};
-Physical Surface(3) = {7};
-Physical Volume(4) = {1};

benchmarks/3d_large/Falcon/FullSurfaceMeshUniform.geo

-Mesh.RemeshParametrization=7; //0=harmonic_circle, 1=conformal, 2=rbf, 3=harmonic_plane, 4=convex_circle, 5=convex_plane, 6=harmonic square
-Mesh.RemeshAlgorithm=1; //(0) nosplit (1) automatic (2) split only with metis
-Mesh.Algorithm=6; //(1=MeshAdapt, 2=Automatic, 5=Delaunay, 6=Frontal, 7=bamg, 8=delquad)
-
-lc = 0.3;
-
-Mesh.CharacteristicLengthFromPoints = 0;
-Mesh.CharacteristicLengthMin = lc;
-Mesh.CharacteristicLengthMax = lc;
-Mesh.LcIntegrationPrecision = 1.e-5;
-Mesh.MinimumCirclePoints = 50;
-Mesh.CharacteristicLengthExtendFromBoundary = 0;
-Mesh.CharacteristicLengthFromCurvature = 0;
-Mesh.CharacteristicLengthFromPoints = 0;
-
-Merge "FullInitialMeshFalcon.msh";
-
-CreateTopology;
-
-// Make all lines compound:
-ll[] = Line "*";
-For j In {0 : #ll[]-1}
-  Compound Line(1000+j) = ll[j];
-  Physical Line(1000+j) = (1000+j);
-EndFor
-
-// Make all surfaces compound and physical:
-ss[] = Surface "*";
-For i In {0 : #ss[]-1}
-  Compound Surface(i+1000) = ss[i];
-  Physical Surface(i+1) = { i+1000 };
-EndFor
-

benchmarks/3d_large/Falcon/README.txt

-// for the symmetric falcon:
-
-gmsh falcon_surface.script - // generates remeshed falcon boundary mesh(es)
-gmsh falcon_volume.script -  // generates volume mesh in box around falcon(s)
-gmsh partition.script -  // partitions the mesh and saves using new physicals
-gmsh split.script - // splits the partitionned mesh into separate files
-
-// for the full falcon(s)
-
-gmsh CreateFullInitialMesh.geo - // generates the intial falcon boundary mesh
-gmsh fullfalcon_surface.script - // generates remeshed falcon boundary mesh(es)
-gmsh fullfalcon_volume.script - // generates volume mesh in box around falcon(s)
-gmsh partition.script -  // partitions the mesh and saves using new physicals
-gmsh fullsplit.script - // splits the partitionned mesh into separate files

benchmarks/3d_large/Falcon/SurfaceMeshAnisoCurvature.geo

-Mesh.RemeshParametrization=7; //0=harmonic_circle, 1=conformal, 2=rbf, 3=harmonic_plane, 4=convex_circle, 5=convex_plane, 6=harmonic square
-Mesh.RemeshAlgorithm=1; //(0) nosplit (1) automatic (2) split only with metis
-
-Mesh.Algorithm=7; //(1=MeshAdapt, 2=Automatic, 5=Delaunay, 6=Frontal, 7=bamg, 8=delquad)
-
-Mesh.CharacteristicLengthFromPoints = 0;
-Mesh.CharacteristicLengthMin = 0.005;
-Mesh.CharacteristicLengthMax = 2.0;
-Mesh.LcIntegrationPrecision=1.e-5;
-Mesh.MinimumCirclePoints=50;
-Mesh.CharacteristicLengthExtendFromBoundary=0;
-Mesh.CharacteristicLengthFromCurvature = 1;
-Mesh.CharacteristicLengthFromPoints = 0;
-
-Merge "InitialMeshFalcon.msh";
-
-
-CreateTopology;
-
-
-// Make all lines compound:
-ll[] = Line "*";
-For j In {0 : #ll[]-1}
-  Compound Line(1000+j) = ll[j];
-  Physical Line(1000+j) = (1000+j);
-EndFor
-
-
-
-// Make all surfaces compound and physical:
-ss[] = Surface "*";
-For i In {0 : #ss[]-1}
-  Compound Surface(i+100) = ss[i];
-  Physical Surface(i+1) = { i+100 };
-EndFor
-
-// The trailing edge of the horizontal tail surface
-// should be an attractor
-
-//Point(100000) = {16.858,2.0482,0.018793,0.1};
-Point(100000) = {-1.16411804,1.83824622,0.29113513};
-Field[1] = Attractor;
-Field[1].NodesList = { 100000 };
-
-// Local refinement near the airfoil:
-Field[2] = Threshold;
-Field[2].IField = 1;
-Field[2].DistMin = 0.05;
-Field[2].DistMax = 0.3;
-Field[2].LcMin = 0.05;
-Field[2].LcMax = 1.0;
-Field[2].StopAtDistMax = 1;
-
-Background Field = 2;
-
-
-
-

benchmarks/3d_large/Falcon/SurfaceMeshUniform.geo

-Mesh.RemeshParametrization=7; //0=harmonic_circle, 1=conformal, 2=rbf, 3=harmonic_plane, 4=convex_circle, 5=convex_plane, 6=harmonic square
-Mesh.RemeshAlgorithm=1; //(0) nosplit (1) automatic (2) split only with metis
-Mesh.Algorithm=6; //(1=MeshAdapt, 2=Automatic, 5=Delaunay, 6=Frontal, 7=bamg, 8=delquad)
-
-lc = 0.3;
-
-Mesh.CharacteristicLengthFromPoints = 0;
-Mesh.CharacteristicLengthMin = lc;
-Mesh.CharacteristicLengthMax = lc;
-Mesh.LcIntegrationPrecision = 1.e-5;
-Mesh.MinimumCirclePoints = 50;
-Mesh.CharacteristicLengthExtendFromBoundary = 0;
-Mesh.CharacteristicLengthFromCurvature = 0;
-Mesh.CharacteristicLengthFromPoints = 0;
-
-Merge "InitialMeshFalcon.msh";
-
-CreateTopology;
-
-// Make all lines compound:
-ll[] = Line "*";
-For j In {0 : #ll[]-1}
-  Compound Line(1000+j) = ll[j];
-  Physical Line(1000+j) = (1000+j);
-EndFor
-
-// Make all surfaces compound and physical:
-ss[] = Surface "*";
-For i In {0 : #ss[]-1}
-  Compound Surface(i+100) = ss[i];
-  Physical Surface(i+1) = { i+100 };
-EndFor
-

benchmarks/3d_large/Falcon/falcon_surface.script

-Include "param.geo";
-
-Mesh.RemeshParametrization=7;
-Mesh.RemeshAlgorithm=1;
-Mesh.Algorithm=6;
-
-Mesh.CharacteristicLengthFromPoints = 0;
-Mesh.CharacteristicLengthMin = lc;
-Mesh.CharacteristicLengthMax = lc;
-Mesh.LcIntegrationPrecision = 1.e-5;
-Mesh.MinimumCirclePoints = 50;
-Mesh.CharacteristicLengthExtendFromBoundary = 0;
-
-Merge "InitialMeshFalcon.msh";
-
-CreateTopology;
-
-// Make all lines compound:
-ll[] = Line "*";
-For j In {0 : #ll[]-1}
-  Compound Line(1000+j) = ll[j];
-  Physical Line(1000+j) = (1000+j);
-EndFor
-
-// Make all surfaces compound and physical:
-ss[] = Surface "*";
-For i In {0 : #ss[]-1}
-  Compound Surface(i+100) = ss[i];
-  Physical Surface(i+1) = { i+100 };
-EndFor
-
-Mesh 2;
-Save "tmp.stl";

benchmarks/3d_large/Falcon/falcon_volume.script

-Include "param.geo";
-
-Merge "tmp.stl";
-
-CreateTopology;
-
-z0 = 0.272342;
-zmax = 11;
-xmin = -20;
-xmax = 3;
-ymin = -3.0;
-ymax = 6.0;
-
-Point(20001) = {xmin, ymin, z0, lc2};
-Point(20002) = {xmax, ymin, z0, lc2};
-Point(20003) = {xmax, ymax, z0, lc2};
-Point(20004) = {xmin, ymax, z0, lc2};
-Point(20005) = {xmin, ymin, zmax, lc2};
-Point(20006) = {xmax, ymin, zmax, lc2};
-Point(20007) = {xmax, ymax, zmax, lc2};
-Point(20008) = {xmin, ymax, zmax, lc2};
-
-Line(10) = {20001, 20002};
-Line(11) = {20002, 20003};
-Line(12) = {20003, 20004};
-Line(13) = {20004, 20001};
-Line(14) = {20005, 20006};
-Line(15) = {20006, 20007};
-Line(16) = {20007, 20008};
-Line(17) = {20008, 20005};
-Line(18) = {20001, 20005};
-Line(19) = {20002, 20006};
-Line(20) = {20003, 20007};
-Line(21) = {20004, 20008};
-
-Line Loop(22) = {10,19,-14,-18};
-Line Loop(23) = {11,20,-15,-19};
-Line Loop(24) = {12,21,-16,-20};
-Line Loop(25) = {13,18,-17,-21};
-Line Loop(26) = {14,15,16,17};
-Line Loop(27) = {1,10,11,12,13};
-
-Plane Surface(2) = {22};
-Plane Surface(3) = {23};
-Plane Surface(4) = {24};
-Plane Surface(5) = {25};
-Plane Surface(6) = {26};
-Plane Surface(7) = {27};
-
-Surface Loop(8) = {1,2,3,4,5,6,7};
-Volume(1) = {8};
-
-Physical Surface(1) = {1};
-Physical Surface(2) = {5, 4, 3, 2, 6};
-Physical Surface(3) = {7};
-Physical Volume(4) = {1};
-
-Mesh 3;
-Save "falcon.msh";

benchmarks/3d_large/Falcon/fullfalcon_surface.script

-Include "param.geo";
-
-//Mesh.RemeshParametrization=7;
-//Mesh.RemeshAlgorithm=1;
-Mesh.Algorithm=6;
-
-Mesh.CharacteristicLengthFromPoints = 0;
-Mesh.CharacteristicLengthMin = lc;
-Mesh.CharacteristicLengthMax = lc;
-Mesh.LcIntegrationPrecision = 1.e-5;
-Mesh.MinimumCirclePoints = 50;
-Mesh.CharacteristicLengthExtendFromBoundary = 0;
-Mesh.CharacteristicLengthFromCurvature = 0;
-Mesh.CharacteristicLengthFromPoints = 0;
-
-Merge "FullInitialMeshFalcon.msh";
-
-CreateTopology;
-
-// Make all lines compound:
-ll[] = Line "*";
-For j In {0 : #ll[]-1}
-  Compound Line(1000+j) = ll[j];
-  Physical Line(1000+j) = (1000+j);
-EndFor
-
-// Make all surfaces compound and physical:
-ss[] = Surface "*";
-For i In {0 : #ss[]-1}
-  Compound Surface(i+1000) = ss[i];
-  Physical Surface(i+1) = { i+1000 };
-EndFor
-
-Mesh 2;
-Save "SurfaceMeshUniform.stl";
-Save Sprintf("SurfaceMeshUniform_%g.msh", lc);
-
-If(TWO_PLANES)
-Plugin(NewView).Run;
-x = -Pi/10;
-y = -Pi/10;
-z = -Pi/10;
-A = Cos(x); B = Sin(x); C = Cos(y); D = Sin(y); E = Cos(z); F = Sin(z);
-AD = A * D; BD = B * D;
-Plugin(Transform).A11 = C*E;
-Plugin(Transform).A12 = BD*E+A*F;
-Plugin(Transform).A13 =-AD*E+B*F;
-Plugin(Transform).A21 =-C*F;
-Plugin(Transform).A22 =-BD*F+A*E;
-Plugin(Transform).A23 = AD*F+B*E;
-Plugin(Transform).A31 = D;
-Plugin(Transform).A32 =-B*C;
-Plugin(Transform).A33 = A*C;
-Plugin(Transform).Tx = 11;
-Plugin(Transform).Ty = 1;
-Plugin(Transform).Tz = -1;
-Plugin(Transform).Run;
-Save "SurfaceMeshUniform2.stl";
-EndIf
-
-//Delete All;

benchmarks/3d_large/Falcon/fullfalcon_volume.script

-
-Include "param.geo";
-
-Merge "SurfaceMeshUniform.stl";
-If(TWO_PLANES)
-  Merge "SurfaceMeshUniform2.stl";
-EndIf
-
-zmin = General.MinZ-dd;
-zmax = General.MaxZ+dd;
-xmin = General.MinX-dd;
-xmax = General.MaxX+dd;
-ymin = General.MinY-dd;
-ymax = General.MaxY+dd;
-
-Point(20001) = {xmin, ymin, zmin, lc2};
-Point(20002) = {xmax, ymin, zmin, lc2};
-Point(20003) = {xmax, ymax, zmin, lc2};
-Point(20004) = {xmin, ymax, zmin, lc2};
-Point(20005) = {xmin, ymin, zmax, lc2};
-Point(20006) = {xmax, ymin, zmax, lc2};
-Point(20007) = {xmax, ymax, zmax, lc2};
-Point(20008) = {xmin, ymax, zmax, lc2};
-
-Line(10) = {20001, 20002};
-Line(11) = {20002, 20003};
-Line(12) = {20003, 20004};
-Line(13) = {20004, 20001};
-Line(14) = {20005, 20006};
-Line(15) = {20006, 20007};
-Line(16) = {20007, 20008};
-Line(17) = {20008, 20005};
-Line(18) = {20001, 20005};
-Line(19) = {20002, 20006};
-Line(20) = {20003, 20007};
-Line(21) = {20004, 20008};
-
-Line Loop(22) = {10,19,-14,-18};
-Line Loop(23) = {11,20,-15,-19};
-Line Loop(24) = {12,21,-16,-20};
-Line Loop(25) = {13,18,-17,-21};
-Line Loop(26) = {14,15,16,17};
-Line Loop(27) = {10,11,12,13};
-
-Plane Surface(22) = {22};
-Plane Surface(23) = {23};
-Plane Surface(24) = {24};
-Plane Surface(25) = {25};
-Plane Surface(26) = {26};
-Plane Surface(27) = {27};
-
-Surface Loop(20) = {22,23,24,25,26,27};
-Surface Loop(1) = {1};
-
-If(!TWO_PLANES)
-  Volume(1) = {20,1};
-  Physical Surface(1) = {1};
-EndIf
-If(TWO_PLANES)
-  Surface Loop(2) = {2};
-  Volume(1) = {20,1,2};
-  Physical Surface(1) = {1,2};
-EndIf
-
-Physical Surface(2) = {22:27};
-Physical Volume(4) = {1};
-
-Mesh 3;
-Save "falcon.msh";

benchmarks/3d_large/Falcon/fullsplit.script

-Include "param.geo";
-Merge "part.msh";
-N_DOM=NUM_SLICES;
-
-vol[] = {};         idx_vol[] = {};
-sur_scat[] = {};    idx_sur_scat[] = {};
-sur_inf[] = {};     idx_sur_inf[] = {};
-sigma_left[] = {};  idx_sigma_left[] = {};
-sigma_right[] = {}; idx_sigma_right[] = {};
-For i In {1:N_DOM}
-  idx_vol[] += #vol[];                 vol[] += Physical Volume{4000 + i};
-  idx_sur_scat[] += #sur_scat[];       sur_scat[] += Physical Surface{1000 + i};
-  idx_sur_inf[] += #sur_inf[];         sur_inf[] += Physical Surface{2000 + i};
-  idx_sigma_left[] += #sigma_left[];   sigma_left[] += Physical Surface{3000 + 1000 * (i-2)};
-  idx_sigma_right[] += #sigma_right[]; sigma_right[] += Physical Surface{3000 + 1000 * (i-1)};
-EndFor
-idx_vol[] += #vol[];
-idx_sur_scat[] += #sur_scat[];
-idx_sur_inf[] += #sur_inf[];
-idx_sigma_left[] += #sigma_left[];
-idx_sigma_right[] += #sigma_right[];
-
-Mesh.Binary = 1;
-Geometry.OrientedPhysicals = 0; // partitions have negative tags!!
-
-For i In {1:N_DOM}
-  idom = i-1;
-  Delete Physicals;
-  Physical Volume(4000+i) = vol[{idx_vol[i-1]:idx_vol[i]-1}];
-  Physical Surface(1000+i) = sur_scat[{idx_sur_scat[i-1]:idx_sur_scat[i]-1:1}];
-  Physical Surface(2000+i) = sur_inf[{idx_sur_inf[i-1]:idx_sur_inf[i]-1:1}];
-  Physical Surface(-(4000 + 1000 * (i-2))) = sigma_left[{idx_sigma_left[i-1]:idx_sigma_left[i]-1:1}];
-  Physical Surface(4000 + 1000 * (i-1)) = sigma_right[{idx_sigma_right[i-1]:idx_sigma_right[i]-1:1}];
-  Save Sprintf("falcon_mshcut%g.msh", idom);
-EndFor
-

benchmarks/3d_large/Falcon/param.geo

-lc = 0.25; // surface
-lc2 = lc * 5; // volume
-dd = 2;
-TWO_PLANES = 1;
-NUM_SLICES = 4;

benchmarks/3d_large/Falcon/partition.script

-Include "param.geo";
-
-Merge "falcon.msh";
-
-Plugin(SimplePartition).NumSlices = NUM_SLICES;
-Plugin(SimplePartition).Run;
-
-Mesh.Binary=1;
-Mesh.MshFilePartitioned = 2; // hack to save a physical for each partition
-
-Save "part.msh";

benchmarks/3d_large/Falcon/split.script

-Include "param.geo";
-Merge "part.msh";
-N_DOM=NUM_SLICES;
-
-vol[] = {};         idx_vol[] = {};
-sur_scat[] = {};    idx_sur_scat[] = {};
-sur_inf[] = {};     idx_sur_inf[] = {};
-sur_sym[] = {};     idx_sur_sym[] = {};
-sigma_left[] = {};  idx_sigma_left[] = {};
-sigma_right[] = {}; idx_sigma_right[] = {};
-For i In {1:N_DOM}
-  idx_vol[] += #vol[];                 vol[] += Physical Volume{4000 + i};
-  idx_sur_scat[] += #sur_scat[];       sur_scat[] += Physical Surface{1000 + i};
-  idx_sur_inf[] += #sur_inf[];         sur_inf[] += Physical Surface{2000 + i};
-  idx_sur_sym[] += #sur_sym[];         sur_sym[] += Physical Surface{3000 + i};
-  idx_sigma_left[] += #sigma_left[];   sigma_left[] += Physical Surface{4000 + 1000 * (i-2)};
-  idx_sigma_right[] += #sigma_right[]; sigma_right[] += Physical Surface{4000 + 1000 * (i-1)};
-EndFor
-idx_vol[] += #vol[];
-idx_sur_scat[] += #sur_scat[];
-idx_sur_inf[] += #sur_inf[];
-idx_sur_sym[] += #sur_sym[];
-idx_sigma_left[] += #sigma_left[];
-idx_sigma_right[] += #sigma_right[];
-
-Mesh.Binary = 1;
-Geometry.OrientedPhysicals = 0; // partitions have negative tags!!
-
-For i In {1:N_DOM}
-  idom = i-1;
-  Delete Physicals;
-  Physical Volume(4000+i) = vol[{idx_vol[i-1]:idx_vol[i]-1}];
-  Physical Surface(1000+i) = sur_scat[{idx_sur_scat[i-1]:idx_sur_scat[i]-1:1}];
-  Physical Surface(2000+i) = sur_inf[{idx_sur_inf[i-1]:idx_sur_inf[i]-1:1}];
-  Physical Surface(3000+i) = sur_sym[{idx_sur_sym[i-1]:idx_sur_sym[i]-1:1}];
-  Physical Surface(-(4000 + 1000 * (i-2))) = sigma_left[{idx_sigma_left[i-1]:idx_sigma_left[i]-1:1}];
-  Physical Surface(4000 + 1000 * (i-1)) = sigma_right[{idx_sigma_right[i-1]:idx_sigma_right[i]-1:1}];
-  Save Sprintf("falcon_mshcut%g.msh", idom);
-EndFor
-