fix: gradient constraints in pli weighted by vol

Author: Lachlan Grose

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')