Merge pull request #19 from pymor/fix_amax

fix DealIIVector amax
The vector binding to pyMOR was missing taking the absolute
values of the contained vector before computing
and returning the max+index

Author: renefritze

src/pymor_dealii/pymor/vectorarray.py

@@ -75,8 +75,9 @@ def dofs(self, dof_indices):
         return np.array([self.impl[i] for i in dof_indices])
 
     def amax(self):
-        max_ind = np.argmax(self.impl)
-        return max_ind, self.impl[max_ind]
+        A = np.abs(self.to_numpy())
+        max_ind = np.argmax(A)
+        return max_ind, A[max_ind]
 
     def __add__(self, other):
         return DealIIVector(self.impl + other.impl)