This commit was manufactured by cvs2svn to create tag 'gmsh_1_62'.

benchmarks/bugs/bug_1element.geo

+maille=100;
+Point(1)={0,0,24.5,maille};
+Point(2)={0,0,-24.5,maille};
+Point(3)={0,24.5,0,maille};
+Point(4)={0,-24.5,0,maille};
+Point(5)={0,0,0,maille};
+Point(6)={100,-42.5,17.5,maille};
+Point(7)={100,42.5,17.5,maille};
+Point(8)={100,42.5,-17.5,maille};
+Point(9)={100,-42.5,-17.5,maille};
+Circle(1) = {1,5,3};
+Circle(2) = {3,5,2};
+Circle(3) = {2,5,4};
+Circle(4) = {4,5,1};
+Line(5) = {6,7};
+Line(6) = {7,8};
+Line(7) = {8,9};
+Line(8) = {9,6};
+Line(9) = {1,7};
+Line(10) = {1,6};
+Line(11) = {3,7};
+Line(12) = {3,8};
+Line(13) = {2,8};
+Line(14) = {2,9};
+Line(15) = {4,9};
+Line(16) = {4,6};
+Line Loop(17) = {5,-9,10};
+Plane Surface(18) = {17};
+Line Loop(19) = {11,6,-12};
+Plane Surface(20) = {19};
+Line Loop(21) = {13,7,-14};
+Plane Surface(22) = {21};
+Line Loop(23) = {16,-8,-15};
+Plane Surface(24) = {23};
+Line Loop(25) = {7,8,5,6};
+Line Loop(27) = {1,2,3,4};
+Plane Surface(28) = {27};
+Line Loop(29) = {9,-11,-1};
+Ruled Surface(30) = {29};
+Line Loop(31) = {12,-13,-2};
+Ruled Surface(32) = {31};
+Line Loop(33) = {14,-15,-3};
+Ruled Surface(34) = {33};
+Line Loop(35) = {16,-10,-4};
+Ruled Surface(36) = {35};
+Point(37)={420,-95,17.5,maille};
+Point(38)={420,95,17.5,maille};
+Point(39)={420,95,-17.5,maille};
+Point(40)={420,-95,-17.5,maille};
+Line(37) = {37,38};
+Line(38) = {38,39};
+Line(39) = {39,40};
+Line(40) = {40,37};
+Line(41) = {6,37};
+Line(42) = {7,38};
+Line(43) = {8,39};
+Line(44) = {9,40};
+Line Loop(45) = {37,-42,-5,41};
+Plane Surface(46) = {45};
+Line Loop(47) = {41,-40,-44,8};
+Plane Surface(48) = {47};
+Line Loop(49) = {42,38,-43,-6};
+Plane Surface(50) = {49};
+Line Loop(51) = {44,-39,-43,7};
+Plane Surface(52) = {51};
+Line Loop(53) = {37,38,39,40};
+Point(55)={720,-150,110,5};
+Point(56)={720,150,110,5};
+Point(57)={720,150,-110,5};
+Point(58)={720,-150,-110,5};
+Line(55) = {55,56};
+Line(56) = {56,57};
+Line(57) = {57,58};
+Line(58) = {58,55};
+Line(59) = {37,55};
+Line(60) = {38,56};
+Line(61) = {39,57};
+Line(62) = {40,58};
+Line Loop(63) = {59,-58,-62,40};
+Plane Surface(64) = {63};
+Line Loop(65) = {59,55,-60,-37};
+Plane Surface(66) = {65};
+Line Loop(67) = {62,-57,-61,39};
+Plane Surface(68) = {67};
+Line Loop(69) = {60,56,-61,-38};
+Plane Surface(70) = {69};
+Line Loop(71) = {55,56,57,58};
+Point(73)={880,-210,150,maille};
+Point(74)={880,210,150,maille};
+Point(75)={880,210,-150,maille};
+Point(76)={880,-210,-150,maille};
+Line(73) = {73,74};
+Line(74) = {74,75};
+Line(75) = {75,76};
+Line(76) = {76,73};
+Line(77) = {55,73};
+Line(78) = {56,74};
+Line(79) = {57,75};
+Line(80) = {58,76};
+Line Loop(81) = {77,73,-78,-55};
+Plane Surface(82) = {81};
+Line Loop(83) = {78,74,-79,-56};
+Plane Surface(84) = {83};
+Line Loop(85) = {79,75,-80,-57};
+Plane Surface(86) = {85};
+Line Loop(87) = {80,76,-77,-58};
+Plane Surface(88) = {87};
+Line Loop(89) = {73,74,75,76};
+Point(91)={1220,-330,290,maille};
+Point(92)={1220,330,290,maille};
+Point(93)={1220,330,-290,maille};
+Point(94)={1220,-330,-290,maille};
+Line(91) = {91,92};
+Line(92) = {92,93};
+Line(93) = {93,94};
+Line(94) = {94,91};
+Line(95) = {73,91};
+Line(96) = {76,94};
+Line(97) = {75,93};
+Line(98) = {74,92};
+Line Loop(99) = {95,91,-98,-73};
+Plane Surface(100) = {99};
+Line Loop(101) = {95,-94,-96,76};
+Plane Surface(102) = {101};
+Line Loop(103) = {96,-93,-97,75};
+Plane Surface(104) = {103};
+Line Loop(105) = {98,92,-97,-74};
+Plane Surface(106) = {105};
+Line Loop(107) = {91,92,93,94};
+Point(109)={1360,-390,380,maille};
+Point(110)={1360,390,380,maille};
+Point(111)={1360,390,-380,maille};
+Point(112)={1360,-390,-380,maille};
+Line(109) = {109,110};
+Line(110) = {110,111};
+Line(111) = {111,112};
+Line(112) = {112,109};
+Line(113) = {91,109};
+Line(114) = {92,110};
+Line(115) = {93,111};
+Line(116) = {94,112};
+Line Loop(117) = {109,-114,-91,113};
+Plane Surface(118) = {117};
+Line Loop(119) = {113,-112,-116,94};
+Plane Surface(120) = {119};
+Line Loop(121) = {116,-111,-115,93};
+Plane Surface(122) = {121};
+Line Loop(123) = {114,110,-115,-92};
+Plane Surface(124) = {123};
+Line Loop(125) = {110,111,112,109};
+Plane Surface(126) = {125};
+Surface Loop(127) = {126,-124,-118,-100,102,120,122,104,-106,-84,-82,-88,-86,68,64,-66,-70,-50,-46,-18,-30,-20,-32,-22,52,48,24,-36,-28,-34};
+Volume(128) = {127};

