fix: updating product to prod

Author: lachlangrose

LoopStructural/interpolators/_builders.py

@@ -7,6 +7,7 @@
     FiniteDifferenceInterpolator,
     GeologicalInterpolator,
     DiscreteInterpolator,
+    DiscreteFoldInterpolator,
     StructuredGrid,
     TetMesh,
 )
@@ -31,7 +32,7 @@ def get_interpolator(
     # buffer = bb[1, :]
     origin = bounding_box.with_buffer(buffer).origin
     maximum = bounding_box.with_buffer(buffer).maximum
-    box_vol = np.product(maximum - origin)
+    box_vol = np.prod(maximum - origin)
     if interpolatortype == "PLI":
         if support is None:
             if element_volume is None:
@@ -98,6 +99,23 @@ def get_interpolator(
             "for modelling using FDI"
         )
         return FiniteDifferenceInterpolator(grid)
+    if interpolatortype == "DFI":
+        if element_volume is None:
+            nelements /= 5
+            element_volume = box_vol / nelements
+        # calculate the step vector of a regular cube
+        step_vector = np.zeros(3)
+        step_vector[:] = element_volume ** (1.0 / 3.0)
+        # number of steps is the length of the box / step vector
+        nsteps = np.ceil((maximum - origin) / step_vector).astype(int)
+        # create a structured grid using the origin and number of steps
+
+        mesh = TetMesh(origin=origin, nsteps=nsteps, step_vector=step_vector)
+        logger.info(
+            f"Creating regular tetrahedron mesh with {mesh.ntetra} elements \n"
+            "for modelling using DFI"
+        )
+        return DiscreteFoldInterpolator(mesh, None)
     raise LoopException("No interpolator")
     # fi interpolatortype == "DFI" and dfi is True:
     #     if element_volume is None:

LoopStructural/interpolators/supports/_3d_base_structured.py

@@ -134,11 +134,11 @@ def maximum(self, maximum):
 
     @property
     def n_nodes(self):
-        return np.product(self.nsteps)
+        return np.prod(self.nsteps)
 
     @property
     def n_elements(self):
-        return np.product(self.nsteps_cells)
+        return np.prod(self.nsteps_cells)
 
     def __str__(self):
         return (

LoopStructural/interpolators/supports/_3d_structured_tetra.py

@@ -32,15 +32,15 @@ def __init__(
 
     @property
     def ntetra(self) -> int:
-        return np.product(self.nsteps_cells) * 5
+        return np.prod(self.nsteps_cells) * 5
 
     @property
     def n_elements(self) -> int:
         return self.ntetra
 
     @property
     def n_cells(self) -> int:
-        return np.product(self.nsteps_cells)
+        return np.prod(self.nsteps_cells)
 
     @property
     def barycentre(self) -> np.ndarray:

LoopStructural/interpolators/supports/_3d_unstructured_tetra.py

@@ -60,7 +60,7 @@ def __init__(
         self.minimum -= length * 0.1
         self.maximum += length * 0.1
         if aabb_nsteps == None:
-            box_vol = np.product(self.maximum - self.minimum)
+            box_vol = np.prod(self.maximum - self.minimum)
             element_volume = box_vol / (len(self.elements) / 20)
             # calculate the step vector of a regular cube
             step_vector = np.zeros(3)