couple of test files

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_aniso.geo

+Mesh.CharacteristicLengthFactor = 1;
+Mesh.Algorithm = 2; /*2D mesh algorithm (1=isotropic, 2=anisotropic, 3=triangle)*/
+
+
+CLF = 0.995 ; // MARCHE PAS
+//CLF = 0.999 ; //OK
+
+RRes = 0.001;  dzRes1 = 0.01 ; 
+REl1 = 0.025;  dzEl1 = REl1/25 ;  
+RLay1 = 0.030;  dzLay1 = 0.0003 ;
+RLay2 = 0.050;  dzLay2 = 0.050 ; 
+REl2 = 0.025; dzEl2 = REl2/25 ; 
+
+
+SigCopper = 1/1.0e-8 ;
+SigGel = 1/5000.0 ;
+SigMuscle = 1/5.0 ;
+Res = 1e6 ; 
+SigRes = 1/Res * dzRes1 / Pi / RRes^2 ;
+
+j = 0.1/Pi/0.025^2 ;
+I = j * Pi * REl1^2 ; VElec1Pot = I * Res ;
+
+
+
+
+dztot = 0 ;
+
+
+/* Resistance */
+Point(1) = {0,dztot,0.0,RRes/2*CLF};
+Point(2) = {RRes,dztot,0,RRes/2*CLF};
+dztot = dztot - dzRes1 ; 
+Point(3) = {RRes,dztot,0.0,RRes/2*CLF};
+Point(4) = {0,dztot,0.0,RRes/2*CLF};		// not used
+
+/* electrode 1 */
+lc1 = REl1/50 ; lc2 = REl1/50 ;
+Point(5) = {0,dztot,0.0,lc1*CLF};
+Point(6) = {REl1,dztot,0,lc1*CLF};
+dztot = dztot - dzEl1 ;
+Point(7) = {REl1,dztot,0.0,lc2/10*CLF};
+Point(8) = {0,dztot,0.0,lc2*CLF};		// not used
+
+
+/* layer 1 Gel or Resistive slab */ 
+lc1 = RLay1/250 ; lc2 = RLay1/250 ;
+Point(9) = {0,dztot,0.0,REl1/100*CLF};
+Point(10) = {RLay1,dztot,0,lc2*CLF};
+dztot = dztot - dzLay1 ; dzElecSurf1 = dztot ;
+Point(11) = {RLay1,dztot,0,lc2*CLF};
+/*Point(12) = {0,dztot,0.0,lc1*CLF};		// not used
+
+/* layer 2 muscle */
+lc1 = REl1/10*CLF; lc2 = Fabs(REl1-RLay2)/10*CLF ; lc3 = dzLay2/10*CLF ;
+
+Point(13) = {0.0,dztot,0.0,lc1/20};
+Point(14) = {REl1,dztot,0.0,lc1/20};	/*attractor*/
+Point(15) = {RLay2,dztot,0,lc2/2};
+dztot = dztot - dzLay2/5 ;			/*attractor*/
+Point(16) = {RLay2,dztot,0.0,Sqrt(lc2 * lc3)/2};
+Point(17) = {0,dztot,0.0,Sqrt(lc1 * lc3)/2};
+dztot = dztot - 3*dzLay2/5 ;			/*attractor*/
+Point(18) = {RLay2,dztot,0.0,Sqrt(lc3* lc2)/2 };
+Point(19) = {0,dztot,0.0,Sqrt(lc3* lc1)/2};
+dztot = dztot - dzLay2/5 ; dzElecSurf2 = dztot ;
+Point(20) = {RLay2,dztot,0,lc2/2};	
+/*Point(21) = {0,dztot,0.0,lc1/20};		// not used
+
+/* electrode 2 */
+lc1 = REl2/100 ; 
+Point(22) = {0,dztot,0.0,lc1*CLF};
+Point(23) = {REl2,dztot,0,lc1/5*CLF};
+dztot = dztot - dzEl2 ;
+Point(24) = {REl2,dztot,0.0,1};
+Point(25) = {0,dztot,0.0,1};
+
+/* ElectrodeCenterEdge */
+Point(26) = {0,dzElecSurf1,0.0,1};
+Point(27) = {REl1,dzElecSurf1,0.0,1};
+Point(28) = {0,dzElecSurf2,0.0,1};
+Point(29) = {REl2,dzElecSurf2,0.0,1};
+
+
+Line(1) = {1,2};
+Line(2) = {2,3};
+Line(3) = {3,5};
+Line(4) = {5,1};
+Line(5) = {3,6};
+Line(6) = {6,7};
+Line(7) = {7,9};
+Line(8) = {9,5};
+Line(9) = {7,10};
+Line(10) = {10,11};
+Line(11) = {11,14};
+Line(12) = {14,13};
+Line(13) = {13,9};
+Line(14) = {11,15};
+Line(15) = {15,16};
+Line(16) = {16,18};
+Line(17) = {18,20};
+Line(18) = {20,23};
+Line(19) = {23,22};
+Line(20) = {22,19};
+Line(21) = {19,17};
+Line(22) = {17,13};
+Line(23) = {23,24};
+Line(24) = {24,25};
+Line(25) = {25,22};
+
+Line Loop(26) = {2,3,4,1};
+Plane Surface(27) = {26};
+Line Loop(28) = {5,6,7,8,-3};
+Plane Surface(29) = {28};
+Line Loop(30) = {9,10,11,12,13,-7};
+Plane Surface(31) = {30};
+Line Loop(32) = {15,16,17,18,19,20,21,22,-12,-11,14};
+Plane Surface(33) = {32};
+Physical Point(100) = {26,27,28,29};
+Physical Line(101) = {1};
+Physical Line(102) = {19};
+Physical Surface(1001) = {27};
+Physical Surface(1002) = {29};
+Physical Surface(1003) = {31};
+Physical Surface(1004) = {33};
+