benchmarks/bugs/bug_demi_cercles.geo

+Point(1) = {0.0,0.0,0.0,.08};
+Point(2) = {-.5,0.0,0.0,.08};
+Point(3) = {.5,0.0,0.0,.08};
+Point(4) = {.5,-.5,0.0,.08};
+Point(5) = {-.5,-.5,0.0,.08};
+Circle(1) = {3,1,2} Plane {0,0,1};
+Line(2) = {3,4};
+Line(3) = {4,5};
+Line(4) = {5,2};
+Line Loop(5) = {3,4,-1,2};
+Plane Surface(6) = {5};
+Extrude Surface {6, {0,0,.5}};
+Coherence;

benchmarks/bugs/bug_ruled.geo

+
+// solved by setting the tolerance lower in sys3x3_with_tol (-> again,
+// this arised from a wrong mean plane computation)
+
+//lc = 0.00001;
+lc = 0.001;
+Point(61) = {0.058, -0.005, 0, lc};
+Point(62) = {0.058, -0.005, 0.000625, lc};
+Point(64) = {0.058625, -0.005, 0, lc};
+
+Point(85) = {0.058, -0.006, 0, lc};
+Point(86) = {0.058, -0.006, 0.000625, lc};
+Point(88) = {0.058625, -0.006, 0, lc};
+Line(1) = {86,62};
+Line(2) = {88,64};
+Circle(3) = {62,61,64};
+Circle(4) = {86,85,88};
+Line Loop(5) = {2,-3,-1,4};
+Ruled Surface(6) = {5};

benchmarks/bugs/disk.geo

