debugging.............

Author: pancetta

.travis.yml

@@ -1,4 +1,4 @@
-sudo: required
+sudo: false
 dist: trusty
 
 language: python

pySDC/implementations/problem_classes/HeatEquation_2D_PETSc_forced.py

@@ -1,5 +1,6 @@
 from __future__ import division
 import numpy as np
+from petsc4py import PETSc
 
 from pySDC.core.Problem import ptype
 from pySDC.core.Errors import ParameterError, ProblemError
@@ -31,8 +32,6 @@ def __init__(self, problem_params, dtype_u, dtype_f):
 
         # these parameters will be used later, so assert their existence
 
-        from petsc4py import PETSc
-
         if 'comm' not in problem_params:
             problem_params['comm'] = PETSc.COMM_WORLD
         if 'sol_tol' not in problem_params:
@@ -62,9 +61,6 @@ def __init__(self, problem_params, dtype_u, dtype_f):
         self.dy = 1.0 / (self.params.nvars[1] - 1)
         (self.xs, self.xe), (self.ys, self.ye) = self.init.getRanges()
 
-        self.row = PETSc.Mat.Stencil()
-        self.col = PETSc.Mat.Stencil()
-
         # compute discretization matrix A and identity
         self.A = self.__get_A()
         self.Id = self.__get_Id()
@@ -95,14 +91,16 @@ def __get_A(self):
 
         # fill matrix
         A.zeroEntries()
+        row = PETSc.Mat.Stencil()
+        col = PETSc.Mat.Stencil()
         mx, my = self.init.getSizes()
         (xs, xe), (ys, ye) = self.init.getRanges()
         for j in range(ys, ye):
             for i in range(xs, xe):
-                self.row.index = (i, j)
-                self.row.field = 0
+                row.index = (i, j)
+                row.field = 0
                 if i == 0 or j == 0 or i == mx - 1 or j == my - 1:
-                    A.setValueStencil(self.row, self.row, 1.0)
+                    A.setValueStencil(row, row, 1.0)
                 else:
                     diag = self.params.nu * (-2.0 / self.dx ** 2 - 2.0 / self.dy ** 2)
                     for index, value in [
@@ -112,9 +110,9 @@ def __get_A(self):
                         ((i + 1, j), self.params.nu / self.dx ** 2),
                         ((i, j + 1), self.params.nu / self.dy ** 2),
                     ]:
-                        self.col.index = index
-                        self.col.field = 0
-                        A.setValueStencil(self.row, self.col, value)
+                        col.index = index
+                        col.field = 0
+                        A.setValueStencil(row, col, value)
         A.assemble()
 
         return A
@@ -136,12 +134,13 @@ def __get_Id(self):
 
         # fill matrix
         Id.zeroEntries()
+        row = PETSc.Mat.Stencil()
         (xs, xe), (ys, ye) = self.init.getRanges()
         for j in range(ys, ye):
             for i in range(xs, xe):
-                self.row.index = (i, j)
-                self.row.field = 0
-                Id.setValueStencil(self.row, self.row, 1.0)
+                row.index = (i, j)
+                row.field = 0
+                Id.setValueStencil(row, row, 1.0)
 
         Id.assemble()
 

pySDC/projects/AsympConv/conv_test_to0.py

@@ -56,8 +56,7 @@ def compute_and_plot_specrad(Nnodes, lam):
                 x * lam * np.kron(np.eye(Nsteps), (Qmat - QDmat)) + np.kron(Emat, Nmat))
 
             Prho_list.append(max(abs(np.linalg.eigvals(mat))))
-            # predict_list.append((1 + x) ** (1.0 - 1.0 / (Nnodes * Nsteps)) * x ** (1.0 / (Nnodes * Nsteps)))
-            predict_list.append(x ** (1.0 / (Nnodes * Nsteps)))
+            predict_list.append((1 + x) ** (1.0 - 1.0 / (Nnodes * Nsteps)) * x ** (1.0 / (Nnodes * Nsteps)))
 
             if len(predict_list) > 1:
                 print(x, predict_list[-1], Prho_list[-1], Prho_list[-2] / Prho_list[-1],