LoopStructural 1.6.9__py3-none-any.whl → 1.6.10__py3-none-any.whl

LoopStructural/datatypes/_bounding_box.py

@@ -326,7 +326,9 @@ class BoundingBox:
         if iy == -1:
             return self.origin[ix]
 
-        return self.bb[ix,]
+        return self.bb[
+            ix,
+        ]
 
     def __getitem__(self, name):
         if isinstance(name, str):
@@ -429,6 +431,7 @@ class BoundingBox:
             "maximum": self.maximum.tolist(),
             "nsteps": self.nsteps.tolist(),
         }
+
     @classmethod
     def from_dict(cls, data: dict) -> 'BoundingBox':
         """Create a bounding box from a dictionary

LoopStructural/datatypes/_surface.py

@@ -1,5 +1,5 @@
 from dataclasses import dataclass, field
-from typing import Optional
+from typing import Optional, Union
 import numpy as np
 import io
 from LoopStructural.utils import getLogger
@@ -11,6 +11,8 @@ logger = getLogger(__name__)
 class Surface:
     vertices: np.ndarray = field(default_factory=lambda: np.array([[0, 0, 0]]))
     triangles: np.ndarray = field(default_factory=lambda: np.array([[0, 0, 0]]))
+    colour: Optional[Union[str, np.ndarray]] = field(default_factory=lambda: None)
+
     normals: Optional[np.ndarray] = None
     name: str = 'surface'
     values: Optional[np.ndarray] = None

LoopStructural/interpolators/_finite_difference_interpolator.py

@@ -272,12 +272,7 @@ class FiniteDifferenceInterpolator(DiscreteInterpolator):
             idc[inside, :] = gi[node_idx[inside, :]]
             inside = np.logical_and(~np.any(idc == -1, axis=1), inside)
 
-            (
-                vertices,
-                T,
-                elements,
-                inside_,
-            ) = self.support.get_element_gradient_for_location(
+            (vertices, T, elements, inside_,) = self.support.get_element_gradient_for_location(
                 points[inside, : self.support.dimension]
             )
             # normalise constraint vector and scale element matrix by this
@@ -340,12 +335,7 @@ class FiniteDifferenceInterpolator(DiscreteInterpolator):
             # calculate unit vector for node gradients
             # this means we are only constraining direction of grad not the
             # magnitude
-            (
-                vertices,
-                T,
-                elements,
-                inside_,
-            ) = self.support.get_element_gradient_for_location(
+            (vertices, T, elements, inside_,) = self.support.get_element_gradient_for_location(
                 points[inside, : self.support.dimension]
             )
             # T*=np.product(self.support.step_vector)
@@ -432,12 +422,7 @@ class FiniteDifferenceInterpolator(DiscreteInterpolator):
             vectors[norm > 0, :] /= norm[norm > 0, None]
 
             # normalise element vector to unit vector for dot product
-            (
-                vertices,
-                T,
-                elements,
-                inside_,
-            ) = self.support.get_element_gradient_for_location(
+            (vertices, T, elements, inside_,) = self.support.get_element_gradient_for_location(
                 points[inside, : self.support.dimension]
             )
             T[norm > 0, :, :] /= norm[norm > 0, None, None]

LoopStructural/interpolators/supports/_3d_base_structured.py

@@ -162,9 +162,9 @@ class BaseStructuredSupport(BaseSupport):
         length = self.maximum - origin
         length /= self.step_vector
         self._nsteps = np.ceil(length).astype(np.int64)
-        self._nsteps[self._nsteps == 0] = (
-            3  # need to have a minimum of 3 elements to apply the finite difference mask
-        )
+        self._nsteps[
+            self._nsteps == 0
+        ] = 3  # need to have a minimum of 3 elements to apply the finite difference mask
         if np.any(~(self._nsteps > 0)):
             logger.error(
                 f"Cannot resize the interpolation support. The proposed number of steps is {self._nsteps}, these must be all > 0"

LoopStructural/interpolators/supports/_3d_structured_tetra.py

@@ -166,9 +166,9 @@ class TetMesh(BaseStructuredSupport):
         shared_face_index[:] = -1
         shared_face_index[row.reshape(-1, 3)[:, 0], :] = col.reshape(-1, 3)
 
-        self.shared_elements[np.arange(self.shared_element_relationships.shape[0]), :] = (
-            shared_face_index
-        )
+        self.shared_elements[
+            np.arange(self.shared_element_relationships.shape[0]), :
+        ] = shared_face_index
         # resize
         self.shared_elements = self.shared_elements[: len(self.shared_element_relationships), :]
 

LoopStructural/interpolators/supports/_3d_unstructured_tetra.py

@@ -173,9 +173,9 @@ class UnStructuredTetMesh(BaseSupport):
         shared_face_index[:] = -1
         shared_face_index[row.reshape(-1, 3)[:, 0], :] = col.reshape(-1, 3)
 
-        self.shared_elements[np.arange(self.shared_element_relationships.shape[0]), :] = (
-            shared_face_index
-        )
+        self.shared_elements[
+            np.arange(self.shared_element_relationships.shape[0]), :
+        ] = shared_face_index
         # resize
         self.shared_elements = self.shared_elements[: len(self.shared_element_relationships), :]
         # flag = np.zeros(self.elements.shape[0])

LoopStructural/interpolators/supports/_face_table.py

@@ -63,8 +63,8 @@ def _init_face_table(grid):
     shared_face_index = np.zeros((shared_faces.shape[0], grid.dimension), dtype=int)
     shared_face_index[:] = -1
     shared_face_index[row.reshape(-1, grid.dimension)[:, 0], :] = col.reshape(-1, grid.dimension)
-    grid._shared_elements[np.arange(grid.shared_element_relationships.shape[0]), :] = (
-        shared_face_index
-    )
+    grid._shared_elements[
+        np.arange(grid.shared_element_relationships.shape[0]), :
+    ] = shared_face_index
     # resize
     grid._shared_elements = grid.shared_elements[: len(grid.shared_element_relationships), :]

LoopStructural/modelling/core/geological_model.py

@@ -598,10 +598,12 @@ class GeologicalModel:
                 * self._data.loc[mask, "polarity"].to_numpy()[:, None]
             )
             self._data.drop(["strike", "dip"], axis=1, inplace=True)
-        self._data[['X', 'Y', 'Z', 'val', 'nx', 'ny', 'nz', 'gx', 'gy', 'gz', 'tx', 'ty', 'tz']] = (
-            self._data[
-                ['X', 'Y', 'Z', 'val', 'nx', 'ny', 'nz', 'gx', 'gy', 'gz', 'tx', 'ty', 'tz']
-            ].astype(float)
+        self._data[
+            ['X', 'Y', 'Z', 'val', 'nx', 'ny', 'nz', 'gx', 'gy', 'gz', 'tx', 'ty', 'tz']
+        ] = self._data[
+            ['X', 'Y', 'Z', 'val', 'nx', 'ny', 'nz', 'gx', 'gy', 'gz', 'tx', 'ty', 'tz']
+        ].astype(
+            float
         )
 
     def set_model_data(self, data):
@@ -1293,6 +1295,7 @@ class GeologicalModel:
         fault_trace_anisotropy=0.0,
         fault_dip=90,
         fault_dip_anisotropy=0.0,
+        fault_pitch=None,
         **kwargs,
     ):
         """
