diff --git a/benchmarks/bugs/projection.geo b/benchmarks/bugs/projection.geo new file mode 100644 index 0000000000000000000000000000000000000000..ca302181108ba14784307ca748320fe1a9d0479d --- /dev/null +++ b/benchmarks/bugs/projection.geo @@ -0,0 +1,13 @@ +// +Point(1) = {0,0,0,1}; +Point(2) = {0,10,0,1}; +Line(1) = {2,1}; +Extrude Line {1, {50.0,0.0,0.0}, {1,0,0}, {0,5,0}, Pi}{Recombine;Layers{50,9000,1};}; +// +Point(1000) = {0,0,20,1}; +Point(2000) = {0,10,20,1}; +Line(1000) = {2000,1000}; +Extrude Line {1000, {50.0,0.0,0.0}, {1,0,0}, {0,5,20}, Pi*0.45}; + + + diff --git a/benchmarks/bugs/torus.geo b/benchmarks/bugs/torus.geo new file mode 100644 index 0000000000000000000000000000000000000000..6c275f7cac6e173c21f758b8713973aea696f687 --- /dev/null +++ b/benchmarks/bugs/torus.geo @@ -0,0 +1,49 @@ +/********************************************************************* + * + * Torus.geo Gmsh Program to Construct a Torus + * + *********************************************************************/ + +cm = 1e-02; // Define our scale +Lc1 = 0.5 * cm; // Mesh size +r1 = 4 * cm; // Origin->Disk center radius +r2 = 1 * cm; // Disk radius + +r2c = r2 * Cos( Pi/4 ); // Cos/Sine terms to define disk circumference points +r2s = r2 * Sin( Pi/4 ); + +deg90 = Pi / 2; // Radian representation of 90 Degs + +Point(1) = { 0.0, 0.0, 0.0, Lc1 }; // Origin + +Point(2) = { 0.0, r1 , 0.0 , Lc1 }; // Circle center +Point(3) = { 0.0, r1+r2 , 0.0 , Lc1 }; // 8 points along the circumference +Point(4) = { 0.0, r1+r2c, r2s , Lc1 }; +Point(5) = { 0.0, r1 , r2 , Lc1 }; +Point(6) = { 0.0, r1-r2c, r2s , Lc1 }; +Point(7) = { 0.0, r1-r2 , 0.0 , Lc1 }; +Point(8) = { 0.0, r1-r2c, -r2s, Lc1 }; +Point(9) = { 0.0, r1 , -r2 , Lc1 }; +Point(10) = { 0.0, r1+r2c, -r2s, Lc1 }; + +Circle(1) = { 3, 2, 4 } ; // Arcsegments along the circumference +Circle(2) = { 4, 2, 5 } ; +Circle(3) = { 5, 2, 6 } ; +Circle(4) = { 6, 2, 7 } ; +Circle(5) = { 7, 2, 8 } ; +Circle(6) = { 8, 2, 9 } ; +Circle(7) = { 9, 2, 10 } ; +Circle(8) = { 10, 2, 3 } ; + +// Merge arcsegments into a disk +Line Loop (9) = { 1,2,3,4,5,6,7,8 }; +Plane Surface(10) = {9}; + +// Extrude the disk to form a partial torus + +// Things are fine as long as extrusion angle is <= Pi +/* Extrude Surface {10, {0,0,1}, {0,0,0}, 1* deg90 }; */ + +// Things get somewhat odd for angles > Pi +Extrude Surface {10, {0,0,1}, {0,0,0}, 1.9* deg90 };// { Layers{100,10,1}; }; +