From 2f382640a3db431f33978153205530bb286f66c4 Mon Sep 17 00:00:00 2001 From: Stefan Possanner <spossann@ipp.mpg.de> Date: Tue, 19 Mar 2024 14:12:20 +0100 Subject: [PATCH] Add forgotten flat_eval --- src/gvec_to_python/geometry/domain.py | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/src/gvec_to_python/geometry/domain.py b/src/gvec_to_python/geometry/domain.py index 8ea0710..d0a1e9c 100644 --- a/src/gvec_to_python/geometry/domain.py +++ b/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): -- GitLab