@@ -1379,6 +1382,7 @@ class GeologicalModel:
             major_axis = major_axis / self.scale_factor
         if intermediate_axis:
             intermediate_axis = intermediate_axis / self.scale_factor
+        print(fault_dip, fault_normal_vector)
         fault_frame_builder.create_data_from_geometry(
             fault_frame_data=fault_frame_data,
             fault_center=fault_center,
@@ -1393,6 +1397,7 @@ class GeologicalModel:
             fault_trace_anisotropy=fault_trace_anisotropy,
             fault_dip=fault_dip,
             fault_dip_anisotropy=fault_dip_anisotropy,
+            fault_pitch=fault_pitch,
         )
         if "force_mesh_geometry" not in kwargs:
             fault_frame_builder.set_mesh_geometry(kwargs.get("fault_buffer", 0.2), 0)
@@ -1798,7 +1803,10 @@ class GeologicalModel:
             values = values.loc[~np.logical_or(values == np.inf, values == -np.inf)]
             surfaces.extend(
                 self.get_feature_by_name(u).surfaces(
-                    values.to_list(), self.bounding_box, name=names.loc[values.index].to_list()
+                    values.to_list(),
+                    self.bounding_box,
+                    name=names.loc[values.index].to_list(),
+                    colours=unit_table.loc[unit_table['feature_name'] == u, 'colour'].tolist(),
                 )
             )
 

