LoopStructural 1.6.3__py3-none-any.whl → 1.6.5__py3-none-any.whl

LoopStructural/datatypes/_point.py

@@ -123,12 +123,15 @@ class VectorPoints:
     def from_dict(self, d):
         return VectorPoints(d['locations'], d['vectors'], d['name'], d.get('properties', None))
 
-    def vtk(self, geom='arrow', scale=1.0, scale_function=None, tolerance=0.05):
+    def vtk(self, geom='arrow', scale=1.0, scale_function=None, normalise=True, tolerance=0.05):
         import pyvista as pv
 
         vectors = np.copy(self.vectors)
+        if normalise:
+            norm = np.linalg.norm(vectors, axis=1)
+            vectors[norm > 0, :] /= norm[norm > 0][:, None]
         if scale_function is not None:
-            vectors /= np.linalg.norm(vectors, axis=1)[:, None]
+            # vectors /= np.linalg.norm(vectors, axis=1)[:, None]
             vectors *= scale_function(self.locations)[:, None]
         points = pv.PolyData(self.locations)
         points.point_data.set_vectors(vectors, 'vectors')

LoopStructural/interpolators/supports/_3d_structured_grid.py

@@ -342,9 +342,6 @@ class StructuredGrid(BaseStructuredSupport):
         idc, inside = self.position_to_cell_corners(evaluation_points)
         T = np.zeros((idc.shape[0], 3, 8))
         T[inside, :, :] = self.get_element_gradient_for_location(evaluation_points[inside, :])[1]
-        # indices = np.array([self.position_to_cell_index(evaluation_points)])
-        # idc = self.global_indicies(indices.swapaxes(0,1))
-        # print(idc)
         if np.max(idc[inside, :]) > property_array.shape[0]:
             cix, ciy, ciz = self.position_to_cell_index(evaluation_points)
             if not np.all(cix[inside] < self.nsteps_cells[0]):

LoopStructural/modelling/__init__.py

@@ -25,5 +25,3 @@ except (LoopImportError, ImportError):
     logger.warning(
         "Cannot use LoopProjectfileProcessor: Loop project file cannot be imported, try installing LoopProjectFile"
     )
-# from LoopStructural.modelling.features import StructuralFrame
-# from LoopStructural.modelling.features.fault import FaultSegment

LoopStructural/modelling/core/geological_model.py

@@ -649,7 +649,7 @@ class GeologicalModel:
             for g in stratigraphic_column.keys():
                 for u in stratigraphic_column[g].keys():
                     stratigraphic_column[g][u]["colour"] = cmap_colours[ci, :]
