[Gmsh] Quality of constrained Delaunay triangulations - Automatic removal of trisected elements

Wed Sep 28 13:59:59 CEST 2016

Hello,

For certain geometries GMSH generates triangulations containing interior 
vertices connected to only three neighbours. The three angles meeting at such a 
point are all obtuse. The three associated elements can always be merged to 
improve the mesh quality. It were useful if gmsh automatically improved the 
meshes accordingly, similar to what it does with convex quartered elements when 
supplied with the Remove4Triangles option. In general, it would be really nice 
to have an option that limits the connectivity of interior vertices to just 5, 6 
and 7.

Best regards,

--
Karl Kästner
------------
Phd Student at Wageningen University & Research centre
Hydrology and Quantitative Water Management Group
Droevendaalsesteeg 3, Lumen building, Room D101, 6708 PB Wageningen, The Netherlands
Phone: +31617876970 (Netherlands)
Mail: karl.kastner at wur.nl
-------------- next part --------------
Point(1) = {316149.640172,-50772.708746,0.000000,100.000000};
Point(2) = {316209.731626,-50851.961995,0.000000,100.000000};
Point(3) = {316265.748554,-50934.119710,0.000000,100.000000};
Point(4) = {316322.288007,-51015.909604,0.000000,100.000000};
Point(5) = {316380.499920,-51096.548525,0.000000,100.000000};
Point(6) = {316439.931681,-51176.300974,0.000000,100.000000};
Point(7) = {316496.596275,-51258.040506,0.000000,100.000000};
Point(8) = {316552.890947,-51340.028981,0.000000,100.000000};
Point(9) = {316611.252665,-51420.545412,0.000000,100.000000};
Point(10) = {316673.052378,-51498.474955,0.000000,100.000000};
Point(11) = {316731.983113,-51578.589580,0.000000,100.000000};
Point(12) = {316790.012779,-51659.346118,0.000000,100.000000};
Point(13) = {316848.588958,-51739.724012,0.000000,100.000000};
Point(14) = {316909.802527,-51818.013943,0.000000,100.000000};
Point(15) = {316987.531851,-51879.238556,0.000000,100.000000};
Point(16) = {317084.141433,-51891.883270,0.000000,100.000000};
Point(17) = {317169.942578,-51844.006962,0.000000,100.000000};
Point(18) = {317235.134664,-51769.051239,0.000000,100.000000};
Point(19) = {317298.935933,-51692.751871,0.000000,100.000000};
Point(20) = {317362.631654,-51616.367055,0.000000,100.000000};
Point(21) = {317424.310205,-51538.360080,0.000000,100.000000};
Point(22) = {317490.757925,-51464.358006,0.000000,100.000000};
Point(23) = {317557.490306,-51390.609461,0.000000,100.000000};
Point(24) = {317620.446927,-51313.625041,0.000000,100.000000};
Point(25) = {317686.749037,-51239.498060,0.000000,100.000000};
Point(26) = {317753.947671,-51166.170805,0.000000,100.000000};
Point(27) = {317817.341702,-51089.564773,0.000000,100.000000};
Point(28) = {317882.494722,-51014.426737,0.000000,100.000000};
Point(29) = {317948.633085,-50940.144704,0.000000,100.000000};
Point(30) = {318010.697003,-50862.440969,0.000000,100.000000};
Point(31) = {318073.912015,-50785.653789,0.000000,100.000000};
Point(32) = {318134.036722,-50706.435960,0.000000,100.000000};
Point(33) = {318193.769952,-50626.926546,0.000000,100.000000};
Point(34) = {318251.928086,-50546.243623,0.000000,100.000000};
Point(35) = {318305.798546,-50462.665500,0.000000,100.000000};
Point(36) = {318358.749037,-50378.483268,0.000000,100.000000};
Point(37) = {318415.731561,-50297.151756,0.000000,100.000000};
Point(38) = {318480.387501,-50221.577805,0.000000,100.000000};
Point(39) = {318543.209875,-50144.491482,0.000000,100.000000};
Point(40) = {318608.546031,-50069.499936,0.000000,100.000000};
Point(41) = {318670.707015,-49991.870376,0.000000,100.000000};
Point(42) = {318735.445444,-49916.365134,0.000000,100.000000};
Point(43) = {318798.540666,-49839.499857,0.000000,100.000000};
Point(44) = {318863.240522,-49763.999535,0.000000,100.000000};
Point(45) = {318936.624242,-49773.250119,0.000000,100.000000};
Point(46) = {319013.428606,-49836.449139,0.000000,100.000000};
Point(47) = {319085.663491,-49899.897398,0.000000,100.000000};
Point(48) = {319017.671060,-49972.479466,0.000000,100.000000};
Point(49) = {318954.715082,-50049.452098,0.000000,100.000000};
Point(50) = {318890.105287,-50125.065941,0.000000,100.000000};
Point(51) = {318827.054916,-50201.959164,0.000000,100.000000};
Point(52) = {318765.545076,-50280.111152,0.000000,100.000000};
Point(53) = {318704.567345,-50358.677352,0.000000,100.000000};
Point(54) = {318642.486956,-50436.375751,0.000000,100.000000};
Point(55) = {318589.129170,-50520.218816,0.000000,100.000000};
Point(56) = {318551.004342,-50611.818277,0.000000,100.000000};
Point(57) = {318530.364577,-50708.935419,0.000000,100.000000};
Point(58) = {318523.279999,-50808.123649,0.000000,100.000000};
Point(59) = {318519.141424,-50907.500763,0.000000,100.000000};
Point(60) = {318511.227301,-51006.620507,0.000000,100.000000};
Point(61) = {318496.121554,-51104.923209,0.000000,100.000000};
Point(62) = {318480.204650,-51203.094951,0.000000,100.000000};
Point(63) = {318466.823001,-51301.627385,0.000000,100.000000};
Point(64) = {318456.268660,-51400.526600,0.000000,100.000000};
Point(65) = {318444.604822,-51499.278529,0.000000,100.000000};
Point(66) = {318440.710134,-51598.647918,0.000000,100.000000};
Point(67) = {318440.877222,-51698.109473,0.000000,100.000000};
Point(68) = {318425.131260,-51795.946364,0.000000,100.000000};
Point(69) = {318399.771909,-51892.114090,0.000000,100.000000};
Point(70) = {318372.941477,-51987.885418,0.000000,100.000000};
Point(71) = {318342.187251,-52082.452475,0.000000,100.000000};
Point(72) = {318304.918701,-52174.640877,0.000000,100.000000};
Point(73) = {318263.984845,-52265.288028,0.000000,100.000000};
Point(74) = {318222.559656,-52355.687483,0.000000,100.000000};
Point(75) = {318182.234711,-52446.575551,0.000000,100.000000};
Point(76) = {318144.147866,-52538.447350,0.000000,100.000000};
Point(77) = {318110.433312,-52632.017700,0.000000,100.000000};
Point(78) = {318078.614790,-52726.250342,0.000000,100.000000};
Point(79) = {318000.502763,-52729.243026,0.000000,100.000000};
Point(80) = {318003.365795,-52645.590287,0.000000,100.000000};
Point(81) = {318037.173154,-52552.063642,0.000000,100.000000};
Point(82) = {318074.649660,-52459.963605,0.000000,100.000000};
Point(83) = {318117.190399,-52370.064543,0.000000,100.000000};
Point(84) = {318161.268617,-52280.904690,0.000000,100.000000};
Point(85) = {318204.069696,-52191.127691,0.000000,100.000000};
Point(86) = {318243.249351,-52099.720061,0.000000,100.000000};
Point(87) = {318282.361859,-52008.275060,0.000000,100.000000};
Point(88) = {318318.854429,-51915.752516,0.000000,100.000000};
Point(89) = {318348.465532,-51820.915805,0.000000,100.000000};
Point(90) = {318362.358318,-51722.676983,0.000000,100.000000};
Point(91) = {318350.679497,-51624.072548,0.000000,100.000000};
Point(92) = {318311.205516,-51533.636004,0.000000,100.000000};
Point(93) = {318232.007096,-51476.101606,0.000000,100.000000};
Point(94) = {318133.580715,-51467.632517,0.000000,100.000000};
Point(95) = {318037.794744,-51492.847180,0.000000,100.000000};
Point(96) = {317948.722869,-51536.926428,0.000000,100.000000};
Point(97) = {317865.116479,-51590.680042,0.000000,100.000000};
Point(98) = {317786.354264,-51651.338181,0.000000,100.000000};
Point(99) = {317714.127155,-51719.692316,0.000000,100.000000};
Point(100) = {317644.769235,-51790.980359,0.000000,100.000000};
Point(101) = {317576.407589,-51863.227606,0.000000,100.000000};
Point(102) = {317507.887405,-51935.322122,0.000000,100.000000};
Point(103) = {317435.784932,-52003.765219,0.000000,100.000000};
Point(104) = {317348.912227,-52050.976076,0.000000,100.000000};
Point(105) = {317251.469239,-52069.512877,0.000000,100.000000};
Point(106) = {317152.108398,-52072.333621,0.000000,100.000000};
Point(107) = {317053.173482,-52063.128969,0.000000,100.000000};
Point(108) = {316959.823829,-52030.440739,0.000000,100.000000};
Point(109) = {316879.707409,-51971.824423,0.000000,100.000000};
Point(110) = {316810.650963,-51900.522476,0.000000,100.000000};
Point(111) = {316750.977060,-51820.966765,0.000000,100.000000};
Point(112) = {316686.540586,-51745.216651,0.000000,100.000000};
Point(113) = {316620.510307,-51670.834887,0.000000,100.000000};
Point(114) = {316558.351084,-51593.217990,0.000000,100.000000};
Point(115) = {316495.030949,-51516.516929,0.000000,100.000000};
Point(116) = {316433.837019,-51438.122264,0.000000,100.000000};
Point(117) = {316370.707597,-51361.264755,0.000000,100.000000};
Point(118) = {316309.783021,-51282.660125,0.000000,100.000000};
Point(119) = {316247.382623,-51205.210060,0.000000,100.000000};
Point(120) = {316188.842445,-51124.819031,0.000000,100.000000};
Point(121) = {316131.014642,-51043.905389,0.000000,100.000000};
Point(122) = {316070.724694,-50964.801591,0.000000,100.000000};
Point(123) = {316009.007702,-50886.845244,0.000000,100.000000};
Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 5};
Line(5) = {5, 6};
Line(6) = {6, 7};
Line(7) = {7, 8};
Line(8) = {8, 9};
Line(9) = {9, 10};
Line(10) = {10, 11};
Line(11) = {11, 12};
Line(12) = {12, 13};
Line(13) = {13, 14};
Line(14) = {14, 15};
Line(15) = {15, 16};
Line(16) = {16, 17};
Line(17) = {17, 18};
Line(18) = {18, 19};
Line(19) = {19, 20};
Line(20) = {20, 21};
Line(21) = {21, 22};
Line(22) = {22, 23};
Line(23) = {23, 24};
Line(24) = {24, 25};
Line(25) = {25, 26};
Line(26) = {26, 27};
Line(27) = {27, 28};
Line(28) = {28, 29};
Line(29) = {29, 30};
Line(30) = {30, 31};
Line(31) = {31, 32};
Line(32) = {32, 33};
Line(33) = {33, 34};
Line(34) = {34, 35};
Line(35) = {35, 36};
Line(36) = {36, 37};
Line(37) = {37, 38};
Line(38) = {38, 39};
Line(39) = {39, 40};
Line(40) = {40, 41};
Line(41) = {41, 42};
Line(42) = {42, 43};
Line(43) = {43, 44};
Line(44) = {44, 45};
Line(45) = {45, 46};
Line(46) = {46, 47};
Line(47) = {47, 48};
Line(48) = {48, 49};
Line(49) = {49, 50};
Line(50) = {50, 51};
Line(51) = {51, 52};
Line(52) = {52, 53};
Line(53) = {53, 54};
Line(54) = {54, 55};
Line(55) = {55, 56};
Line(56) = {56, 57};
Line(57) = {57, 58};
Line(58) = {58, 59};
Line(59) = {59, 60};
Line(60) = {60, 61};
Line(61) = {61, 62};
Line(62) = {62, 63};
Line(63) = {63, 64};
Line(64) = {64, 65};
Line(65) = {65, 66};
Line(66) = {66, 67};
Line(67) = {67, 68};
Line(68) = {68, 69};
Line(69) = {69, 70};
Line(70) = {70, 71};
Line(71) = {71, 72};
Line(72) = {72, 73};
Line(73) = {73, 74};
Line(74) = {74, 75};
Line(75) = {75, 76};
Line(76) = {76, 77};
Line(77) = {77, 78};
Line(78) = {78, 79};
Line(79) = {79, 80};
Line(80) = {80, 81};
Line(81) = {81, 82};
Line(82) = {82, 83};
Line(83) = {83, 84};
Line(84) = {84, 85};
Line(85) = {85, 86};
Line(86) = {86, 87};
Line(87) = {87, 88};
Line(88) = {88, 89};
Line(89) = {89, 90};
Line(90) = {90, 91};
Line(91) = {91, 92};
Line(92) = {92, 93};
Line(93) = {93, 94};
Line(94) = {94, 95};
Line(95) = {95, 96};
Line(96) = {96, 97};
Line(97) = {97, 98};
Line(98) = {98, 99};
Line(99) = {99, 100};
Line(100) = {100, 101};
Line(101) = {101, 102};
Line(102) = {102, 103};
Line(103) = {103, 104};
Line(104) = {104, 105};
Line(105) = {105, 106};
Line(106) = {106, 107};
Line(107) = {107, 108};
Line(108) = {108, 109};
Line(109) = {109, 110};
Line(110) = {110, 111};
Line(111) = {111, 112};
Line(112) = {112, 113};
Line(113) = {113, 114};
Line(114) = {114, 115};
Line(115) = {115, 116};
Line(116) = {116, 117};
Line(117) = {117, 118};
Line(118) = {118, 119};
Line(119) = {119, 120};
Line(120) = {120, 121};
Line(121) = {121, 122};
Line(122) = {122, 123};
Line(123) = {123, 1};
Line Loop(1) = {1:123};
Point(124) = {317414.577311,-51906.840017,0.000000,100.000000};
Point(125) = {317495.777759,-51862.221335,0.000000,100.000000};
Point(126) = {317557.284095,-51786.133661,0.000000,100.000000};
Point(127) = {317618.004967,-51709.476259,0.000000,100.000000};
Point(128) = {317688.453813,-51641.668768,0.000000,100.000000};
Point(129) = {317762.030189,-51577.143125,0.000000,100.000000};
Point(130) = {317838.386260,-51516.026949,0.000000,100.000000};
Point(131) = {317923.764148,-51468.463694,0.000000,100.000000};
Point(132) = {318015.935764,-51435.999938,0.000000,100.000000};
Point(133) = {318112.251593,-51419.635113,0.000000,100.000000};
Point(134) = {318210.018130,-51420.881335,0.000000,100.000000};
Point(135) = {318307.402836,-51427.580614,0.000000,100.000000};
Point(136) = {318381.348845,-51369.368487,0.000000,100.000000};
Point(137) = {318410.397240,-51276.553013,0.000000,100.000000};
Point(138) = {318421.437129,-51179.327272,0.000000,100.000000};
Point(139) = {318434.467150,-51082.348898,0.000000,100.000000};
Point(140) = {318445.329888,-50985.098497,0.000000,100.000000};
Point(141) = {318455.175825,-50887.731516,0.000000,100.000000};
Point(142) = {318459.051172,-50790.526064,0.000000,100.000000};
Point(143) = {318390.152380,-50734.727694,0.000000,100.000000};
Point(144) = {318315.414706,-50792.871272,0.000000,100.000000};
Point(145) = {318254.297901,-50869.299127,0.000000,100.000000};
Point(146) = {318192.460271,-50945.147392,0.000000,100.000000};
Point(147) = {318133.102146,-51022.940166,0.000000,100.000000};
Point(148) = {318070.122293,-51097.841640,0.000000,100.000000};
Point(149) = {318008.428825,-51173.793778,0.000000,100.000000};
Point(150) = {317944.534578,-51247.884758,0.000000,100.000000};
Point(151) = {317876.774916,-51318.494670,0.000000,100.000000};
Point(152) = {317807.691041,-51387.812713,0.000000,100.000000};
Point(153) = {317738.203869,-51456.727911,0.000000,100.000000};
Point(154) = {317668.279971,-51525.200052,0.000000,100.000000};
Point(155) = {317597.820266,-51593.120255,0.000000,100.000000};
Point(156) = {317527.145382,-51660.815894,0.000000,100.000000};
Point(157) = {317457.573886,-51729.634993,0.000000,100.000000};
Point(158) = {317395.914500,-51805.449453,0.000000,100.000000};
Point(159) = {317392.211569,-51894.425165,0.000000,100.000000};
Line(124) = {124, 125};
Line(125) = {125, 126};
Line(126) = {126, 127};
Line(127) = {127, 128};
Line(128) = {128, 129};
Line(129) = {129, 130};
Line(130) = {130, 131};
Line(131) = {131, 132};
Line(132) = {132, 133};
Line(133) = {133, 134};
Line(134) = {134, 135};
Line(135) = {135, 136};
Line(136) = {136, 137};
Line(137) = {137, 138};
Line(138) = {138, 139};
Line(139) = {139, 140};
Line(140) = {140, 141};
Line(141) = {141, 142};
Line(142) = {142, 143};
Line(143) = {143, 144};
Line(144) = {144, 145};
Line(145) = {145, 146};
Line(146) = {146, 147};
Line(147) = {147, 148};
Line(148) = {148, 149};
Line(149) = {149, 150};
Line(150) = {150, 151};
Line(151) = {151, 152};
Line(152) = {152, 153};
Line(153) = {153, 154};
Line(154) = {154, 155};
Line(155) = {155, 156};
Line(156) = {156, 157};
Line(157) = {157, 158};
Line(158) = {158, 159};
Line(159) = {159, 124};
Line Loop(2) = {124:159};
Plane Surface(1) = {1:2};
Mesh.Remove4Triangles = 1;
-------------- next part --------------
A non-text attachment was scrubbed...
Name: kubu-bifuraction-scale-00001-maxres-00100-minres-00100-ratio-00001.msh
Type: model/mesh
Size: 22596 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20160928/07def849/attachment-0001.msh>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: kubu-bifurcation-trisected-element.png
Type: image/png
Size: 65877 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20160928/07def849/attachment-0001.png>