+
+// proj. problem with angle=pi/2
+
+angle=Pi/2;
+
+mpi     = 0.0254;
+
+dia1    = 3.50 * mpi;
+rad1    = dia1 / 2.0;
+
+dia2    = 75.00 * mpi;
+rad2    = dia2 / 2.0;
+
+thick   = 1.0 * mpi;
+length1 = 1.5 * dia1;
+
+x0      = 0.0;
+y0      = 0.0;
+z0      = 0.0;
+
+y1  = y0 - 0.50 * length1 ;
+yr1 = y1 - 0.25 * length1 ;
+
+y2  = y1 + length1;
+yr2 = y2 + 0.25 * length1;
+
+lcart1  = 0.05 * rad1;
+// lcart2  = 0.100 * rad2;
+// lcart2  = 0.050 * rad2;
+lcart2  = 0.025 * rad2;
+
+Printf ( " %f %f ", lcart1, lcart2 );
+
+p01 = newp; Point(p01) = { x0 + rad1,  yr1,   z0, lcart1};
+p02 = newp; Point(p02) = { x0 + rad2,  yr1,   z0, lcart2};
+p03 = newp; Point(p03) = { x0 + rad2,  yr2,   z0, lcart2};
+p04 = newp; Point(p04) = { x0 + rad1,  yr2,   z0, lcart1};
+
+l01  = newreg; Line (l01) = { p01, p02 };
+l02  = newreg; Line (l02) = { p02, p03 };
+l03  = newreg; Line (l03) = { p03, p04 };
+l04  = newreg; Line (l04) = { p01, p04 };
+
+a1 = newreg; s1 = newreg; Line Loop ( a1 ) = { l01, l02, l03, -l04 } ;  Plane Surface ( s1 ) = { a1 };
+s11[] = Extrude Surface { s1, { 0, 1, 0 }, { 0.0, 0.0, 0.0 }, -angle  }; ;

benchmarks/bugs/nested_tramsforms.geo

+a = 0.004;
+b = 0.002;
+lc = 0.001;
+Point(1) = { 0, 0, 0, lc };
+Point(100) = { a, 0, 0, lc };
+Point(101) = { 0, b, 0, lc };
+Point(102) = {-a, 0, 0, lc };
+Point(103) = { 0,-b, 0, lc };
+Ellipse(100) = { 100, 1, 100, 101 };
+Ellipse(101) = { 101, 1, 102, 102 };
+Ellipse(102) = { 102, 1, 102, 103 };
+Ellipse(103) = { 103, 1, 100, 100 };
+Line Loop(4) = { 100, 101, 102, 103 };
+Plane Surface(4) = { 4 };
+
+
+
+to_rad = 3.141592654/180;
+
+Translate { -0.01, 0.01, 0 } 
+{ 
+
+    Rotate { {0, 0, 1}, {0, 0, 0}, 120*to_rad } 
+    { 
+      Duplicata { Surface {4}; } 
+    }
+
+}

benchmarks/bugs/p19-bug.geo