-
+                    ci += 1
         self.stratigraphic_column = stratigraphic_column
 
     def create_and_add_foliation(
@@ -1809,7 +1809,7 @@ class GeologicalModel:
         grid = self.bounding_box.structured_grid(name=name)
 
         grid.cell_properties['stratigraphy'] = self.evaluate_model(
-            self.bounding_box.cell_centers(), scale=False
+            self.rescale(self.bounding_box.cell_centers())
         )
         return grid, self.stratigraphic_ids()
 

LoopStructural/modelling/features/_base_geological_feature.py

@@ -291,11 +291,24 @@ class BaseFeature(metaclass=ABCMeta):
             if self.model is None:
                 raise ValueError("Must specify bounding box")
             bounding_box = self.model.bounding_box
-        callable = lambda xyz: self.evaluate_value(self.model.scale(xyz))
-        isosurfacer = LoopIsosurfacer(bounding_box, callable=callable)
-        if name is None and self.name is not None:
-            name = self.name
-        return isosurfacer.fit(value, name)
+        regions = self.regions
+        try:
+            self.regions = [
+                r for r in self.regions if r.name != self.name and r.parent.name != self.name
+            ]
+            callable = lambda xyz: self.evaluate_value(self.model.scale(xyz))
+            isosurfacer = LoopIsosurfacer(bounding_box, callable=callable)
+            if name is None and self.name is not None:
+                name = self.name
+            surfaces = isosurfacer.fit(value, name)
+        except Exception as e:
+            logger.error(f"Failed to create surface for {self.name} at value {value}")
+            logger.error(e)
+            surfaces = []
+        finally:
+            self.regions = regions
+
+        return surfaces
 
     def scalar_field(self, bounding_box=None):
         """Create a scalar field for the feature
@@ -341,10 +354,10 @@ class BaseFeature(metaclass=ABCMeta):
             if self.model is None:
                 raise ValueError("Must specify bounding box")
             bounding_box = self.model.bounding_box
-        grid = bounding_box.vtk()
-        points = grid.points
+        points = bounding_box.cell_centers()
         value = self.evaluate_gradient(points)
-
+        if self.model is not None:
+            points = self.model.rescale(points)
         return VectorPoints(points, value, self.name)
 
     @abstractmethod

LoopStructural/modelling/features/_structural_frame.py

@@ -115,7 +115,7 @@ class StructuralFrame(BaseFeature):
         """
         return self.features[i]
 
-    def evaluate_value(self, evaluation_points, ignore_regions=False):
+    def evaluate_value(self, pos, ignore_regions=False):
         """
         Evaluate the value of the structural frame for the points.
         Can optionally only evaluate one coordinate
@@ -129,14 +129,14 @@ class StructuralFrame(BaseFeature):
         -------
 
         """
-        v = np.zeros(evaluation_points.shape)  # create new 3d array of correct length
+        v = np.zeros(pos.shape)  # create new 3d array of correct length
         v[:] = np.nan
-        v[:, 0] = self.features[0].evaluate_value(evaluation_points, ignore_regions=ignore_regions)
-        v[:, 1] = self.features[1].evaluate_value(evaluation_points, ignore_regions=ignore_regions)
-        v[:, 2] = self.features[2].evaluate_value(evaluation_points, ignore_regions=ignore_regions)
+        v[:, 0] = self.features[0].evaluate_value(pos, ignore_regions=ignore_regions)
+        v[:, 1] = self.features[1].evaluate_value(pos, ignore_regions=ignore_regions)
+        v[:, 2] = self.features[2].evaluate_value(pos, ignore_regions=ignore_regions)
         return v
 
-    def evaluate_gradient(self, evaluation_points, i=None, ignore_regions=False):
+    def evaluate_gradient(self, pos, i=None, ignore_regions=False):
         """
         Evaluate the gradient of the structural frame.
         Can optionally only evaluate the ith coordinate
@@ -151,19 +151,11 @@ class StructuralFrame(BaseFeature):
 
         """
         if i is not None:
-            return self.features[i].support.evaluate_gradient(
-                evaluation_points, ignore_regions=ignore_regions
-            )
+            return self.features[i].support.evaluate_gradient(pos, ignore_regions=ignore_regions)
         return (
-            self.features[0].support.evaluate_gradient(
-                evaluation_points, ignore_regions=ignore_regions
-            ),
-            self.features[1].support.evaluate_gradient(
-                evaluation_points, ignore_regions=ignore_regions
-            ),
-            self.features[2].support.evaluate_gradient(
-                evaluation_points, ignore_regions=ignore_regions
-            ),
+            self.features[0].support.evaluate_gradient(pos, ignore_regions=ignore_regions),
+            self.features[1].support.evaluate_gradient(pos, ignore_regions=ignore_regions),
+            self.features[2].support.evaluate_gradient(pos, ignore_regions=ignore_regions),
         )
 
     def get_data(self, value_map: Optional[dict] = None) -> List[Union[ValuePoints, VectorPoints]]:

LoopStructural/modelling/features/_unconformity_feature.py

@@ -19,7 +19,7 @@ class UnconformityFeature(GeologicalFeature):
         # just don't link the regions
         GeologicalFeature.__init__(
             self,
-            name=f"{feature.name}_unconformity",
+            name=f"__{feature.name}_unconformity",
             faults=feature.faults,
             regions=[],  # feature.regions.copy(),  # don't want to share regionsbetween unconformity and # feature.regions,
             builder=feature.builder,
@@ -75,9 +75,9 @@ class UnconformityFeature(GeologicalFeature):
             true if above the unconformity, false if below
         """
         if self.sign:
-            return self.evaluate_value(pos) < self.value
+            return self.evaluate_value(pos) <= self.value
         if not self.sign:
-            return self.evaluate_value(pos) > self.value
+            return self.evaluate_value(pos) >= self.value
 
     def __call__(self, pos) -> np.ndarray:
         return self.evaluate(pos)

LoopStructural/modelling/features/fault/_fault_segment.py

@@ -114,6 +114,26 @@ class FaultSegment(StructuralFrame):
     def displacementfeature(self):
         return FaultDisplacementFeature(self, self.faultfunction, name=self.name, model=self.model)
 
+    def fault_ellipsoid(self, **kwargs):
+        try:
+            import pyvista as pv
+
+            fault_ellipsoid = pv.PolyData(
+                self.model.rescale(self.fault_centre[None, :], inplace=False)
+            )
+            fault_ellipsoid["norm"] = self.builder.fault_normal_vector[None, :]
+
+            geom = pv.ParametricEllipsoid(
+                self.fault_minor_axis,
+                self.fault_major_axis,
+                self.fault_intermediate_axis,
+            )
+            ellipsoid = fault_ellipsoid.glyph(geom=geom, **kwargs)
+            return ellipsoid
+        except ImportError:
+            logger.error("pyvista not installed")
+            return None
+
     def set_fault_offset(self, offset: float):
         self.fault_offset = offset
 
@@ -177,7 +197,7 @@ class FaultSegment(StructuralFrame):
         return np.abs(v) > threshold
         # return np.all(np.logical_and(v > -1,v<1),axis=1)
 
-    def evaluate_value(self, locations):
+    def evaluate_value(self, locations, ignore_regions=False):
         """
         Return the value of the fault surface scalar field
 
@@ -202,7 +222,10 @@ class FaultSegment(StructuralFrame):
         #     except:
         #         logger.error("nan slicing")
         # v[mask] = self.__getitem__(0).evaluate_value(locations[mask, :])
-        return self.__getitem__(0).evaluate_value(locations)
+        return self.__getitem__(0).evaluate_value(locations, ignore_regions=ignore_regions)
+
+    def ellipsoid(self):
+        pass
 
     def mean(self):
         return self.__getitem__(0).mean()

LoopStructural/modelling/input/process_data.py

@@ -165,8 +165,19 @@ class ProcessInputData:
             unit_id = 1
             val = self._stratigraphic_value()
             for name, sg in self._stratigraphic_order:
+                # set the oldest unit to be the basement.
+                # this has no observed basal contact and then
+                # top of the unit is the 0 isovalue.
+                # this is the minimum of the next unit
                 stratigraphic_column[name] = {}
-                for i, g in enumerate(reversed(sg)):
+                stratigraphic_column[name][sg[-1]] = {
+                    "max": 0,
+                    "min": -np.inf,
+                    "id": unit_id,
+                    "colour": self.colours[sg[-1]],
+                }
+                # iterate through the remaining units (in reverse)
+                for g in reversed(sg[:-1]):
                     if g in self.thicknesses:
                         stratigraphic_column[name][g] = {
                             "max": val[g] + self.thicknesses[g],
@@ -174,10 +185,6 @@ class ProcessInputData:
                             "id": unit_id,
                             "colour": self.colours[g],
                         }
-                        if i == 0:
-                            stratigraphic_column[name][g]["min"] = 0
-                        if i == len(sg) - 1:
-                            stratigraphic_column[name][g]["max"] = np.inf
 
                     unit_id += 1
             # add faults into the column
@@ -438,7 +445,14 @@ class ProcessInputData:
         stratigraphic_value = {}
         for _name, sg in self.stratigraphic_order:
             value = 0.0  # reset for each supergroup
-            for g in reversed(sg):
+            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
+            for g in reversed(
+                sg[:-1]
+            ):  # don't add the last unit as we never see the base of this unit.
+                # It should be "basement"
                 if g not in self.thicknesses:
                     logger.warning(f"No thicknesses for {g}")
                     stratigraphic_value[g] = np.nan
@@ -469,30 +483,13 @@ class ProcessInputData:
 
     @property
     def contacts(self):
-        return self._contacts
-
-    @contacts.setter
-    def contacts(self, contacts):
-        """Function to convert input contact to loopstructural input
-
-        either uses the thickness values or assigns unique ids given
-        the units named in stratigraphic order
-
-        Returns
-        -------
-        DataFrame
-            data frame with x,y,y,val/interface,feature_name
-        """
-        if contacts is None:
-            return
-        contacts = contacts.copy()
-        self._update_feature_names(contacts)
+        contacts = self._contacts.copy()
         if self._use_thickness:
             contacts["val"] = np.nan
             for k, v in self._stratigraphic_value().items():
                 contacts.loc[contacts["name"] == k, "val"] = v
 
-            self._contacts = contacts.loc[
+            contacts = contacts.loc[
                 ~np.isnan(contacts["val"]), ["X", "Y", "Z", "feature_name", "val"]
             ]
         if not self._use_thickness:
@@ -500,10 +497,28 @@ class ProcessInputData:
             interface_val = 0
             for k in self._stratigraphic_value().keys():
                 contacts.loc[contacts["name"] == k, "interface"] = interface_val
-            self._contacts = contacts.loc[
+            contacts = contacts.loc[
                 ~np.isnan(contacts["interface"]),
                 ["X", "Y", "Z", "feature_name", "interface"],
             ]
+        return contacts
+
+    @contacts.setter
+    def contacts(self, contacts):
+        """Function to convert input contact to loopstructural input
+
+        either uses the thickness values or assigns unique ids given
+        the units named in stratigraphic order
+
+        Returns
+        -------
+        DataFrame
+            data frame with x,y,y,val/interface,feature_name
+        """
+        if contacts is None:
+            return
+        self._contacts = contacts.copy()
+        self._update_feature_names(self._contacts)
 
     @property
     def contact_orientations(self):

LoopStructural/modelling/input/project_file.py

@@ -34,24 +34,48 @@ class LoopProjectfileProcessor(ProcessInputData):
                 projectfile["stratigraphicLog"].ThicknessMedian,
             )
         )
-        fault_properties = self.projectfile.faultLog
-        fault_properties.rename(
-            columns={
-                "avgDisplacement": "displacement",
-                "influenceDistance": "minor_axis",
-                "verticalRadius": "intermediate_axis",
-                "horizontalRadius": "major_axis",
-                "name": "fault_name",
-            },
-            inplace=True,
-        )
-        fault_locations = fault_properties.reset_index()[["fault_name", "eventId"]].merge(
-            fault_locations, on="eventId"
-        )
-        fault_orientations = fault_properties.reset_index()[["fault_name", "eventId"]].merge(
-            fault_orientations, on="eventId"
-        )
-        fault_properties.set_index("fault_name", inplace=True)
+        fault_properties = None
+        fault_edges = None
+        fault_edge_properties = None
+        if self.projectfile.faultLog.shape[0] > 0:
+
+            fault_properties = self.projectfile.faultLog
+            fault_properties.rename(
+                columns={
+                    "avgDisplacement": "displacement",
+                    "influenceDistance": "minor_axis",
+                    "verticalRadius": "intermediate_axis",
+                    "horizontalRadius": "major_axis",
+                    "name": "fault_name",
+                },
+                inplace=True,
+            )
+            fault_locations = fault_properties.reset_index()[["fault_name", "eventId"]].merge(
+                fault_locations, on="eventId"
+            )
+            fault_orientations = fault_properties.reset_index()[["fault_name", "eventId"]].merge(
+                fault_orientations, on="eventId"
+            )
+            fault_properties.set_index("fault_name", inplace=True)
+            for i in fault_relationships.index:
+                fault_relationships.loc[i, "Fault1"] = faultLog.loc[
+                    fault_relationships.loc[i, "eventId1"], "name"
+                ]
+                fault_relationships.loc[i, "Fault2"] = faultLog.loc[
+                    fault_relationships.loc[i, "eventId2"], "name"
+                ]
+            fault_edges = []
+            fault_edge_properties = []
+            for i in fault_relationships.index:
+                fault_edges.append(
+                    (fault_relationships.loc[i, "Fault1"], fault_relationships.loc[i, "Fault2"])
+                )
+                fault_edge_properties.append(
+                    {
+                        "type": fault_relationships.loc[i, "type"],
+                        "angle": fault_relationships.loc[i, "angle"],
+                    }
+                )
         colours = dict(
             zip(
                 self.projectfile.stratigraphicLog.name,
@@ -65,25 +89,6 @@ class LoopProjectfileProcessor(ProcessInputData):
             )
         )
 
-        for i in fault_relationships.index:
-            fault_relationships.loc[i, "Fault1"] = faultLog.loc[
-                fault_relationships.loc[i, "eventId1"], "name"
-            ]
-            fault_relationships.loc[i, "Fault2"] = faultLog.loc[
-                fault_relationships.loc[i, "eventId2"], "name"
-            ]
-        fault_edges = []
-        fault_edge_properties = []
-        for i in fault_relationships.index:
-            fault_edges.append(
-                (fault_relationships.loc[i, "Fault1"], fault_relationships.loc[i, "Fault2"])
-            )
-            fault_edge_properties.append(
-                {
-                    "type": fault_relationships.loc[i, "type"],
-                    "angle": fault_relationships.loc[i, "angle"],
-                }
-            )
         super().__init__(
             contacts=contacts,
             contact_orientations=orientations,
@@ -91,8 +96,8 @@ class LoopProjectfileProcessor(ProcessInputData):
                 ("sg", list(self.projectfile.stratigraphicLog.name))
             ],  # needs to be updated,
             thicknesses=thicknesses,
-            fault_orientations=fault_orientations,
-            fault_locations=fault_locations,
+            fault_orientations=fault_orientations if fault_orientations.shape[0] > 0 else None,
+            fault_locations=fault_locations if fault_locations.shape[0] > 0 else None,
             fault_properties=fault_properties,
             fault_edges=fault_edges,  # list(fault_graph.edges),
             colours=colours,

LoopStructural/utils/__init__.py

@@ -36,3 +36,4 @@ import numpy as np
 rng = np.random.default_rng()
 
 from ._surface import LoopIsosurfacer, surface_list
+from .colours import random_colour

LoopStructural/utils/_surface.py

@@ -1,6 +1,7 @@
 from __future__ import annotations
 
 from typing import Optional, Union, Callable, List
+from collections.abc import Iterable
 import numpy as np
 import numpy.typing as npt
 from LoopStructural.utils.logging import getLogger
@@ -86,11 +87,14 @@ class LoopIsosurfacer:
         all_values = self.callable(self.bounding_box.regular_grid(local=False))
         ## set value to mean value if its not specified
         if values is None:
-            values = [(np.nanmax(all_values) - np.nanmin(all_values)) / 2]
-        if isinstance(values, list):
+            values = [((np.nanmax(all_values) - np.nanmin(all_values)) / 2) + np.nanmin(all_values)]
+        if isinstance(values, Iterable):
             isovalues = values
         elif isinstance(values, float):
             isovalues = [values]
+        if isinstance(values, int) and values == 0:
+            values = 0.0  # assume 0 isosurface is meant to be a float
+
         elif isinstance(values, int) and values < 1:
             raise ValueError(
                 "Number of isosurfaces must be greater than 1. Either use a positive integer or provide a list or float for a specific isovalue."

LoopStructural/utils/colours.py

@@ -1,7 +1,7 @@
 from LoopStructural.utils import rng
 
 
-def random_colour(n: int, cmap='tab20'):
+def random_colour(n: int = 1, cmap='tab20'):
     """
     Generate a list of random colours
 

LoopStructural/version.py

@@ -1 +1 @@
-__version__ = "1.6.3"
+__version__ = "1.6.5"

{LoopStructural-1.6.3.dist-info → LoopStructural-1.6.5.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: LoopStructural
-Version: 1.6.3
+Version: 1.6.5
 Summary: 3D geological modelling
 Author-email: Lachlan Grose <lachlan.grose@monash.edu>
 License: MIT
@@ -23,27 +23,27 @@ Classifier: Programming Language :: Python :: 3.12
 Requires-Python: >=3.8
 Description-Content-Type: text/markdown
 License-File: LICENSE
-Requires-Dist: numpy >=1.18
+Requires-Dist: numpy>=1.18
 Requires-Dist: pandas
 Requires-Dist: scipy
 Requires-Dist: scikit-image
 Requires-Dist: scikit-learn
 Requires-Dist: tqdm
 Provides-Extra: all
-Requires-Dist: loopstructural[export,inequalities,visualisation] ; extra == 'all'
+Requires-Dist: loopstructural[export,inequalities,visualisation]; extra == "all"
+Provides-Extra: visualisation
+Requires-Dist: matplotlib; extra == "visualisation"
+Requires-Dist: pyvista; extra == "visualisation"
+Requires-Dist: loopstructuralviusualisation>0.1.4; extra == "visualisation"
 Provides-Extra: export
-Requires-Dist: geoh5py ; extra == 'export'
-Requires-Dist: pyevtk ; extra == 'export'
-Requires-Dist: dill ; extra == 'export'
-Provides-Extra: inequalities
-Requires-Dist: loopsolver ; extra == 'inequalities'
+Requires-Dist: geoh5py; extra == "export"
+Requires-Dist: pyevtk; extra == "export"
+Requires-Dist: dill; extra == "export"
 Provides-Extra: jupyter
-Requires-Dist: pyvista[all] ; extra == 'jupyter'
-Requires-Dist: tqdm ; extra == 'jupyter'
-Provides-Extra: visualisation
-Requires-Dist: matplotlib ; extra == 'visualisation'
-Requires-Dist: pyvista ; extra == 'visualisation'
-Requires-Dist: loopstructuralviusualisation >0.1.4 ; extra == 'visualisation'
+Requires-Dist: pyvista[all]; extra == "jupyter"
+Requires-Dist: tqdm; extra == "jupyter"
+Provides-Extra: inequalities
+Requires-Dist: loopsolver; extra == "inequalities"
 
 # LoopStructural: Loop3D Implicit Geological Modelling
 

{LoopStructural-1.6.3.dist-info → LoopStructural-1.6.5.dist-info}/RECORD

@@ -1,5 +1,5 @@
 LoopStructural/__init__.py,sha256=nlJ0csTZ_xlTZrWvTKiTp3ER_XRQN6iqQ5sR94RD5QA,1331
-LoopStructural/version.py,sha256=X-erfDySXG8dR_oX8H_L2S8f9Y_VDgdltPb-jH2ruKg,22
+LoopStructural/version.py,sha256=LEtHSENWZLwKaCJasxjOIQSt5Y_OahmKZYacXCmyQtI,22
 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
@@ -30,7 +30,7 @@ LoopStructural/datasets/data/geological_map_data/stratigraphic_orientations.csv,
 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=nTtvNj0UiQWX-61tBH6iWS8YDzHBIqTKWHUxStkKBfU,14732
-LoopStructural/datatypes/_point.py,sha256=rYolPcd134DCj5kvGkTG9_T2S68XOPaBuLfq6-Ud7MU,6140
+LoopStructural/datatypes/_point.py,sha256=ayUVA3EuqAc3NFctiKycGehG2HPNj6ZAMbLANjnq1n4,6291
 LoopStructural/datatypes/_structured_grid.py,sha256=AwSzkoAySU8ZHGbbGTyUoU3fCXuBVQAOD_6th6v0J1o,3665
 LoopStructural/datatypes/_surface.py,sha256=3RLA6-zkWE1L3IvXdrq92Dd5asnI9Lf6ebcNmrMnn98,6445
 LoopStructural/export/exporters.py,sha256=BniZu-PqQvHqCU6GIuJQ5FPzI9Dx_T6rI8EW1pykois,17209
@@ -58,7 +58,7 @@ LoopStructural/interpolators/supports/_2d_structured_grid.py,sha256=TbAV1iw2gmZ4
 LoopStructural/interpolators/supports/_2d_structured_tetra.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 LoopStructural/interpolators/supports/_3d_base_structured.py,sha256=keZnVSK0_wmXCLkWicyIlIV7Wsy2p67OC8B2AOUIMDM,15531
 LoopStructural/interpolators/supports/_3d_p2_tetra.py,sha256=CqGVJRUMxbPQZDhhopNt_s9gVhMqh4YbjQyDZonoyxc,11574
-LoopStructural/interpolators/supports/_3d_structured_grid.py,sha256=iHrZIXQk80vAxkkTzmBRLz3yA6MJEZFKm1RaRrkNlog,17446
+LoopStructural/interpolators/supports/_3d_structured_grid.py,sha256=0H5wh5pBcxqw3crFeXdpnaKMfbSmJNOt-v9-IQa8tSM,17286
 LoopStructural/interpolators/supports/_3d_structured_tetra.py,sha256=5zUNtvEXDvbCHZCu6Fz9WjGbnrMaq-sYJqNUufyLcq8,26505
 LoopStructural/interpolators/supports/_3d_unstructured_tetra.py,sha256=vc9ZffCYk7m1Ae5-H_m1P5OCHgBa9jztrklugMPPiHc,23500
 LoopStructural/interpolators/supports/__init__.py,sha256=V0JjixoBIUZVAo5MmqARR67xDOoQwnb4G3SXeOMRSyQ,1603
@@ -66,18 +66,18 @@ LoopStructural/interpolators/supports/_aabb.py,sha256=Z-kH_u6c6izak0aHG3Uo14PEKQ
 LoopStructural/interpolators/supports/_base_support.py,sha256=5Cwq16EhA2rx4-DFUmoM5z6DRnDfRLkachBpGyzShdM,2420
 LoopStructural/interpolators/supports/_face_table.py,sha256=Hyj4Io63NkPRN8ab9uDHyec-2Kb8BLY_xBF6STNlvBw,3095
 LoopStructural/interpolators/supports/_support_factory.py,sha256=c86NnM4azWhS2ajPApePap0sFI82mZC8siVAU1fCOn4,1175
-LoopStructural/modelling/__init__.py,sha256=IaIRrSvTezC1_RHZpGiVuxZ9tv8XZ4yRLHC2AgXi0kE,787
+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=8TFQho3QolAMwb21uWsr_Ys5rSWHOyrXSVBgQPOnoS4,66092
+LoopStructural/modelling/core/geological_model.py,sha256=GS4VaQBdaWDn9W3hngjeUSgI83A5oDR_IjlcRQTwT_U,66120
 LoopStructural/modelling/features/__init__.py,sha256=qTl-jTau-nnSZ4G4Np6JkBcXB01ViXhouUo_-ImoDvk,877
 LoopStructural/modelling/features/_analytical_feature.py,sha256=7-wFD_h6XDE89jn_trMkLeL98nHAj_5mj802mWTIp7E,2825
-LoopStructural/modelling/features/_base_geological_feature.py,sha256=SP-lWDJDdXZY0UAht33sJkOpW-lGN9NZqP2RuZDgyxg,11102
+LoopStructural/modelling/features/_base_geological_feature.py,sha256=NUF0PbkS5hxKfTCAjlC9nJ0airmLUnkcSYniYG1rpjg,11620
 LoopStructural/modelling/features/_cross_product_geological_feature.py,sha256=IOcJBgQeGBOgwymcyyEdsQHHbVs23sYPpl8xavBaPZE,2781
 LoopStructural/modelling/features/_geological_feature.py,sha256=Migi29Qr_hjIhNYuhclSK4I_MyETwa27lHeClhDKyu0,11221
 LoopStructural/modelling/features/_lambda_geological_feature.py,sha256=pEfLN9HcYGluxRB7LJpcKdYqXP_Mkyo1pLZ3kZZ4mvA,2599
 LoopStructural/modelling/features/_region.py,sha256=TB4qnoTDQM2VgRjgyODN839fKe3kuRYLllJj0xnDKXo,478
-LoopStructural/modelling/features/_structural_frame.py,sha256=I6W-YkQiRADX7xcKMFzDefYGDlXWVAp5nAVRFVcqNow,5062
-LoopStructural/modelling/features/_unconformity_feature.py,sha256=-0shSjYdT8IuNDNDYnqJ2niFjFFiBt-Izc-oMBYa0MQ,2405
+LoopStructural/modelling/features/_structural_frame.py,sha256=Jl89dL_GsGWvtXibAr4zitVChKw9sog0-Rs9QMHw8q4,4802
+LoopStructural/modelling/features/_unconformity_feature.py,sha256=6SAgPkjipRR0Wo4aXblKyT2YuIphZc6ys1rbLXt4Umw,2409
 LoopStructural/modelling/features/builders/__init__.py,sha256=Gqld1C-PcaXfJ8vpkWMDCmehmd3hZNYQk1knPtl59Bk,266
 LoopStructural/modelling/features/builders/_base_builder.py,sha256=ImlKZs4sFSrtAAp1NEpQD35AKF_mLEHr7jj2DzSyAd4,3561
 LoopStructural/modelling/features/builders/_fault_builder.py,sha256=EyDLyovrnxqSk11zWFnViok7lXMt_E-pG2QqgRgvOLk,25653
@@ -87,7 +87,7 @@ LoopStructural/modelling/features/builders/_structural_frame_builder.py,sha256=1
 LoopStructural/modelling/features/fault/__init__.py,sha256=4u0KfYzmoO-ddFGo9qd9ov0gBoLqBiPAUsaw5zhEOAQ,189
 LoopStructural/modelling/features/fault/_fault_function.py,sha256=5IzVyvv1tS5Z6geMB8FdTZpGrBYbyu38cOH03NACm0o,12797
 LoopStructural/modelling/features/fault/_fault_function_feature.py,sha256=1F_xGGRtLZ7mBnXutEGkwWoJImSEPUuxstf0Nu-7k5w,2452
-LoopStructural/modelling/features/fault/_fault_segment.py,sha256=cAnf8cKYJd-3Mb8Q6F1ZVAuez2jkHZhRN4yAz_aMNfI,16947
+LoopStructural/modelling/features/fault/_fault_segment.py,sha256=upZtkxdOPkR1bxum8MO8xUbJo_e3HRvl3JTJD1Y-L50,17723
 LoopStructural/modelling/features/fold/__init__.py,sha256=mOq_K5IHjE3FoK7dhYpljUMGj0gvEX2Fyqtzjd4EYfQ,176
 LoopStructural/modelling/features/fold/_fold.py,sha256=bPnnLUSiF4uoMRg8aHoOSTPRgaM0JyLoRQPu5_A-J3w,5448
 LoopStructural/modelling/features/fold/_fold_rotation_angle_feature.py,sha256=CXLbFRQ3CrTMAcHmfdbKcmSvvLs9_6TLe0Wqi1pK2tg,892
@@ -101,8 +101,8 @@ 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=Y3wwUigqc9tqaX4EOzAOAigepzbtO6VWFhl_YmF00nA,25441
-LoopStructural/modelling/input/project_file.py,sha256=BYV9KBkAFHubP3QOWbnPjOPpaoyhh4BuZU1I194wJUs,4221
+LoopStructural/modelling/input/process_data.py,sha256=dmA6F-l8nvdH3Ckp5OaAjl3NJoWgm3PuELdXtBRGrz8,26149
+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
 LoopStructural/modelling/intrusions/geometric_scaling_functions.py,sha256=PK3qf0TiK-WYIBGG7fYhTD1hwlUN0s75BK8d53SLYuQ,3209
@@ -110,10 +110,10 @@ LoopStructural/modelling/intrusions/intrusion_builder.py,sha256=1cJjPyRUf3ZDkpwg
 LoopStructural/modelling/intrusions/intrusion_feature.py,sha256=5DFQ1uccX7Rp3PMc_pidw-HixDtKyiRhfLkIUye_Gcc,15455
 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=K6OwdnC8RE-sDhG4GIbVyoJNUD0WnM1VGG6os1l_Woc,870
-LoopStructural/utils/_surface.py,sha256=8pUn11Pe3ZlYoHkjy06v3QWJ7UJORQ7yWsgwFQptkxs,5530
+LoopStructural/utils/__init__.py,sha256=-tyTlOr939RMgya5LW4HY-lELviTgClPAZt8L0AAE4s,905
+LoopStructural/utils/_surface.py,sha256=fPJIHVk8gs6Ws3kEW1L4gTZS8znZTNrEm8WVdtm0ZII,5721
 LoopStructural/utils/_transformation.py,sha256=7iDPMIBoZ73ZGa_slbC9XzdMQ-ONOxvF7UkY4dJaw54,2443
-LoopStructural/utils/colours.py,sha256=NOYTN-fh-XTfQdyqRVzO5bdt0kJQrnxH25ZxEnqp8ZQ,548
+LoopStructural/utils/colours.py,sha256=39MKnLiama0Ubhlxjhau8KtUGqBrluW0vw25qhuyJS0,552
 LoopStructural/utils/config.py,sha256=ITGOtZTo2_QBwXkG_0AFANfE90J9siCXLzxypVmg9QA,414
 LoopStructural/utils/dtm_creator.py,sha256=-yqGG0wyEJfTCCDghz058wull1q3zGFASjeu8oDgYnk,535
 LoopStructural/utils/exceptions.py,sha256=SJboJ7ncMqVX-ib7MMizClwMrFZRHQhjZr2eCnVwnQE,500
@@ -127,8 +127,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.3.dist-info/LICENSE,sha256=ZqGeNFOgmYevj7Ld7Q-kR4lAxWXuBRUdUmPC6XM_py8,1071
-LoopStructural-1.6.3.dist-info/METADATA,sha256=AQ55C8hbLHvryNgxtgqAnp--xB9ZjKWmrBhg3nGGfCQ,5793
-LoopStructural-1.6.3.dist-info/WHEEL,sha256=OVMc5UfuAQiSplgO0_WdW7vXVGAt9Hdd6qtN4HotdyA,91
-LoopStructural-1.6.3.dist-info/top_level.txt,sha256=QtQErKzYHfg6ddxTQ1NyaTxXBVM6qAqrM_vxEPyXZLg,15
-LoopStructural-1.6.3.dist-info/RECORD,,
+LoopStructural-1.6.5.dist-info/LICENSE,sha256=ZqGeNFOgmYevj7Ld7Q-kR4lAxWXuBRUdUmPC6XM_py8,1071
+LoopStructural-1.6.5.dist-info/METADATA,sha256=69VHE2kjFbIx1loVcB7tV9ExJw14OgGDMka_633-Mzs,5781
+LoopStructural-1.6.5.dist-info/WHEEL,sha256=PZUExdf71Ui_so67QXpySuHtCi3-J3wvF4ORK6k_S8U,91
+LoopStructural-1.6.5.dist-info/top_level.txt,sha256=QtQErKzYHfg6ddxTQ1NyaTxXBVM6qAqrM_vxEPyXZLg,15
+LoopStructural-1.6.5.dist-info/RECORD,,

{LoopStructural-1.6.3.dist-info → LoopStructural-1.6.5.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (75.2.0)
+Generator: setuptools (75.6.0)
 Root-Is-Purelib: true
 Tag: py3-none-any