cleanup

Author: benk-mira

simpeg_drivers/plate_simulation/models/parametric.py

@@ -47,123 +47,6 @@ def mask(self, mesh: Octree) -> np.ndarray:
         """
 
 
-# class Plate(Parametric):
-#     """
-#     Define a rotated rectangular block in 3D space
-#
-#     :param params: Parameters describing the plate.
-#     :param surface: Surface object representing the plate.
-#     """
-#
-#     def __init__(
-#         self,
-#         params: PlateOptions,
-#         workspace: Workspace | None = None,
-#     ):
-#         self.params = params
-#         self._workspace = workspace
-#         super().__init__(self._create_surface())
-#
-#     def _create_surface(self) -> Surface:
-#         """
-#         Create a surface object from a plate object.
-#
-#         :param workspace: Workspace object to create the surface in.
-#         :param out_group: Output group to store the surface.
-#         """
-#         with fetch_active_workspace(self.workspace) as ws:
-#             surface = Surface.create(
-#                 ws,
-#                 vertices=self.vertices,
-#                 cells=self.triangles,
-#                 name=self.params.name,
-#             )
-#
-#         return surface
-#
-#     @property
-#     def surface(self):
-#         """
-#         A surface object representing the plate.
-#         """
-#         return self._surface
-#
-#     @property
-#     def triangles(self) -> np.ndarray:
-#         """Triangulation of the block."""
-#         return np.vstack(
-#             [
-#                 [0, 2, 1],
-#                 [1, 2, 3],
-#                 [0, 1, 4],
-#                 [4, 1, 5],
-#                 [1, 3, 5],
-#                 [5, 3, 7],
-#                 [2, 6, 3],
-#                 [3, 6, 7],
-#                 [0, 4, 2],
-#                 [2, 4, 6],
-#                 [4, 5, 6],
-#                 [6, 5, 7],
-#             ]
-#         )
-#
-#     @property
-#     def vertices(self) -> np.ndarray:
-#         """Vertices for triangulation of a rectangular prism in 3D space."""
-#
-#         u_1 = self.params.geometry.origin[0] - (
-#             self.params.geometry.strike_length / 2.0
-#         )
-#         u_2 = self.params.geometry.origin[0] + (
-#             self.params.geometry.strike_length / 2.0
-#         )
-#         v_1 = self.params.geometry.origin[1] - (self.params.geometry.dip_length / 2.0)
-#         v_2 = self.params.geometry.origin[1] + (self.params.geometry.dip_length / 2.0)
-#         w_1 = self.params.geometry.origin[2] - (self.params.geometry.width / 2.0)
-#         w_2 = self.params.geometry.origin[2] + (self.params.geometry.width / 2.0)
-#
-#         vertices = np.array(
-#             [
-#                 [u_1, v_1, w_1],
-#                 [u_2, v_1, w_1],
-#                 [u_1, v_2, w_1],
-#                 [u_2, v_2, w_1],
-#                 [u_1, v_1, w_2],
-#                 [u_2, v_1, w_2],
-#                 [u_1, v_2, w_2],
-#                 [u_2, v_2, w_2],
-#             ]
-#         )
-#
-#         return self._rotate(vertices)
-#
-#     @property
-#     def workspace(self) -> Workspace:
-#         if self._workspace is None:
-#             self._workspace = Workspace()
-#
-#         return self._workspace
-#
-#     def _rotate(self, vertices: np.ndarray) -> np.ndarray:
-#         """Rotate vertices and adjust for reference point."""
-#         theta = -1 * self.params.geometry.direction
-#         phi = -1 * self.params.geometry.dip
-#         rotated_vertices = rotate_xyz(vertices, self.params.geometry.origin, theta, phi)
-#
-#         return rotated_vertices
-#
-#     def mask(self, mesh: Octree) -> np.ndarray:
-#         rotations = [
-#             z_rotation_matrix(np.deg2rad(self.params.geometry.direction)),
-#             x_rotation_matrix(np.deg2rad(self.params.geometry.dip)),
-#         ]
-#         rotated_centers = rotate_points(
-#             mesh.centroids, origin=self.params.geometry.origin, rotations=rotations
-#         )
-#         return inside_plate(rotated_centers, self.params.geometry)
-
-
 class Body(Parametric):
     """
     Represents a closed surface in the model.