hotfix for sp.diags

Author: pancetta

projects/parallelSDC/GeneralizedFisher_1D_FD_implicit_Jac.py

@@ -20,7 +20,7 @@ def eval_jacobian(self, u):
 
         # noinspection PyTypeChecker
         dfdu = self.A[1:-1, 1:-1] + sp.diags(self.params.lambda0 ** 2 - self.params.lambda0 ** 2 *
-                                             (self.params.nu + 1) * u.values ** self.params.nu)
+                                             (self.params.nu + 1) * u.values ** self.params.nu, offsets=0)
 
         return dfdu
 

pySDC/implementations/problem_classes/GeneralizedFisher_1D_FD_implicit.py

@@ -106,8 +106,7 @@ def solve_system(self, rhs, factor, u0, t):
                 break
 
             # assemble dg
-            dg = sp.eye(self.params.nvars) - factor * (self.A[1:-1, 1:-1] + sp.diags(lambda0 ** 2 - lambda0 ** 2 *
-                                                                                     (nu + 1) * u.values ** nu))
+            dg = sp.eye(self.params.nvars) - factor * (self.A[1:-1, 1:-1] + sp.diags(lambda0 ** 2 - lambda0 ** 2 * (nu + 1) * u.values ** nu, offsets=0))
 
             # newton update: u1 = u0 - g/dg
             u.values -= spsolve(dg, g)