LoopStructural/modelling/features/_base_geological_feature.py

@@ -274,6 +274,7 @@ class BaseFeature(metaclass=ABCMeta):
         value: Union[float, int, List[Union[float, int]]],
         bounding_box=None,
         name: Optional[Union[List[str], str]] = None,
+        colours: Optional[Union[str, np.ndarray]] = None,
     ) -> surface_list:
         """Find the surfaces of the geological feature at a given value
 
@@ -305,7 +306,7 @@ class BaseFeature(metaclass=ABCMeta):
             isosurfacer = LoopIsosurfacer(bounding_box, callable=callable)
             if name is None and self.name is not None:
                 name = self.name
-            surfaces = isosurfacer.fit(value, name)
+            surfaces = isosurfacer.fit(value, name, colours=colours)
         except Exception as e:
             logger.error(f"Failed to create surface for {self.name} at value {value}")
             logger.error(e)
@@ -333,9 +334,7 @@ class BaseFeature(metaclass=ABCMeta):
                 raise ValueError("Must specify bounding box")
             bounding_box = self.model.bounding_box
         grid = bounding_box.structured_grid(name=self.name)
-        value = self.evaluate_value(
-            bounding_box.regular_grid(local=False, order='F')
-        )
+        value = self.evaluate_value(bounding_box.regular_grid(local=False, order='F'))
         if self.model is not None:
 
             value = self.evaluate_value(

LoopStructural/modelling/features/builders/_fault_builder.py

@@ -91,6 +91,7 @@ class FaultBuilder(StructuralFrameBuilder):
         fault_trace_anisotropy=1.0,
         fault_dip=90,
         fault_dip_anisotropy=1.0,
+        fault_pitch=None,
     ):
         """Generate the required data for building a fault frame for a fault with the
         specified parameters
@@ -204,6 +205,7 @@ class FaultBuilder(StructuralFrameBuilder):
                 fault_frame_data["coord"] == 1, ~np.isnan(fault_frame_data["gz"])
             )
             fault_slip_data = fault_frame_data.loc[slip_mask, ["gx", "gy", "gz"]]
+
             if len(fault_slip_data) == 0:
                 logger.warning(
                     "There is no slip vector data for the fault, using vertical slip vector\n\
@@ -211,6 +213,11 @@ class FaultBuilder(StructuralFrameBuilder):
                 )
                 strike_vector, dip_vector = get_vectors(fault_normal_vector[None, :])
                 fault_slip_vector = dip_vector[:, 0]
+                if fault_pitch is not None:
+                    print('using pitch')
+                    rotm = rotation(fault_normal_vector[None,:],[fault_pitch])
+                    print(rotm.shape,fault_slip_vector.shape)
+                    fault_slip_vector = np.einsum("ijk,k->ij", rotm, fault_slip_vector)[0,:]
                 logger.info(f"Estimated fault slip vector: {fault_slip_vector}")
             else:
                 fault_slip_vector = fault_slip_data.mean(axis=0).to_numpy()