+
+Mesh.CharacteristicLengthFactor = 2;
+
+fact     = 0.4 ;
+rondelle = fact * 0.01;
+iris     = fact * 0.004;
+size     = fact * 0.01;
+
+larg = 86.36e-3 / 2.0 ;
+long = 45.0e-3 ;
+l    = 16.e-3 / 2.0 ;
+re   = 90.e-3 / 2.0 ;
+re2  = re + 3.e-3 ;
+ri   = 8.e-3 / 2.0 ;
+ll   = 60.0e-3 / 2.0;
+a    = larg - ll;
+c45  = 0.5^0.5  ;
+
+hg   = 21.59e-3 ;
+hcav = 42.5e-3 ;
+
+x2   = long;
+y1   = ri * c45;
+y2   = ri ;
+y3   = l ;
+ss1  = l / re ;
+xtemp = re * ( ( 1.0 - ss1 * ss1 ) ^ 0.5 ) ;
+ss2  = ( 1.0 - ( xtemp / re2 ) * ( xtemp / re2 ) ) ^ 0.5 ;
+y4   = ss2 * xtemp ;
+ss3  = ll / re2;
+xtemp2 = re2 *  ( ( 1.0 - ss3 * ss3 ) ^ 0.5 ) ;
+x1   = long - ( xtemp - xtemp2 ) ;
+x6   = x1 + xtemp ;
+x3   = x6 - re * c45 ;
+x4   = x6 - ri ;
+x5   = x6 - ri * c45 ;
+x7   = x6 + ri * c45 ;
+x8   = x6 + ri ;
+x9   = x6 + re * c45 ;
+x10  = x6 + re ;
+y5   = ll ;
+y6   = re * c45 ;
+y7   = re ;
+y8   = larg ;
+
+
+Point(1) = {0.0,0.0,0.0,size};
+Point(2) = {0.0,y8,0.0,size};
+Point(3) = {x2,y8,0.0,size};
+Point(4) = {x2,y5,0.0,size};
+Point(5) = {x1,y4,0.0,size};
+Point(6) = {x1,y3,0.0,iris};
+Point(7) = {x1,0.0,0.0,iris};
+Point(8) = {x3,y6,0.0,size};
+Point(9) = {x6,y7,0.0,size};
+Point(10) = {x9,y6,0.0,size};
+Point(11) = {x10,0.0,0.0,size};
+Point(12) = {x8,0.0,0.0,rondelle};
+Point(13) = {x6,0.0,0.0,rondelle};
+Point(14) = {x4,0.0,0.0,rondelle};
+Point(15) = {x7,y1,0.0,rondelle};
+Point(16) = {x6,y2,0.0,rondelle};
+Point(17) = {x5,y1,0.0,rondelle};
+
+Line(1) = {1,2};
+Line(2) = {2,3};
+Line(3) = {3,4};
+Line(4) = {5,6};
+Line(5) = {1,7};
+Line(6) = {7,14};
+Line(10) = {6,7};
+Line(11) = {14,13};
+Line(12) = {13,12};
+Circle(16) = {11,13,10};
+Circle(17) = {10,13,9};
+Circle(18) = {9,13,8};
+Circle(19) = {8,13,6};
+Circle(20) = {4,13,5};
+Line(21) = {11,12};
+Circle(22) = {17,13,14};
+Circle(23) = {16,13,17};
+Circle(24) = {15,13,16};
+Circle(25) = {12,13,15};
+
+Line Loop(26) = {-5,1,2,3,20,4,10};
+Plane Surface(27) = {26};
+
+Line Loop(28) = {-6,-10,-19,-18,-17,-16,21,25,24,23,22};
+Plane Surface(29) = {28};
+
+Line Loop(30) = {11,12,25,24,23,22};
+Plane Surface(31) = {30};
+
+Extrude Surface{27, {0,0,hg} };
+Coherence;
+
+Extrude Surface{29, {0,0,hg} };
+Coherence;
+
+Extrude Surface{31, {0,0,hg} };
+Coherence;
+
+Extrude Surface {105, {0,0,hcav-hg} };
+Coherence;
+
+Extrude Surface {126, {0,0,hcav-hg} };
+Coherence;
+
+//Characteristic Length {58,71} = 0.01; 

benchmarks/bugs/piece-bad.geo

