[ci] add test cgals2m

Author: fangq

README.md

@@ -199,7 +199,7 @@ tracked in https://github.com/NeuroJSON/pyiso2mesh/issues/1
 | > Iso2mesh main function backend| |         | ✅ `readasc.m` | ⭕️ tested |
 | ✅ `binsurface.m` | ✅ tested |             | ⭕️ `readinr.m` | ⭕️ tested |
 | ✅ `cgalv2m.m` | ✅ tested |                | ✅ `readmedit.m` | ⭕️ tested |
-| ✅ `cgals2m.m` | ⭕️ tested |                | ✅ `readoff.m` | ✅ tested |
+| ✅ `cgals2m.m` | ✅ tested |                | ✅ `readoff.m` | ✅ tested |
 | ✅ `vol2restrictedtri.m` | ✅ tested |      | ⭕️ `readsmf.m` | ⭕️ tested |
 | ✅ `surf2volz.m` | ✅ tested |              | ✅ `readtetgen.m` | ✅ tested |
 | ✅ `mesh2mask.m` | ⭕️ tested |              | ✅ `deletemeshfile.m` | ✅ tested |

iso2mesh/core.py

@@ -897,7 +897,7 @@ def cgals2m(v, f, opt, maxvol, **kwargs):
         If it's a scalar, it only specifies radbound.
     maxvol : float
         Target maximum tetrahedral element volume.
-    *args : Additional arguments
+    kwargs : Additional arguments
 
     Returns:
     node : ndarray
@@ -920,7 +920,7 @@ def cgals2m(v, f, opt, maxvol, **kwargs):
     if not isinstance(opt, dict):
         ssize = opt
 
-    if isinstance(opt, dict) and len(opt) == 1:
+    if isinstance(opt, dict):
         ssize = opt.get("radbound", ssize)
         ang = opt.get("angbound", ang)
         approx = opt.get("distbound", approx)

test/run_test.py

@@ -1180,6 +1180,15 @@ def test_cgalv2m(self):
         self.assertEqual(removedupelem(fc[:, :3]).shape[0], 0)
         self.assertAlmostEqual(sum(elemvolume(no[:, :3], el[:, :4])), 0.7455, 3)
 
+    def test_cgals2m(self):
+        node, elem, face = cgals2m(self.no[:, :3], self.fc[:, :3], 3, 50)
+        self.assertAlmostEqual(
+            sum(elemvolume(node[:, :3], face[:, :3])) * 0.001, 2.393432610273361, 3
+        )
+        self.assertAlmostEqual(
+            sum(elemvolume(node[:, :3], elem[:, :4])) * 0.001, 3.817319008857012, 3
+        )
+
     def test_v2s_label(self):
         no, fc, _, _ = v2s(self.mask, 0.5, 0.5)
         self.assertAlmostEqual(
@@ -1265,8 +1274,8 @@ def test_surf2vol_fill(self):
         self.assertAlmostEqual(np.sum(vol.astype(np.float32)) * 0.0001, 3.5536, 2)
 
     def test_s2v(self):
-        vol = s2v(self.no, self.fc, 80, fill=1)
-        self.assertAlmostEqual(np.sum(vol.astype(np.float32)) * 0.0001, 31.4066, 2)
+        vol = s2v(self.no, self.fc, 100, fill=1)
+        self.assertAlmostEqual(np.sum(vol.astype(np.float32)) * 0.0001, 60.827, 2)
 
 
 @unittest.skipIf(