From 9f6c64a4f9eb45a702e2c0321972704df225f6d8 Mon Sep 17 00:00:00 2001
From: "Hochsteger, Matthias" <mhochsteger@cerbsim.com>
Date: Fri, 27 Sep 2024 01:00:50 +0200
Subject: [PATCH] Fix splitimprove edge

---
 libsrc/meshing/improve3.cpp |  12 +-
 libsrc/meshing/improve3.hpp |  11 --
 tests/pytest/results.json   | 356 ++++++++++++++++++------------------
 3 files changed, 182 insertions(+), 197 deletions(-)

diff --git a/libsrc/meshing/improve3.cpp b/libsrc/meshing/improve3.cpp
index 620a650b..ac6fbb4d 100644
--- a/libsrc/meshing/improve3.cpp
+++ b/libsrc/meshing/improve3.cpp
@@ -526,8 +526,6 @@ double MeshOptimize3d :: SplitImproveEdge (Table<ElementIndex,PointIndex> & elem
 
   double bad2 = pf.Func (px);
 
-  mesh[ptmp] = Point<3>(pnew);
-
   for (int k = 0; k < hasbothpoints.Size(); k++)
     {
       Element & oldel = mesh[hasbothpoints[k]];
@@ -543,11 +541,12 @@ double MeshOptimize3d :: SplitImproveEdge (Table<ElementIndex,PointIndex> & elem
           if (newel2[l] == pi1) newel2[l] = ptmp;
         }
 
-      if (!mesh.LegalTet (oldel)) bad1 += GetLegalPenalty();
-      if (!mesh.LegalTet (newel1)) bad2 += GetLegalPenalty();
-      if (!mesh.LegalTet (newel2)) bad2 += GetLegalPenalty();
+      if (!mesh.LegalTet (oldel)) return 0.0;
+      if (!mesh.LegalTet (newel1)) return 0.0;
+      if (!mesh.LegalTet (newel2)) return 0.0;
     }
 
+  if(bad2 >= 1e24) return 0.0;
   d_badness = bad2-bad1;
   if(check_only)
       return d_badness;
@@ -573,9 +572,6 @@ double MeshOptimize3d :: SplitImproveEdge (Table<ElementIndex,PointIndex> & elem
               if (newel2[l] == pi1) newel2[l] = pinew;
             }
 
-          if( newel1.Volume(mesh.Points()) < 0.0 || newel2.Volume(mesh.Points()) < 0.0)
-            return 0.0;
-
           oldel.Touch();
           oldel.Delete();
 
diff --git a/libsrc/meshing/improve3.hpp b/libsrc/meshing/improve3.hpp
index e35ac924..ceca4c26 100644
--- a/libsrc/meshing/improve3.hpp
+++ b/libsrc/meshing/improve3.hpp
@@ -4,10 +4,6 @@
 namespace netgen
 {
 
-extern double CalcTotalBad (const Mesh::T_POINTS & points, 
-			    const Array<Element, ElementIndex> & elements,
-			    const MeshingParameters & mp);
-
 
 ///
 class MeshOptimize3d
@@ -64,13 +60,6 @@ public:
   }
 
 
-  double CalcTotalBad (const Mesh::T_POINTS & points, 
-		       const Array<Element, ElementIndex> & elements)
-  {
-    return netgen::CalcTotalBad (points, elements, mp);
-  }
-
-
   double GetLegalPenalty()
   {
     return goal == OPT_LEGAL ? 1e15 : 1e6;
diff --git a/tests/pytest/results.json b/tests/pytest/results.json
index c455838d..a466bb45 100644
--- a/tests/pytest/results.json
+++ b/tests/pytest/results.json
@@ -124,8 +124,8 @@
         },
         {
             "angles_tet": [
-                14.935,
-                158.04
+                13.199,
+                159.74
             ],
             "angles_trig": [
                 19.228,
@@ -134,8 +134,8 @@
             "ne1d": 136,
             "ne2d": 204,
             "ne3d": 298,
-            "quality_histogram": "[0, 0, 0, 2, 4, 7, 10, 10, 9, 12, 13, 23, 22, 27, 43, 47, 35, 25, 8, 1]",
-            "total_badness": 498.40723523
+            "quality_histogram": "[0, 0, 0, 2, 3, 4, 9, 12, 12, 19, 23, 12, 26, 30, 36, 38, 35, 26, 11, 0]",
+            "total_badness": 500.45563074
         },
         {
             "angles_tet": [
@@ -293,33 +293,33 @@
         },
         {
             "angles_tet": [
-                12.894,
-                161.49
+                3.7664,
+                169.51
             ],
             "angles_trig": [
-                14.156,
-                150.69
+                12.979,
+                139.05
             ],
             "ne1d": 32,
             "ne2d": 208,
-            "ne3d": 477,
-            "quality_histogram": "[0, 0, 0, 1, 4, 9, 20, 22, 19, 29, 46, 33, 45, 51, 50, 47, 41, 32, 20, 8]",
-            "total_badness": 817.64355327
+            "ne3d": 268,
+            "quality_histogram": "[0, 1, 3, 3, 6, 18, 39, 44, 21, 10, 14, 10, 8, 24, 12, 20, 14, 12, 9, 0]",
+            "total_badness": 632.29269253
         },
         {
             "angles_tet": [
-                8.5678,
-                164.24
+                11.644,
+                159.22
             ],
             "angles_trig": [
-                12.521,
-                143.3
+                15.078,
+                146.22
             ],
             "ne1d": 48,
             "ne2d": 420,
-            "ne3d": 626,
-            "quality_histogram": "[0, 0, 0, 10, 8, 19, 21, 14, 37, 71, 91, 75, 76, 62, 48, 37, 18, 18, 18, 3]",
-            "total_badness": 1193.0076867
+            "ne3d": 524,
+            "quality_histogram": "[0, 0, 0, 0, 5, 17, 10, 10, 36, 63, 83, 67, 66, 54, 50, 29, 12, 14, 8, 0]",
+            "total_badness": 962.09837603
         },
         {
             "angles_tet": [
@@ -471,9 +471,9 @@
             ],
             "ne1d": 262,
             "ne2d": 652,
-            "ne3d": 1850,
-            "quality_histogram": "[0, 2, 4, 19, 46, 79, 87, 113, 85, 41, 56, 75, 136, 174, 213, 238, 192, 160, 101, 29]",
-            "total_badness": 3376.5621892
+            "ne3d": 1831,
+            "quality_histogram": "[0, 2, 4, 16, 42, 82, 88, 111, 82, 36, 58, 73, 138, 172, 211, 238, 191, 161, 97, 29]",
+            "total_badness": 3327.3635018
         },
         {
             "angles_tet": [
@@ -516,9 +516,9 @@
             ],
             "ne1d": 262,
             "ne2d": 652,
-            "ne3d": 1653,
-            "quality_histogram": "[0, 0, 2, 15, 35, 67, 71, 93, 82, 34, 41, 56, 104, 138, 185, 235, 222, 147, 101, 25]",
-            "total_badness": 2904.5831463
+            "ne3d": 1641,
+            "quality_histogram": "[0, 0, 3, 13, 35, 76, 82, 101, 77, 29, 36, 54, 101, 133, 179, 232, 215, 148, 102, 25]",
+            "total_badness": 2926.5758078
         },
         {
             "angles_tet": [
@@ -753,18 +753,18 @@
         },
         {
             "angles_tet": [
-                10.212,
-                162.36
+                2.8131,
+                175.86
             ],
             "angles_trig": [
-                9.0545,
-                133.12
+                10.214,
+                139.62
             ],
             "ne1d": 44,
             "ne2d": 142,
-            "ne3d": 358,
-            "quality_histogram": "[0, 0, 1, 3, 8, 9, 27, 25, 43, 45, 41, 38, 32, 29, 23, 19, 7, 5, 3, 0]",
-            "total_badness": 752.80521692
+            "ne3d": 344,
+            "quality_histogram": "[1, 0, 3, 6, 6, 20, 30, 33, 37, 44, 40, 32, 20, 22, 23, 12, 8, 3, 4, 0]",
+            "total_badness": 801.60189009
         },
         {
             "angles_tet": [
@@ -854,9 +854,9 @@
             ],
             "ne1d": 48,
             "ne2d": 100,
-            "ne3d": 102,
-            "quality_histogram": "[0, 0, 0, 2, 1, 4, 12, 15, 11, 11, 10, 5, 4, 3, 6, 10, 8, 0, 0, 0]",
-            "total_badness": 224.10472313
+            "ne3d": 94,
+            "quality_histogram": "[0, 0, 0, 2, 1, 4, 12, 13, 8, 8, 6, 6, 4, 4, 6, 10, 8, 2, 0, 0]",
+            "total_badness": 203.19030185
         },
         {
             "angles_tet": [
@@ -922,18 +922,18 @@
         },
         {
             "angles_tet": [
-                4.1029,
-                170.7
+                5.2584,
+                171.14
             ],
             "angles_trig": [
-                6.1932,
-                162.02
+                13.512,
+                152.22
             ],
             "ne1d": 0,
             "ne2d": 156,
-            "ne3d": 617,
-            "quality_histogram": "[0, 6, 38, 56, 76, 81, 63, 54, 37, 40, 32, 31, 21, 29, 21, 18, 10, 1, 2, 1]",
-            "total_badness": 2089.5884695
+            "ne3d": 264,
+            "quality_histogram": "[0, 4, 11, 25, 53, 58, 33, 23, 23, 12, 7, 5, 4, 4, 1, 0, 1, 0, 0, 0]",
+            "total_badness": 995.21455744
         },
         {
             "angles_tet": [
@@ -1014,7 +1014,7 @@
         },
         {
             "angles_tet": [
-                18.758,
+                18.757,
                 150.89
             ],
             "angles_trig": [
@@ -1023,9 +1023,9 @@
             ],
             "ne1d": 86,
             "ne2d": 336,
-            "ne3d": 459,
-            "quality_histogram": "[0, 0, 0, 0, 0, 0, 2, 2, 10, 21, 41, 41, 53, 60, 65, 58, 48, 41, 11, 6]",
-            "total_badness": 690.5842012
+            "ne3d": 453,
+            "quality_histogram": "[0, 0, 0, 0, 0, 0, 2, 2, 9, 17, 38, 40, 52, 65, 61, 58, 49, 42, 10, 8]",
+            "total_badness": 676.43193266
         },
         {
             "angles_tet": [
@@ -1038,9 +1038,9 @@
             ],
             "ne1d": 130,
             "ne2d": 794,
-            "ne3d": 1383,
-            "quality_histogram": "[0, 0, 0, 0, 1, 1, 12, 28, 56, 58, 69, 92, 155, 145, 190, 170, 169, 125, 94, 18]",
-            "total_badness": 2076.8283951
+            "ne3d": 1385,
+            "quality_histogram": "[0, 0, 0, 0, 1, 2, 12, 28, 58, 58, 77, 88, 158, 144, 191, 155, 170, 136, 90, 17]",
+            "total_badness": 2086.8296552
         },
         {
             "angles_tet": [
@@ -1252,9 +1252,9 @@
             ],
             "ne1d": 276,
             "ne2d": 544,
-            "ne3d": 603,
-            "quality_histogram": "[0, 0, 0, 0, 0, 0, 1, 8, 26, 43, 52, 60, 73, 62, 74, 62, 77, 43, 15, 7]",
-            "total_badness": 934.85096826
+            "ne3d": 598,
+            "quality_histogram": "[0, 0, 0, 0, 0, 0, 1, 8, 26, 43, 52, 58, 68, 66, 78, 63, 75, 42, 12, 6]",
+            "total_badness": 928.42108853
         }
     ],
     "fichera.geo": [
@@ -1367,8 +1367,8 @@
         },
         {
             "angles_tet": [
-                4.8415,
-                163.34
+                4.0123,
+                173.21
             ],
             "angles_trig": [
                 8.9881,
@@ -1376,13 +1376,13 @@
             ],
             "ne1d": 298,
             "ne2d": 502,
-            "ne3d": 598,
-            "quality_histogram": "[0, 0, 1, 4, 10, 14, 29, 39, 46, 40, 56, 63, 61, 50, 43, 60, 37, 26, 14, 5]",
-            "total_badness": 1132.0903047
+            "ne3d": 600,
+            "quality_histogram": "[0, 1, 1, 4, 10, 14, 27, 40, 49, 43, 56, 62, 63, 49, 43, 59, 36, 24, 14, 5]",
+            "total_badness": 1150.0886518
         },
         {
             "angles_tet": [
-                12.929,
+                11.259,
                 152.0
             ],
             "angles_trig": [
@@ -1391,9 +1391,9 @@
             ],
             "ne1d": 370,
             "ne2d": 758,
-            "ne3d": 980,
-            "quality_histogram": "[0, 0, 0, 0, 0, 8, 8, 23, 33, 29, 49, 79, 112, 130, 144, 142, 95, 67, 49, 12]",
-            "total_badness": 1495.506642
+            "ne3d": 970,
+            "quality_histogram": "[0, 0, 0, 0, 1, 10, 9, 24, 34, 28, 49, 77, 110, 126, 142, 142, 93, 64, 49, 12]",
+            "total_badness": 1491.287287
         },
         {
             "angles_tet": [
@@ -1427,7 +1427,7 @@
         },
         {
             "angles_tet": [
-                18.283,
+                18.259,
                 142.44
             ],
             "angles_trig": [
@@ -1436,9 +1436,9 @@
             ],
             "ne1d": 1862,
             "ne2d": 18540,
-            "ne3d": 97376,
-            "quality_histogram": "[0, 0, 0, 0, 0, 0, 1, 1, 18, 74, 251, 659, 1770, 4375, 8501, 14080, 20129, 22819, 18498, 6200]",
-            "total_badness": 117815.89866
+            "ne3d": 97369,
+            "quality_histogram": "[0, 0, 0, 0, 0, 0, 1, 1, 18, 73, 250, 659, 1775, 4369, 8501, 14086, 20127, 22811, 18498, 6200]",
+            "total_badness": 117806.4728
         }
     ],
     "lense.in2d": [
@@ -1637,9 +1637,9 @@
             ],
             "ne1d": 5886,
             "ne2d": 45882,
-            "ne3d": 154702,
-            "quality_histogram": "[0, 0, 0, 0, 0, 1, 8, 66, 222, 737, 2094, 5786, 10315, 16832, 23407, 25576, 26164, 23123, 16303, 4068]",
-            "total_badness": 202905.19298
+            "ne3d": 154731,
+            "quality_histogram": "[0, 0, 0, 0, 0, 1, 8, 66, 223, 738, 2099, 5794, 10308, 16876, 23404, 25577, 26145, 23113, 16306, 4073]",
+            "total_badness": 202956.41935
         },
         {
             "angles_tet": [
@@ -1652,9 +1652,9 @@
             ],
             "ne1d": 2746,
             "ne2d": 10428,
-            "ne3d": 22177,
-            "quality_histogram": "[0, 0, 0, 1, 2, 10, 37, 141, 314, 667, 1312, 2125, 2625, 3088, 2860, 2703, 2437, 2076, 1429, 350]",
-            "total_badness": 32293.784552
+            "ne3d": 22190,
+            "quality_histogram": "[0, 0, 0, 1, 2, 10, 37, 144, 313, 674, 1311, 2118, 2650, 3080, 2854, 2712, 2432, 2069, 1426, 357]",
+            "total_badness": 32325.6857
         },
         {
             "angles_tet": [
@@ -1667,9 +1667,9 @@
             ],
             "ne1d": 4106,
             "ne2d": 23238,
-            "ne3d": 56933,
-            "quality_histogram": "[0, 0, 0, 0, 31, 64, 210, 334, 666, 1194, 2196, 3631, 5728, 7499, 8338, 8581, 7869, 6084, 3580, 928]",
-            "total_badness": 80610.044111
+            "ne3d": 56938,
+            "quality_histogram": "[0, 0, 0, 0, 31, 62, 209, 334, 665, 1195, 2198, 3630, 5732, 7498, 8341, 8581, 7870, 6084, 3580, 928]",
+            "total_badness": 80610.741802
         }
     ],
     "manyholes2.geo": [
@@ -1684,9 +1684,9 @@
             ],
             "ne1d": 10202,
             "ne2d": 41054,
-            "ne3d": 97037,
-            "quality_histogram": "[0, 0, 0, 0, 3, 14, 76, 230, 646, 1720, 3702, 6469, 9523, 12171, 13058, 13512, 13508, 12168, 8074, 2163]",
-            "total_badness": 134237.53167
+            "ne3d": 97025,
+            "quality_histogram": "[0, 0, 0, 0, 3, 15, 78, 230, 659, 1718, 3706, 6461, 9551, 12155, 13121, 13548, 13473, 12121, 8033, 2153]",
+            "total_badness": 134274.87401
         }
     ],
     "matrix.geo": [
@@ -1707,8 +1707,8 @@
         },
         {
             "angles_tet": [
-                6.4945,
-                166.83
+                5.0793,
+                171.79
             ],
             "angles_trig": [
                 8.2262,
@@ -1716,14 +1716,14 @@
             ],
             "ne1d": 106,
             "ne2d": 314,
-            "ne3d": 876,
-            "quality_histogram": "[0, 0, 5, 33, 47, 50, 72, 68, 97, 92, 81, 76, 63, 51, 44, 32, 29, 28, 7, 1]",
-            "total_badness": 2078.9282016
+            "ne3d": 843,
+            "quality_histogram": "[0, 1, 7, 36, 52, 60, 74, 75, 96, 90, 86, 69, 55, 39, 34, 27, 22, 18, 2, 0]",
+            "total_badness": 2121.4782246
         },
         {
             "angles_tet": [
-                6.3225,
-                170.81
+                3.2223,
+                172.77
             ],
             "angles_trig": [
                 10.133,
@@ -1731,9 +1731,9 @@
             ],
             "ne1d": 132,
             "ne2d": 588,
-            "ne3d": 1733,
-            "quality_histogram": "[0, 0, 2, 14, 30, 79, 151, 126, 80, 114, 115, 148, 155, 191, 164, 139, 98, 75, 43, 9]",
-            "total_badness": 3419.5579224
+            "ne3d": 1721,
+            "quality_histogram": "[0, 1, 5, 17, 33, 72, 137, 128, 86, 121, 130, 155, 144, 192, 159, 125, 92, 74, 41, 9]",
+            "total_badness": 3443.2721678
         },
         {
             "angles_tet": [
@@ -1968,8 +1968,8 @@
         },
         {
             "angles_tet": [
-                12.454,
-                162.07
+                8.8923,
+                165.12
             ],
             "angles_trig": [
                 10.412,
@@ -1977,24 +1977,24 @@
             ],
             "ne1d": 160,
             "ne2d": 234,
-            "ne3d": 390,
-            "quality_histogram": "[0, 0, 0, 5, 7, 13, 24, 18, 26, 26, 32, 36, 45, 41, 22, 33, 30, 21, 10, 1]",
-            "total_badness": 745.1543892
+            "ne3d": 372,
+            "quality_histogram": "[0, 0, 0, 7, 6, 15, 27, 17, 25, 30, 34, 36, 43, 38, 18, 24, 27, 17, 7, 1]",
+            "total_badness": 738.96263535
         },
         {
             "angles_tet": [
-                13.042,
-                161.3
+                10.154,
+                164.25
             ],
             "angles_trig": [
-                13.898,
-                141.37
+                15.73,
+                145.3
             ],
             "ne1d": 232,
             "ne2d": 494,
-            "ne3d": 1188,
-            "quality_histogram": "[0, 0, 0, 0, 8, 21, 41, 30, 61, 76, 112, 110, 137, 141, 109, 124, 85, 89, 38, 6]",
-            "total_badness": 2016.0253302
+            "ne3d": 1001,
+            "quality_histogram": "[0, 0, 0, 2, 16, 39, 50, 69, 58, 62, 85, 92, 106, 102, 99, 91, 59, 49, 18, 4]",
+            "total_badness": 1873.6285579
         },
         {
             "angles_tet": [
@@ -2046,22 +2046,22 @@
         {
             "angles_tet": [
                 1.2182,
-                170.22
+                165.85
             ],
             "angles_trig": [
                 1.8555,
-                148.54
+                160.04
             ],
             "ne1d": 892,
             "ne2d": 2220,
-            "ne3d": 6353,
-            "quality_histogram": "[4, 9, 34, 34, 43, 54, 51, 62, 96, 128, 246, 338, 551, 714, 899, 958, 936, 708, 385, 103]",
-            "total_badness": 9830.4969118
+            "ne3d": 6361,
+            "quality_histogram": "[4, 14, 29, 34, 42, 56, 52, 60, 94, 131, 241, 341, 537, 709, 899, 953, 944, 728, 387, 106]",
+            "total_badness": 9865.9905616
         },
         {
             "angles_tet": [
-                1.3667,
-                170.92
+                0.11753,
+                179.78
             ],
             "angles_trig": [
                 2.5035,
@@ -2069,14 +2069,14 @@
             ],
             "ne1d": 572,
             "ne2d": 742,
-            "ne3d": 959,
-            "quality_histogram": "[2, 21, 50, 69, 85, 99, 78, 82, 76, 76, 57, 54, 65, 52, 36, 25, 15, 12, 4, 1]",
-            "total_badness": 3180.1990985
+            "ne3d": 931,
+            "quality_histogram": "[6, 25, 55, 67, 80, 92, 69, 76, 70, 67, 59, 55, 63, 53, 38, 24, 16, 11, 5, 0]",
+            "total_badness": 3693.358071
         },
         {
             "angles_tet": [
                 1.1094,
-                172.05
+                172.08
             ],
             "angles_trig": [
                 3.4703,
@@ -2084,9 +2084,9 @@
             ],
             "ne1d": 724,
             "ne2d": 1340,
-            "ne3d": 2199,
-            "quality_histogram": "[3, 12, 31, 52, 50, 61, 53, 72, 84, 126, 178, 184, 234, 254, 262, 206, 180, 102, 48, 7]",
-            "total_badness": 4492.7256936
+            "ne3d": 2183,
+            "quality_histogram": "[3, 13, 33, 52, 47, 59, 55, 66, 82, 121, 178, 184, 232, 255, 254, 206, 182, 107, 48, 6]",
+            "total_badness": 4468.3231535
         },
         {
             "angles_tet": [
@@ -2101,7 +2101,7 @@
             "ne2d": 2330,
             "ne3d": 6152,
             "quality_histogram": "[3, 9, 26, 42, 39, 53, 62, 67, 77, 120, 158, 230, 393, 629, 835, 972, 1018, 839, 453, 127]",
-            "total_badness": 9316.7375577
+            "total_badness": 9316.7379917
         },
         {
             "angles_tet": [
@@ -2114,9 +2114,9 @@
             ],
             "ne1d": 1554,
             "ne2d": 5646,
-            "ne3d": 23496,
-            "quality_histogram": "[2, 5, 10, 10, 20, 55, 63, 51, 78, 136, 198, 405, 840, 1576, 2705, 3890, 4717, 4583, 3248, 904]",
-            "total_badness": 30261.471634
+            "ne3d": 23451,
+            "quality_histogram": "[2, 5, 10, 10, 20, 55, 63, 51, 78, 135, 197, 409, 833, 1584, 2686, 3853, 4682, 4599, 3278, 901]",
+            "total_badness": 30195.253166
         },
         {
             "angles_tet": [
@@ -2129,9 +2129,9 @@
             ],
             "ne1d": 2992,
             "ne2d": 22730,
-            "ne3d": 206422,
-            "quality_histogram": "[4, 6, 13, 10, 14, 19, 26, 63, 87, 165, 481, 1163, 3333, 7807, 16722, 29390, 42845, 50235, 40650, 13389]",
-            "total_badness": 249062.5796
+            "ne3d": 206391,
+            "quality_histogram": "[4, 6, 13, 10, 14, 19, 26, 63, 87, 166, 484, 1162, 3348, 7835, 16694, 29391, 42810, 50220, 40636, 13403]",
+            "total_badness": 249038.12593
         }
     ],
     "revolution.geo": [
@@ -2167,18 +2167,18 @@
         },
         {
             "angles_tet": [
-                17.869,
-                147.84
+                17.637,
+                150.4
             ],
             "angles_trig": [
-                17.577,
-                134.89
+                17.573,
+                134.96
             ],
             "ne1d": 240,
             "ne2d": 1584,
-            "ne3d": 3411,
-            "quality_histogram": "[0, 0, 0, 0, 0, 0, 2, 11, 53, 110, 240, 368, 463, 499, 461, 406, 324, 266, 163, 45]",
-            "total_badness": 5043.717387
+            "ne3d": 3395,
+            "quality_histogram": "[0, 0, 0, 0, 0, 0, 2, 12, 54, 120, 224, 382, 446, 488, 441, 388, 346, 275, 178, 39]",
+            "total_badness": 5015.8391018
         },
         {
             "angles_tet": [
@@ -2321,18 +2321,18 @@
     "shaft.geo": [
         {
             "angles_tet": [
-                5.9923,
-                168.9
+                9.5959,
+                166.59
             ],
             "angles_trig": [
-                10.45,
-                155.97
+                13.405,
+                152.57
             ],
             "ne1d": 708,
             "ne2d": 1656,
-            "ne3d": 2645,
-            "quality_histogram": "[0, 1, 14, 31, 44, 52, 49, 56, 80, 150, 244, 358, 301, 220, 224, 284, 254, 169, 90, 24]",
-            "total_badness": 4677.5910875
+            "ne3d": 2471,
+            "quality_histogram": "[0, 0, 0, 7, 32, 18, 13, 36, 58, 141, 235, 358, 321, 223, 218, 274, 249, 175, 90, 23]",
+            "total_badness": 4010.0976507
         },
         {
             "angles_tet": [
@@ -2341,28 +2341,28 @@
             ],
             "angles_trig": [
                 10.094,
-                120.04
+                123.18
             ],
             "ne1d": 410,
             "ne2d": 542,
-            "ne3d": 717,
-            "quality_histogram": "[0, 0, 0, 2, 1, 4, 3, 9, 19, 28, 37, 57, 65, 82, 74, 80, 80, 88, 62, 26]",
-            "total_badness": 1058.1207822
+            "ne3d": 615,
+            "quality_histogram": "[0, 0, 0, 3, 1, 4, 8, 16, 37, 54, 65, 67, 67, 71, 59, 58, 50, 43, 9, 3]",
+            "total_badness": 1035.1361164
         },
         {
             "angles_tet": [
-                13.382,
-                157.07
+                1.088,
+                177.95
             ],
             "angles_trig": [
-                14.807,
-                149.61
+                9.6019,
+                144.68
             ],
             "ne1d": 510,
             "ne2d": 912,
-            "ne3d": 1820,
-            "quality_histogram": "[0, 0, 0, 2, 7, 27, 46, 83, 98, 117, 134, 154, 181, 195, 223, 189, 172, 107, 66, 19]",
-            "total_badness": 3049.4355671
+            "ne3d": 1258,
+            "quality_histogram": "[3, 0, 0, 0, 16, 51, 82, 87, 46, 39, 68, 93, 108, 122, 155, 141, 136, 64, 36, 11]",
+            "total_badness": 2595.8155702
         },
         {
             "angles_tet": [
@@ -2370,14 +2370,14 @@
                 166.59
             ],
             "angles_trig": [
-                13.661,
-                152.23
+                13.405,
+                152.57
             ],
             "ne1d": 708,
             "ne2d": 1656,
-            "ne3d": 2553,
-            "quality_histogram": "[0, 0, 0, 1, 8, 3, 11, 32, 83, 164, 227, 330, 314, 265, 249, 303, 271, 182, 88, 22]",
-            "total_badness": 4003.23251
+            "ne3d": 2391,
+            "quality_histogram": "[0, 0, 0, 5, 31, 10, 9, 21, 66, 148, 209, 332, 268, 233, 234, 277, 242, 190, 90, 26]",
+            "total_badness": 3820.2725527
         },
         {
             "angles_tet": [
@@ -2535,8 +2535,8 @@
         },
         {
             "angles_tet": [
-                6.4483,
-                169.86
+                5.8125,
+                172.01
             ],
             "angles_trig": [
                 9.6618,
@@ -2544,9 +2544,9 @@
             ],
             "ne1d": 30,
             "ne2d": 100,
-            "ne3d": 242,
-            "quality_histogram": "[0, 0, 6, 15, 20, 43, 45, 12, 21, 10, 21, 8, 12, 4, 4, 6, 4, 7, 3, 1]",
-            "total_badness": 732.39744461
+            "ne3d": 236,
+            "quality_histogram": "[0, 1, 6, 15, 19, 40, 50, 12, 19, 9, 16, 8, 12, 4, 4, 6, 4, 7, 3, 1]",
+            "total_badness": 725.69720862
         },
         {
             "angles_tet": [
@@ -2888,18 +2888,18 @@
         },
         {
             "angles_tet": [
-                2.7899,
-                172.79
+                2.2624,
+                174.01
             ],
             "angles_trig": [
-                5.1656,
-                163.81
+                9.0702,
+                156.78
             ],
             "ne1d": 0,
             "ne2d": 648,
-            "ne3d": 3545,
-            "quality_histogram": "[7, 92, 302, 447, 533, 448, 365, 262, 259, 198, 148, 117, 83, 62, 67, 57, 45, 23, 21, 9]",
-            "total_badness": 14086.983422
+            "ne3d": 1316,
+            "quality_histogram": "[4, 69, 112, 154, 262, 273, 143, 70, 53, 27, 33, 27, 23, 20, 14, 11, 7, 6, 3, 5]",
+            "total_badness": 5967.2368224
         },
         {
             "angles_tet": [
@@ -2981,7 +2981,7 @@
         {
             "angles_tet": [
                 2.8238,
-                174.1
+                175.21
             ],
             "angles_trig": [
                 7.7605,
@@ -2989,9 +2989,9 @@
             ],
             "ne1d": 390,
             "ne2d": 516,
-            "ne3d": 1330,
-            "quality_histogram": "[0, 1, 5, 14, 16, 42, 76, 118, 126, 151, 160, 134, 120, 106, 86, 80, 53, 28, 11, 3]",
-            "total_badness": 2746.9153651
+            "ne3d": 1332,
+            "quality_histogram": "[0, 2, 5, 14, 16, 42, 76, 118, 126, 153, 159, 133, 121, 106, 87, 80, 53, 27, 11, 3]",
+            "total_badness": 2769.107096
         },
         {
             "angles_tet": [
@@ -3257,32 +3257,32 @@
         {
             "angles_tet": [
                 18.044,
-                148.8
+                142.71
             ],
             "angles_trig": [
-                21.63,
-                119.03
+                25.435,
+                118.39
             ],
             "ne1d": 68,
             "ne2d": 98,
-            "ne3d": 120,
-            "quality_histogram": "[0, 0, 0, 0, 0, 1, 1, 4, 8, 14, 6, 12, 4, 13, 10, 11, 15, 17, 3, 1]",
-            "total_badness": 193.19570482
+            "ne3d": 106,
+            "quality_histogram": "[0, 0, 0, 0, 0, 1, 0, 4, 6, 5, 5, 13, 5, 15, 7, 9, 15, 18, 3, 0]",
+            "total_badness": 164.37650945
         },
         {
             "angles_tet": [
-                12.305,
-                161.08
+                11.996,
+                161.96
             ],
             "angles_trig": [
-                15.704,
-                148.08
+                16.697,
+                138.95
             ],
             "ne1d": 102,
             "ne2d": 234,
-            "ne3d": 445,
-            "quality_histogram": "[0, 0, 0, 0, 5, 11, 17, 23, 47, 28, 34, 35, 45, 31, 29, 41, 44, 30, 19, 6]",
-            "total_badness": 795.41761043
+            "ne3d": 317,
+            "quality_histogram": "[0, 0, 0, 6, 7, 26, 26, 33, 44, 25, 11, 12, 6, 19, 16, 18, 42, 16, 8, 2]",
+            "total_badness": 688.54847792
         },
         {
             "angles_tet": [