+/* 
+   Gmsh demo file (C) 2000 C. Geuzaine, J.-F. Remacle
+*/
+
+r1 = .1;
+l1 = 1.;
+l2 = .8;
+l3 = .1;
+r2 = 1.1;
+lc = .05;
+lc2 = .05;
+rint = .2;
+rext = .3;
+
+Point(1) = {0.0,0.0,0.0,lc};
+Point(2) = {l1,0.0,0.0,lc2};
+Point(3) = {l1-r1,0.0,0.0,lc2};
+Point(4) = {l1,r1,0.0,lc2};
+Point(5) = {l1,-r1,0.0,lc2};
+Point(6) = {l1+l2,r1,0.0,lc};
+Point(7) = {l1+l2,-r1,0.0,lc};
+Point(8) = {l1+l2,-r1-l3,0.0,lc};
+Point(9) = {l1+l2,r1+l3,0.0,lc};
+
+Line(1) = {4,6};
+Line(2) = {6,9};
+Line(3) = {7,8};
+Line(4) = {5,7};
+Circle(5) = {4,2,3};
+Circle(6) = {3,2,5};
+
+r = 2*3.14159/5;
+Point(10) = { (l1 + r2) * Cos(r/2) , (l1 + r2) * Sin(r/2), 0.0, lc};
+
+// Remember, all rotations are specified by the axis direction
+// ({0,0,1}), an axis point ({0,0,0}) and a rotation angle (r)
+
+i = 0 ;
+
+For(1:4)
+
+  i+=1 ; 
+
+  Rotate {{0.0,0.0,1.0},{0.0,0.0,0.0}, i*r} {
+    Duplicata {
+      Line{1}; Line{2}; Line{3}; 
+      Line{4}; Line{5}; Line{6}; 
+      Point{10};
+    }
+  }
+
+EndFor
+
+Point(newp) = {rint,0,0,lc};
+Point(newp) = {rext,0,0,lc};
+Point(newp) = {-rint,0,0,lc};
+Point(newp) = {-rext,0,0,lc};
+Point(newp) = {0,rint,0,lc};
+Point(newp) = {0,rext,0,lc};
+Point(newp) = {0,-rint,0,lc};
+Point(newp) = {0,-rext,0,lc};
+
+Circle(31) = {8,118,97};
+Circle(32) = {20,10,9};
+Circle(33) = {47,37,16};
+Circle(34) = {74,64,43};
+Circle(35) = {101,91,70};
+Circle(36) = {119,1,123};
+Circle(37) = {123,1,121};
+Circle(38) = {121,1,125};
+Circle(39) = {125,1,119};
+Circle(40) = {124,1,122};
+Circle(41) = {122,1,126};
+Circle(42) = {126,1,120};
+Circle(43) = {120,1,124};
+
+Line Loop(44) = {37,38,39,36};
+Line Loop(46) = {40,41,42,43};
+Plane Surface(47) = {46,44};
+
+Line Loop(48) = {-2,-1,5,6,4,3,31,-26,-25,29,30,28,27,35,-20,-19,
+  23,24,22,21,34,-14,-13,17,18,16,15,33,-8,-7,11,12,10,9,32};
+Plane Surface(49) = {48,46};
+
+Extrude Surface {47, {0,0,0.2}};
+Extrude Surface {49, {0,0,0.2}};
+Extrude Surface {91, {0,0,0.2}};
+Extrude Surface {47, {0,0,-0.2}};
+
+Physical Volume(1) = {1,2,3,4};
+Physical Surface(2) = {1:1000};

benchmarks/bugs/transfo_neg_curves.geo

+Point(1) = {2,5,0,0.1};
+Point(2) = {5,5,0,0.1};
+Point(3) = {2,10,0,0.1};
+Point(4) = {2,0,0,0.1};
+Point(5) = {7,5,0,0.1};
+Circle(3) = {4,1,5};
+Circle(4) = {5,1,3};
+
+// this used to fail (due to the neg curves in the line loop) in
+// gmsh <= 1.61
+
+Ellipse(1) = {4,1,1,2};
+Ellipse(2) = {2,1,1,3};
+Line Loop(5) = {3,4,-2,-1};
+Plane Surface(6) = {5};
+Symmetry { 1.0,0.0,0.0,-2 } {  Surface{6}; }
+//Rotate { {0.0,0,1}, {2,5,0.0}, Pi/7 } { Surface{6}; }
+
+
+
+// this worked:
+/*
+Ellipse(1) = {2,1,1,4};
+Ellipse(2) = {3,1,1,2};
+Line Loop(5) = {3,4,2,1};
+Plane Surface(6) = {5};
+Symmetry { 1.0,0.0,0.0,-2 } {  Surface{6}; }
+//Rotate { {0.0,0,1}, {2,5,0.0}, Pi/7 } { Surface{6}; }
+
+*/

benchmarks/bugs/xytouv2.geo

