Add forgotten flat_eval

src/gvec_to_python/geometry/domain.py

@@ -1110,21 +1110,21 @@ class GVEC_domain:
         '''
         return self.Pmap(*self.Lmap(s, u, v, flat_eval=flat_eval), flat_eval=flat_eval)
     
-    def dPL(self, s, u, v):
+    def dPL(self, s, u, v, flat_eval=False):
         '''Jacobian matrix of the composition P ° L.'''
 
-        tmp1 = self.dL(s, u, v)  # Jacobian of the Lmap
-        tmp2 = self.dP(*self.Lmap(s, u, v))  # Jacobian of the Pmap
+        tmp1 = self.dL(s, u, v, flat_eval=flat_eval)  # Jacobian of the Lmap
+        tmp2 = self.dP(*self.Lmap(s, u, v, flat_eval=flat_eval), flat_eval=flat_eval)  # Jacobian of the Pmap
 
         return tmp2 @ tmp1
 
-    def dPL_inv(self, s, u, v):
+    def dPL_inv(self, s, u, v, flat_eval=False):
         '''Inverse Jacobian matrix of the composition P ° L.'''
-        return np.linalg.inv(self.dPL(s, u, v))
+        return np.linalg.inv(self.dPL(s, u, v, flat_eval=flat_eval))
 
-    def det_dPL(self, s, u, v):
+    def det_dPL(self, s, u, v, flat_eval=False):
         '''Jacobian determinant of the composition P ° L.'''
-        return np.linalg.det(self.dPL(s, u, v))
+        return np.linalg.det(self.dPL(s, u, v, flat_eval=flat_eval))
 
     @staticmethod
     def prepare_args(s, u, v, sparse=False, flat_eval=False):