fix: number of constraints warning message

Loop3D · Jan 19, 2022 · 9abcc04 · 9abcc04
1 parent f77f915
commit 9abcc04
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 5 deletions.
diff --git a/LoopStructural/interpolators/finite_difference_interpolator.py b/LoopStructural/interpolators/finite_difference_interpolator.py
@@ -177,7 +177,7 @@ def add_vaue_constraints(self, w=1.0):
                 w=w * points[inside, 4],
                 name="value",
             )
-            if np.sum(inside)>0:
+            if np.sum(inside)<=0:
                 logger.warning(f"{self.propertyname}: {np.sum(~inside)} value constraints not added: outside of model bounding box")
     def add_interface_constraints(
         self, w=1.0
@@ -295,7 +295,7 @@ def add_gradient_constraints(self, w=1.0):
             self.add_constraints_to_least_squares(
                 A, B, idc[inside, :], w=w * self.vol, name="gradient"
             )
-            if np.sum(inside)>0:
+            if np.sum(inside)<=0:
                 logger.warning(f"{self.propertyname}: {np.sum(~inside)} norm constraints not added: outside of model bounding box")
     def add_norm_constraints(self, w=1.0):
         """
@@ -356,7 +356,7 @@ def add_norm_constraints(self, w=1.0):
                 w=w * self.vol,
                 name="norm",
             )
-            if np.sum(inside)>0:
+            if np.sum(inside)<=0:
                 logger.warning(f"{self.propertyname}: {np.sum(~inside)} norm constraints not added: outside of model bounding box")
 
     def add_gradient_orthogonal_constraints(self, points, vector, w=1.0, B=0):
@@ -408,7 +408,7 @@ def add_gradient_orthogonal_constraints(self, points, vector, w=1.0, B=0):
             self.add_constraints_to_least_squares(
                 A, B, idc[inside, :], w=w * self.vol, name="gradient orthogonal"
             )
-            if np.sum(inside)>0:
+            if np.sum(inside)<=0:
                 logger.warning(f"{self.propertyname}: {np.sum(~inside)} gradient constraints not added: outside of model bounding box")
 
     def add_regularisation(self, operator, w=0.1):

diff --git a/LoopStructural/modelling/fault/fault_builder.py b/LoopStructural/modelling/fault/fault_builder.py
@@ -263,7 +263,6 @@ def set_mesh_geometry(self, buffer, rotation):
         length = np.max(self.maximum - self.origin)
         # for builder in self.builders:
         # all three coordinates share the same support
-        print(self.maximum,self.origin)
         self.builders[0].set_interpolation_geometry(
             self.origin - length * buffer, self.maximum + length * buffer, rotation
         )