LoopStructural/modelling/input/process_data.py

@@ -299,9 +299,9 @@ class ProcessInputData:
             pts = self.fault_locations.loc[
                 self.fault_locations["feature_name"] == fname, ["X", "Y", "Z"]
             ]
-            fault_properties.loc[fname, ["centreEasting", "centreNorthing", "centreAltitude"]] = (
-                np.nanmean(pts, axis=0)
-            )
+            fault_properties.loc[
+                fname, ["centreEasting", "centreNorthing", "centreAltitude"]
+            ] = np.nanmean(pts, axis=0)
         if (
             "avgNormalEasting" not in fault_properties.columns
             or "avgNormalNorthing" not in fault_properties.columns
@@ -449,9 +449,9 @@ class ProcessInputData:
         for _name, sg in self.stratigraphic_order:
             value = 0.0  # reset for each supergroup
             if sg[0] not in self.thicknesses or self.thicknesses[sg[0]] <= 0:
-                self.thicknesses[sg[0]] = (
-                    np.inf
-                )  # make the top unit infinite as it should extend to the top of the model
+                self.thicknesses[
+                    sg[0]
+                ] = np.inf  # make the top unit infinite as it should extend to the top of the model
             for g in reversed(
                 sg[:-1]
             ):  # don't add the last unit as we never see the base of this unit.

LoopStructural/utils/_surface.py

@@ -64,6 +64,7 @@ class LoopIsosurfacer:
         values: Optional[Union[list, int, float]],
         name: Optional[Union[List[str], str]] = None,
         local=False,
+        colours: Optional[List] = None,
     ) -> surface_list:
         """Extract isosurfaces from the interpolator
 
@@ -120,7 +121,9 @@ class LoopIsosurfacer:
             names = [name] * len(isovalues)
         if isinstance(name, list):
             names = name
-        for name, isovalue in zip(names, isovalues):
+        if colours is None:
+            colours = [None] * len(isovalues)
+        for name, isovalue, colour in zip(names, isovalues, colours):
             try:
                 step_vector = (self.bounding_box.maximum - self.bounding_box.origin) / (
                     np.array(self.bounding_box.nsteps) - 1
@@ -150,6 +153,7 @@ class LoopIsosurfacer:
                     normals=normals,
                     name=f"{name}_{isovalue}",
                     values=values,
+                    colour=colour,
                 )
             )
         return surfaces

LoopStructural/version.py

@@ -1 +1 @@
-__version__ = "1.6.9"
+__version__ = "1.6.10"

{loopstructural-1.6.9.dist-info → loopstructural-1.6.10.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: LoopStructural
-Version: 1.6.9
+Version: 1.6.10
 Summary: 3D geological modelling
 Author-email: Lachlan Grose <lachlan.grose@monash.edu>
 License: MIT

{loopstructural-1.6.9.dist-info → loopstructural-1.6.10.dist-info}/RECORD

@@ -1,5 +1,5 @@
 LoopStructural/__init__.py,sha256=fg_Vm1aMDYIf_CffTFopLsTx21u6deLaI7JMVpRYdOI,1378
-LoopStructural/version.py,sha256=cqdgnlJ01g3Nm9hqIB0XVUpI__WAWqmD626VQ6wd9sA,22
+LoopStructural/version.py,sha256=sWIpFBcO73DLgomFrOsnrVoNZ_1WQXvqT1hgrs4ja8E,23
 LoopStructural/datasets/__init__.py,sha256=ylb7fzJU_DyQ73LlwQos7VamqkDSGITbbnoKg7KAOmE,677
 LoopStructural/datasets/_base.py,sha256=FB_D5ybBYHoaNbycdkpZcRffzjrrL1xp9X0k-pyob9Y,7618
 LoopStructural/datasets/_example_models.py,sha256=Zg33IeUyh4C-lC0DRMLqCDP2IrX8L-gNV1WxJwBGjzM,113
@@ -29,10 +29,10 @@ LoopStructural/datasets/data/geological_map_data/stratigraphic_order.csv,sha256=
 LoopStructural/datasets/data/geological_map_data/stratigraphic_orientations.csv,sha256=RysyqUAIjY6iIDUfTh11n9QUQWXB_qxKnZeN_DqNzlY,26745
 LoopStructural/datasets/data/geological_map_data/stratigraphic_thickness.csv,sha256=pnSmG-wL8-kxuoHo_pgpJrfTmsZOzc8L0vxpBRh3r8A,355
 LoopStructural/datatypes/__init__.py,sha256=lVg64DnynMm58qvYTjLrcyWH7vk2ngr9JGMo5FaiALI,160
-LoopStructural/datatypes/_bounding_box.py,sha256=6yuqWqUzlWF1k0I09FgaBnb5bDitMQaXjHbwdrIizok,18352
+LoopStructural/datatypes/_bounding_box.py,sha256=NwERtFHvrlQqsUNFZm4JgPEjCrxs31GlaDC7DRG6bO4,18375
 LoopStructural/datatypes/_point.py,sha256=sCZiWRdnbIH_P9b52btNrotrttWlE-e4Mr2BBwfAAZY,7710
 LoopStructural/datatypes/_structured_grid.py,sha256=mc-UM1Gh_BjHFItuPE4FF5wvGzJnSqF2MTx_xvrwcTk,5088
-LoopStructural/datatypes/_surface.py,sha256=IgXKcwZw3mF2DIImO7LT-HUapqABa8LU0WsxXTiqEHo,6562
+LoopStructural/datatypes/_surface.py,sha256=V3HSvPO-awPMlIG22A1ttUCn1BUcdN-29CiDj6CtF-w,6653
 LoopStructural/export/exporters.py,sha256=BniZu-PqQvHqCU6GIuJQ5FPzI9Dx_T6rI8EW1pykois,17209
 LoopStructural/export/file_formats.py,sha256=0xKyYSW4Jv_4jsXwusg-WO6PNUhZKd6HdWSqGSaPve8,232
 LoopStructural/export/geoh5.py,sha256=jLFKC5EB0azT3PgJPtkJzi3_CG28RLgP2FuENAGCQMI,4313
@@ -43,7 +43,7 @@ LoopStructural/interpolators/_api.py,sha256=EC4ogG2uPq-z_pgNGd_eTieTl92eaZ-rjyoF
 LoopStructural/interpolators/_builders.py,sha256=B49KsxB8RRN6IHDfGT43nXWe_Av1SVVT8vm2Nh1oEiQ,6758
 LoopStructural/interpolators/_discrete_fold_interpolator.py,sha256=eDe0R1lcQ0AuMcv7zlpu5c-soCv7AybIqQAuN2vFE3M,6542
 LoopStructural/interpolators/_discrete_interpolator.py,sha256=i_joZ8HOf_s6Q2L8gHFnhkdtgyED1SjATxRsRd1HxRU,26038
-LoopStructural/interpolators/_finite_difference_interpolator.py,sha256=mZ89FWQZ5RbzhL9UYH8VzWME-28dy331KXtYtpqepHo,18351
+LoopStructural/interpolators/_finite_difference_interpolator.py,sha256=SUU_7dmfhKFh3xfwG_w-yCAzWDMoVNC0qHxDW17ajgc,18117
 LoopStructural/interpolators/_geological_interpolator.py,sha256=hcQuyv1zYakJ7mcDFlLj-YarjnMQvlP6pVbK1KuxBWs,11195
 LoopStructural/interpolators/_interpolator_builder.py,sha256=Z8bhmco5aSQX19A8It2SB_rG61wnlyshWfp3ivm8rU0,4586
 LoopStructural/interpolators/_interpolator_factory.py,sha256=fbjebXSe5IgTol1tnBlnsw9gD426v-TGkX3gquIg7LI,2782
@@ -57,22 +57,22 @@ LoopStructural/interpolators/supports/_2d_p1_unstructured.py,sha256=okcy56nyjued
 LoopStructural/interpolators/supports/_2d_p2_unstructured.py,sha256=TeBVtT1PMV7CKzmnFZke37acMoFxouer20cskS7pVoE,10422
 LoopStructural/interpolators/supports/_2d_structured_grid.py,sha256=Pt9fiXyTS-RTd3mxXr3EUQfB6DhKChHQ5zbWub54nW0,16347
 LoopStructural/interpolators/supports/_2d_structured_tetra.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-LoopStructural/interpolators/supports/_3d_base_structured.py,sha256=PYIkgWCTHk0OHn-T9wpF4ZFYnKyvXNZgUyqO7YcoFjU,16481
+LoopStructural/interpolators/supports/_3d_base_structured.py,sha256=XsIGQcv0u43ld90lL7K8f3kwmMJuiSnHem-0LrrJZlY,16479
 LoopStructural/interpolators/supports/_3d_p2_tetra.py,sha256=CqGVJRUMxbPQZDhhopNt_s9gVhMqh4YbjQyDZonoyxc,11574
 LoopStructural/interpolators/supports/_3d_structured_grid.py,sha256=x9NoZRsl58iowcObavgb0nY_C335BmcIYgec9REsFpU,17366
-LoopStructural/interpolators/supports/_3d_structured_tetra.py,sha256=5zUNtvEXDvbCHZCu6Fz9WjGbnrMaq-sYJqNUufyLcq8,26505
-LoopStructural/interpolators/supports/_3d_unstructured_tetra.py,sha256=_peXMTMxctuWNOL74AHxzw0b_1sP5glvbJigIvIkK9I,23867
+LoopStructural/interpolators/supports/_3d_structured_tetra.py,sha256=8poHH0UIMWLpDPHZZG39dSSYiu7vBz-GBHqRpGovFCU,26503
+LoopStructural/interpolators/supports/_3d_unstructured_tetra.py,sha256=TjPuvmm51Ct2zxNeeYr2LjBT90PC1CW5rDV1CIBvBpg,23865
 LoopStructural/interpolators/supports/__init__.py,sha256=V0JjixoBIUZVAo5MmqARR67xDOoQwnb4G3SXeOMRSyQ,1603
 LoopStructural/interpolators/supports/_aabb.py,sha256=Z-kH_u6c6izak0aHG3Uo14PEKQeZmYlevLDC32Q06xk,3208
 LoopStructural/interpolators/supports/_base_support.py,sha256=pYzsmeBu4kLaD9ZKsz_dfjVpfuAd00xENqOQC9Xw5QY,2501
-LoopStructural/interpolators/supports/_face_table.py,sha256=Hyj4Io63NkPRN8ab9uDHyec-2Kb8BLY_xBF6STNlvBw,3095
+LoopStructural/interpolators/supports/_face_table.py,sha256=CjBcGeIdIwPo0IV_-uT6SK-Re0qSeSzPAtuoztmvQGg,3093
 LoopStructural/interpolators/supports/_support_factory.py,sha256=XNAxnr-JS3KEhdsoZeJ-VaLTJwlvxgBuRMCqYrCDW18,1485
 LoopStructural/modelling/__init__.py,sha256=oW7dz6c8K1A0VcW7-mVcyqcENUrtybCb3eVUNXFvMfA,656
 LoopStructural/modelling/core/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-LoopStructural/modelling/core/geological_model.py,sha256=TILvSGV95WtPJArIorYIu27pH3Faf7oTWCh2xj6qqq0,66172
+LoopStructural/modelling/core/geological_model.py,sha256=bhzdOIbFAHVCnUwJe32GEN2nf1D1cT4KH-JXirZpDfY,66430
 LoopStructural/modelling/features/__init__.py,sha256=Vf-qd5EDBtJ1DpuXXyCcw2-wf6LWPRW5wzxDEO3vOc8,939
 LoopStructural/modelling/features/_analytical_feature.py,sha256=U_g86LgQhYY2359rdsDqpvziYwqrWkc5EdvhJARiUWo,3597
-LoopStructural/modelling/features/_base_geological_feature.py,sha256=z7z6v8MwYRVH5OWOi1vJqTt8vVOwrFysryI71L4cj-4,12122
+LoopStructural/modelling/features/_base_geological_feature.py,sha256=985aHq8z5riS7uS04mGBPd8UcPXFv0l_fPAdlE6Ydcc,12175
 LoopStructural/modelling/features/_cross_product_geological_feature.py,sha256=GIyCHUdE6F-bse2e4puG9V2f7qRtDVfby5PRe2BboD4,3021
 LoopStructural/modelling/features/_geological_feature.py,sha256=u6pbKj9BujX1Ijj5eVdhwGDNjrIAI16CpiAn5n8g3RY,11279
 LoopStructural/modelling/features/_lambda_geological_feature.py,sha256=ieZiiMd8HCeVEIl3B406tsysT7iUEGcEUJXALjhjkeM,2924
@@ -82,7 +82,7 @@ LoopStructural/modelling/features/_structural_frame.py,sha256=e3QmNHLwuZc5PX3rLa
 LoopStructural/modelling/features/_unconformity_feature.py,sha256=2Bx0BI38YLdcNvDWuP9E1pKFN4orEUq9aC8b5xG1UVk,2362
 LoopStructural/modelling/features/builders/__init__.py,sha256=Gqld1C-PcaXfJ8vpkWMDCmehmd3hZNYQk1knPtl59Bk,266
 LoopStructural/modelling/features/builders/_base_builder.py,sha256=N3txGC98V08A8-k2TLdoIWgWLfblZ91kaTvciPq_QVM,3750
-LoopStructural/modelling/features/builders/_fault_builder.py,sha256=EyDLyovrnxqSk11zWFnViok7lXMt_E-pG2QqgRgvOLk,25653
+LoopStructural/modelling/features/builders/_fault_builder.py,sha256=3rNfRmMcMCrXMyc_hjltU30wTh041NWD6rShNVI2lCE,25999
 LoopStructural/modelling/features/builders/_folded_feature_builder.py,sha256=1_0BVTzcvmFl6K3_lX-jF0tiMFPmS8j6vPeSLn9MbrE,6607
 LoopStructural/modelling/features/builders/_geological_feature_builder.py,sha256=jn2BiZlzXyWl0_TrsajpFR2wegGOpbuO5yFu2FamuYA,22014
 LoopStructural/modelling/features/builders/_structural_frame_builder.py,sha256=ms3-fuFpDEarjzYU5W499TquOIlTwHPUibVxIypfmWY,8019
@@ -103,7 +103,7 @@ LoopStructural/modelling/features/fold/fold_function/_trigo_fold_rotation_angle.
 LoopStructural/modelling/input/__init__.py,sha256=HhJM3V5b-8_64LiRbF3Bd1pjWhJlcknxMSMPRrqZ0-I,153
 LoopStructural/modelling/input/fault_network.py,sha256=0uxl7lOySdhMhNXoiOkuiHIXqAz1Ls0j-W65cmdQoP8,2348
 LoopStructural/modelling/input/map2loop_processor.py,sha256=T7Fgqd7FNJWylLKvfIniRZBMRMeAoP8iU330-WYU8Fg,7031
-LoopStructural/modelling/input/process_data.py,sha256=9xWZtFWbE5NjJzULSZDJt-qccV0GuZ-oNMi3nfxPT-M,26293
+LoopStructural/modelling/input/process_data.py,sha256=iEq5bj4xqhYUfTd8cBIr03_PIeKBYb4mmKTpRqEbmpU,26289
 LoopStructural/modelling/input/project_file.py,sha256=WhJkMfDK9uE7MK7HK-YK6ZOBAdwLX5P7ThZgXj444Eg,4604
 LoopStructural/modelling/intrusions/__init__.py,sha256=EpZK3cHJwGQhPUYIwKCKu8vkNdt_nOgWF0zfhiqDYDA,712
 LoopStructural/modelling/intrusions/geom_conceptual_models.py,sha256=jwTlhYySUj7z4DEnJoi4AINZB_N3-SW6ONRFL66OsW0,3665
@@ -113,7 +113,7 @@ LoopStructural/modelling/intrusions/intrusion_feature.py,sha256=ESjtikHFJQzUnowb
 LoopStructural/modelling/intrusions/intrusion_frame_builder.py,sha256=Q1TPHxREcrO7Rw71nUfACZHfYnISLjqlgkUNTPT324k,40143
 LoopStructural/modelling/intrusions/intrusion_support_functions.py,sha256=wodakheMD62WJyoKnyX8UO-C1pje0I-5kHQEoDqShzo,13951
 LoopStructural/utils/__init__.py,sha256=OJqNSu40SYJeC26IhoBBXDqQOogWjMGA-YokKVRrwMs,924
-LoopStructural/utils/_surface.py,sha256=qXn0WHiMlYx-e_eRSoRcBIIrAlM8WUsi0cthiUQ_ptM,5974
+LoopStructural/utils/_surface.py,sha256=vTVIkhjYC85CPyirtA4aMPi295fpo40QHzYOqeUS354,6140
 LoopStructural/utils/_transformation.py,sha256=peuLPH3BJ5DxnPbOuNKcqK4eXhAXdbT540L1OIsO3v0,5404
 LoopStructural/utils/colours.py,sha256=-KRf1MXKx4L8TXnwyiunmKAX4tfy0qG68fRadyfn_bM,1163
 LoopStructural/utils/config.py,sha256=ITGOtZTo2_QBwXkG_0AFANfE90J9siCXLzxypVmg9QA,414
@@ -129,8 +129,8 @@ LoopStructural/utils/regions.py,sha256=LvcOCPudF4u95-GKBOZqXVxOEcR3cOFgFpcs5x43s
 LoopStructural/utils/typing.py,sha256=29uVSTZdzXXH-jdlaYyBWZ1gQ2-nlZ2-XoVgG_PXNFY,157
 LoopStructural/utils/utils.py,sha256=2Z4zVE6G752-SPmM29zebk82bROJxEwi_YiiJjcVED4,2438
 LoopStructural/visualisation/__init__.py,sha256=5BDgKor8-ae6DrS7IZybJ3Wq_pTnCchxuY4EgzA7v1M,318
-loopstructural-1.6.9.dist-info/licenses/LICENSE,sha256=ZqGeNFOgmYevj7Ld7Q-kR4lAxWXuBRUdUmPC6XM_py8,1071
-loopstructural-1.6.9.dist-info/METADATA,sha256=em1vuQTzofTNruIoPIFthFCELOwH0KwlpUZpZnXNUDo,6453
-loopstructural-1.6.9.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
-loopstructural-1.6.9.dist-info/top_level.txt,sha256=QtQErKzYHfg6ddxTQ1NyaTxXBVM6qAqrM_vxEPyXZLg,15
-loopstructural-1.6.9.dist-info/RECORD,,
+loopstructural-1.6.10.dist-info/licenses/LICENSE,sha256=ZqGeNFOgmYevj7Ld7Q-kR4lAxWXuBRUdUmPC6XM_py8,1071
+loopstructural-1.6.10.dist-info/METADATA,sha256=SwyCi3Rxn8nYQVe4GcIwYo6LkYXQyxDRjH6a_idDpow,6454
+loopstructural-1.6.10.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
+loopstructural-1.6.10.dist-info/top_level.txt,sha256=QtQErKzYHfg6ddxTQ1NyaTxXBVM6qAqrM_vxEPyXZLg,15
+loopstructural-1.6.10.dist-info/RECORD,,