Multiple SED in ChromaticTransformation by absolute value of determinant

Author: FedericoBerlfein

galsim/chromatic.py

@@ -1903,9 +1903,9 @@ def sed(self):
             @np.vectorize
             def detjac(w):
                 return np.linalg.det(np.asarray(self._jac(w)).reshape(2,2))
-            sed *= detjac
+            sed *= np.abs(detjac)
         elif self._jac is not None:
-            sed *= np.linalg.det(np.asarray(self._jac).reshape(2,2))
+            sed *= np.abs(np.linalg.det(np.asarray(self._jac).reshape(2,2)))
 
         return sed