fix: gradient constraints in pli weighted by vol

Loop3D · Nov 8, 2021 · 9b487fd · 9b487fd
1 parent 78b9eab
commit 9b487fd
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/LoopStructural/interpolators/piecewiselinear_interpolator.py b/LoopStructural/interpolators/piecewiselinear_interpolator.py
@@ -208,7 +208,7 @@ def add_gradient_constraints(self, w=1.0):
             w*=points[:,6]
             self.add_constraints_to_least_squares(A[outside, :],
                                                   B[outside], idc[outside, :],
-                                                  w=w*vol[outside],
+                                                  w=w,
                                                   name = 'gradient_strike')
             A = np.einsum('ji,ijk->ik', dip_vector, element_gradients)
             # A *= vol[:, None]
@@ -388,6 +388,6 @@ def add_gradient_orthogonal_constraints(self, points, vector, w=1.0,
             B = np.zeros(idc.shape[0])+B
             outside = ~np.any(idc == -1, axis=1)
             self.add_constraints_to_least_squares(A[outside, :],
-                                                  B[outside], idc[outside, :], w=w*vol[outside],name='gradient orthogonal')
+                                                  B[outside], idc[outside, :], w=w,name='gradient orthogonal')