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

Author: pancetta

.travis.yml

@@ -16,9 +16,9 @@ before_install:
   - conda update --yes conda
 
 install:
-  - conda install --yes python=${TRAVIS_PYTHON_VERSION} pip numpy future scipy matplotlib nose pep8 coverage mpi4py sympy dill numba
-  - conda install --yes sphinx
-  - conda install --yes -c conda-forge petsc=3.9.0 petsc4py=3.9.0
+  - conda install --yes python=${TRAVIS_PYTHON_VERSION} pip numpy future scipy matplotlib nose pep8 coverage mpi4py sympy dill numba sphinx
+#  - conda install --yes sphinx
+  - conda install --yes -c conda-forge petsc4py
   - conda list
   - pip install coloredlogs
   - pip install travis-sphinx

pySDC/implementations/problem_classes/HeatEquation_2D_PETSc_forced.py

@@ -1,6 +1,5 @@
 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
@@ -32,6 +31,8 @@ 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:
@@ -61,6 +62,9 @@ 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()
@@ -91,16 +95,14 @@ 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):
-                row.index = (i, j)
-                row.field = 0
+                self.row.index = (i, j)
+                self.row.field = 0
                 if i == 0 or j == 0 or i == mx - 1 or j == my - 1:
-                    A.setValueStencil(row, row, 1.0)
+                    A.setValueStencil(self.row, self.row, 1.0)
                 else:
                     diag = self.params.nu * (-2.0 / self.dx ** 2 - 2.0 / self.dy ** 2)
                     for index, value in [
@@ -110,9 +112,9 @@ def __get_A(self):
                         ((i + 1, j), self.params.nu / self.dx ** 2),
                         ((i, j + 1), self.params.nu / self.dy ** 2),
                     ]:
-                        col.index = index
-                        col.field = 0
-                        A.setValueStencil(row, col, value)
+                        self.col.index = index
+                        self.col.field = 0
+                        A.setValueStencil(self.row, self.col, value)
         A.assemble()
 
         return A
@@ -134,13 +136,12 @@ 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):
-                row.index = (i, j)
-                row.field = 0
-                Id.setValueStencil(row, row, 1.0)
+                self.row.index = (i, j)
+                self.row.field = 0
+                Id.setValueStencil(self.row, self.row, 1.0)
 
         Id.assemble()
 

pySDC/implementations/transfer_classes/TransferPETScDMDA.py

@@ -1,5 +1,4 @@
 from __future__ import division
-from petsc4py import PETSc
 
 from pySDC.core.SpaceTransfer import space_transfer
 from pySDC.core.Errors import TransferError
@@ -26,7 +25,7 @@ def __init__(self, fine_prob, coarse_prob, params):
         super(mesh_to_mesh_petsc_dmda, self).__init__(fine_prob, coarse_prob, params)
 
         # set interpolation type (no effect as far as I can tell)
-        self.coarse_prob.init.setInterpolationType(PETSc.DMDA.InterpolationType.Q1)
+        # self.coarse_prob.init.setInterpolationType(PETSc.DMDA.InterpolationType.Q1)
         # define interpolation (only accurate for constant functions)
         self.interp, _ = self.coarse_prob.init.createInterpolation(self.fine_prob.init)
         # define restriction as injection (tranpose of interpolation does not work)

pySDC/projects/AsympConv/conv_test_to0.py

@@ -56,7 +56,8 @@ 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((1 + x) ** (1.0 - 1.0 / (Nnodes * Nsteps)) * x ** (1.0 / (Nnodes * Nsteps)))
+            predict_list.append(x ** (1.0 / (Nnodes * Nsteps)))
 
             if len(predict_list) > 1:
                 print(x, predict_list[-1], Prho_list[-1], Prho_list[-2] / Prho_list[-1],