Improve test coverage

Author: ArzhelaR

mesh_model/mesh_analysis/quadmesh_analysis.py

@@ -56,105 +56,6 @@ def set_scores(self) -> None:
                 n.set_score(s)
             i += 1
 
-    def set_geometric_quality(self) -> None:
-        for d_info in self.mesh.active_darts():
-            d_id = d_info[0]
-            d = Dart(self.mesh, d_id)
-            d.set_quality(self.get_dart_geometric_quality(d))
-
-    def get_dart_geometric_quality(self, d: Dart, m=None) -> int:
-        """
-        Calculate the geometric quality of the surrounding of a dart and his twin dart.
-            * quality = -1: boundary dart
-            * quality = 0: convex quad
-            * quality = 1: concave quad
-            * quality = 2: triangular quad
-            * quality = 3: crossed quad
-            * quality = 4: half flat concave quad
-            * quality = 6: flat quad
-            * quality = 5: one adjacent triangular face is flat
-
-        :param d: dart
-        :return: geometric quality
-        """
-        if d.get_beta(2) is None:
-            return -1 # boundary dart
-
-        d2, d1, d11, d21, d211, A, B, C, D = self.mesh.active_triangles(d)
-
-        u1 = np.array([B.x() - A.x(), B.y() - A.y()])  # d vector (unused)
-        u2 = np.array([C.x() - B.x(), C.y() - B.y()])  # vect(BC)
-        u3 = np.array([A.x() - C.x(), A.y() - C.y()])  # vect(CA)
-        u4 = np.array([D.x() - A.x(), D.y() - A.y()])  # vect(AD)
-        u5 = np.array([B.x() - D.x(), B.y() - D.y()])  # vect(DB)
-
-        # Checking for near-zero vectors (close to (0,0))
-        if (np.linalg.norm(u2) < 1e-8 or
-                np.linalg.norm(u3) < 1e-8 or
-                np.linalg.norm(u4) < 1e-8 or
-                np.linalg.norm(u5) < 1e-8):
-            raise ValueError("near zero vector") # Quad invalid because one side is almost zero
-
-        # Calculate cross product at each node
-        cp_A = self.cross_product(-1 * u3, u4)
-        cp_D = self.cross_product(-1 * u4, u5)
-        cp_B = self.cross_product(-1 * u5, u2)
-        cp_C = self.cross_product(-1 * u2, u3)
-
-        # Counts how many angles are oriented clockwise. If the angle is clockwise oriented, signe(cp) return 1
-        sum_cp = self.signe(cp_A) + self.signe(cp_B) + self.signe(cp_C) + self.signe(cp_D)
-
-
-        zero_count = sum(-1e-8 < cp < 1e-8 for cp in [cp_A, cp_B, cp_C, cp_D])
-
-        if zero_count == 0: # not angles of 180° or 360° / no coolinear vectors
-            if sum_cp == 0:
-                return 0 #convexe
-            elif sum_cp == 1:
-                return 1 # concave
-            elif sum_cp == 2:
-                return 3 # crossed
-        elif zero_count == 1:
-            if sum_cp == 0:
-                return 2 # triangular quad
-            elif sum_cp == 1:
-                print(cp_A, cp_B, cp_C, cp_D)
-                return 4 # half_face_flat
-            elif sum_cp == 2:
-                return 5 # half face flat
-        elif zero_count >= 2:
-            return 6 # full flat face
-        else:
-            raise ValueError("Quality configuration doesn't exist")
-
-    def get_adjacent_faces(self, n: Node, d_from: Dart, d_to: Dart) -> list:
-        adj_faces = []
-        d2 = d_from
-        d = None if d2 is None else d_from.get_beta(1)
-        while d != d_to:
-            if d2 is None and d_to is not None:
-                # chercher dans l'autre sens
-                d = d_to
-                adj_faces.append(d.get_face())
-                d1 = d.get_beta(1)
-                d11 = d1.get_beta(1)
-                d = d11.get_beta(2)
-                while d is not None:
-                    adj_faces.append(d.get_face())
-                    d1 = d.get_beta(1)
-                    d11 = d1.get_beta(1)
-                    d = d11.get_beta(2)
-                break
-            elif d2 is None and d_to is None:
-                break
-            elif d2 is not None:
-                d = d2.get_beta(1)
-                adj_faces.append(d.get_face())
-                d2 = d.get_beta(2)
-            else:
-                break
-        return adj_faces
-
 class QuadMeshOldAnalysis(QuadMeshAnalysis):
     """
     Quadmesh old analysis