+// A shell
+rc=0.025;
+Point(1) = {5.908, 2.821, 0, rc};
+Point(2) = {5.936, 2.661, 0, rc};
+Point(3) = {5.956, 2.479, 0, rc};
+Point(4) = {5.986, 1.869, 0, rc};
+Point(5) = {5.996, 0.744, 0, rc};
+Point(6) = {5.856, 2.611, 1.175, rc};
+Point(7) = {5.876, 2.439, 1.191, rc};
+Point(8) = {5.915, 1.798, 1.238, rc};
+Point(9) = {5.903, 2.816, 0.291, rc};
+Point(10) = {5.891, 2.801, 0.583, rc};
+Point(11) = {5.869, 2.776, 0.874, rc};
+Point(12) = {5.838, 2.741, 1.166, rc};
+Point(13) = {5.992, 0.744, 0.321, rc};
+Point(14) = {5.979, 0.744, 0.641, rc};
+Point(15) = {5.957, 0.744, 0.962, rc};
+Point(16) = {5.926, 0.744, 1.282, rc};
+
+Spline(1) = {1,2,3,4,5};
+Spline(2) = {5,13,14,15,16};
+Spline(3) = {16,8,7,6,12};
+Spline(4) = {12,11,10,9,1};
+Line Loop(1)={1,2,3,4};
+Ruled Surface(1) = {1};
+
+Physical Surface(1) = {1};

examples/cube.geo

-lcar1 = .0275;
-
-length = 0.25;
-height = 1.0;
-depth = 0.25;
-
-Point(newp) = {length/2,height/2,depth,lcar1}; /* Point      1 */
-Point(newp) = {length/2,height/2,0,lcar1}; /* Point      2 */
-Point(newp) = {-length/2,height/2,depth,lcar1}; /* Point      3 */
-Point(newp) = {-length/2,-height/2,depth,lcar1}; /* Point      4 */
-Point(newp) = {length/2,-height/2,depth,lcar1}; /* Point      5 */
-Point(newp) = {length/2,-height/2,0,lcar1}; /* Point      6 */
-Point(newp) = {-length/2,height/2,0,lcar1}; /* Point      7 */
-Point(newp) = {-length/2,-height/2,0,lcar1}; /* Point      8 */
-Line(1) = {3,1};
-Line(2) = {3,7};
-Line(3) = {7,2};
-Line(4) = {2,1};
-Line(5) = {1,5};
-Line(6) = {5,4};
-Line(7) = {4,8};
-Line(8) = {8,6};
-Line(9) = {6,5};
-Line(10) = {6,2};
-Line(11) = {3,4};
-Line(12) = {8,7};
-Line Loop(13) = {-6,-5,-1,11};
-Plane Surface(14) = {13};
-Line Loop(15) = {4,5,-9,10};
-Plane Surface(16) = {15};
-Line Loop(17) = {-3,-12,8,10};
-Plane Surface(18) = {17};
-Line Loop(19) = {7,12,-2,11};
-Plane Surface(20) = {19};
-Line Loop(21) = {-4,-3,-2,1};
-Plane Surface(22) = {21};
-Line Loop(23) = {8,9,6,7};
-Plane Surface(24) = {23};
-
-Surface Loop(25) = {14,24,-18,22,16,-20};
-Complex Volume(26) = {25};
-
-n = 21;
-
-/*
-sizex = 4;
-sizey = 6;
-sizez = 4;
-*/
-sizex = n*length;
-sizey = n*height;
-sizez = n*depth;
-
-sizex = 9;
-sizey = 33;
-sizez = 9;
-
-
-Transfinite Line {2,4,9,7} = sizez;
-Transfinite Line {11,5,10,12} = sizey;
-Transfinite Line {3,1,6,8} = sizex;
-
-Transfinite Surface {14} = {5,4,3,1};
-Transfinite Surface {16} = {6,2,1,5};
-Transfinite Surface {18} = {6,2,7,8};
-Transfinite Surface {20} = {3,7,8,4};
-Transfinite Surface {22} = {3,7,2,1};
-Transfinite Surface {24} = {4,8,6,5};
-
-Recombine Surface {14,16,18,20,22,24};
-
-Transfinite Volume {26} = {3,7,2,1,4,8,6,5};
-
-Physical Surface(200) = {14,16,18,20,24};
-
-Physical Volume(100) = {26};
-/*
-Mesh.use_cut_plane = 1;
-Mesh.cut_planea = 0;
-Mesh.cut_planeb = 0;
-Mesh.cut_planec = 1;
-Mesh.cut_planed = -0.125;
-*/
\ No newline at end of file