fix: default regularisation should be 0.1 for both FDI and PLI.

Author: lachlangrose

LoopStructural/interpolators/_finite_difference_interpolator.py

@@ -454,22 +454,28 @@ def add_gradient_orthogonal_constraints(
                 )
             self.up_to_date = False
 
-    def add_regularisation(self, operator, w=0.1):
+
+
+    def assemble_borders(self, operator, w, name='regularisation'):
         """
+        Adds a constraint to the border of the model to force the value to be equal to the value at the border
 
         Parameters
         ----------
-        operator
-        w
+        operator : Operator
+            operator to use for the regularisation
+        w : double
+            weight of the regularisation
 
         Returns
         -------
 
         """
-        self.assemble_inner(operator, w)
-        # self.assemble_borders()
+        # First get the global indicies of the pairs of neighbours this should be an
+        # Nx27 array for 3d and an Nx9 array for 2d
 
-    # def assemble_borders(self, operator, w):
+        global_indexes = self.support.neighbour_global_indexes()
+        
 
     def assemble_inner(self, operator, w, name='regularisation'):
         """

LoopStructural/interpolators/_p1interpolator.py

@@ -146,7 +146,7 @@ def setup_interpolator(self, **kwargs):
         self.reset()
         for key in kwargs:
             if "regularisation" in kwargs:
-                self.interpolation_weights["cgw"] = 0.1 * kwargs["regularisation"]
+                self.interpolation_weights["cgw"] = kwargs["regularisation"]
             self.up_to_date = False
             self.interpolation_weights[key] = kwargs[key]
         if self.interpolation_weights["cgw"] > 0.0:

LoopStructural/modelling/features/builders/_geological_feature_builder.py

@@ -86,7 +86,7 @@ def __init__(
         self._orthogonal_features = {}
         self._equality_constraints = {}
         # add default parameters
-        self.update_build_arguments({'cpw':1.0,'npw':1.0,'regularisation':1.0,'nelements':self.interpolator.n_elements})
+        self.update_build_arguments({'cpw':1.0,'npw':1.0,'regularisation':.10,'nelements':self.interpolator.n_elements})
     def set_not_up_to_date(self, caller):
         logger.info(
             f"Setting {self.name} to not up to date from an instance of {caller.__class__.__name__}"