fix: removing normalisation from

orthogonal constraint
Loop3D · Sep 21, 2021 · 3216d07 · 3216d07
1 parent f43a583
commit 3216d07
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/LoopStructural/interpolators/finite_difference_interpolator.py b/LoopStructural/interpolators/finite_difference_interpolator.py
@@ -353,10 +353,11 @@ def add_gradient_orthogonal_constraint(self, points, vector, w=1.0,
 
             #normalise element vector to unit vector for dot product
             T = self.support.calcul_T(points[inside, :3])
-            norm = np.linalg.norm(T,axis=1)
-            T/=norm[:,None,:]
+            # norm = np.linalg.norm(T,axis=1)
+            # T/=norm[:,None,:]
             # normalise vector to unit vector for dot product
-            vector[inside,:3] /= np.linalg.norm(vector[inside,:3],axis=1)[:,None]
+            # nro
+            # # vector[inside,:3] /= np.linalg.norm(vector[inside,:3],axis=1)[:,None]
             # dot product of vector and element gradient 
             A = np.einsum('ij,ijk->ik', vector[inside, :3], T)