flake8 compatibility

pancetta · pancetta · commit d047b16bfb68 · 2021-08-10T08:46:12.000+02:00
diff --git a/.flake8 b/.flake8
@@ -0,0 +1,9 @@
+[flake8]
+max-line-length = 120
+ignore = E741, E402, W504, W605
+per-file-ignores =
+    pySDC/tutorial/step_6/C_MPI_parallelization.py:F401,
+    pySDC/projects/Hamiltonian/solar_system.py:F401
+exclude =
+    playgrounds
+    tests
diff --git a/pySDC/helpers/transfer_helper.py b/pySDC/helpers/transfer_helper.py
@@ -266,7 +266,7 @@ def border_padding(grid, l, r, pad_type='mirror'):
 
     assert l < grid.size and r < grid.size
     padded_arr = np.zeros(grid.size + l + r)
-    if pad_type is 'mirror':
+    if pad_type == 'mirror':
         for i in range(l):
             padded_arr[i] = 2 * grid[0] - grid[l - i]
         for j in range(r):
diff --git a/pySDC/helpers/visualization_tools.py b/pySDC/helpers/visualization_tools.py
@@ -65,4 +65,4 @@ def show_residual_across_simulation(stats, fname='residuals.png'):
     ax.set_xticklabels(np.arange(maxiter) + 1, minor=False)
     ax.set_yticklabels(np.arange(maxprocs + 1), minor=False)
 
-    plt.savefig(fname,  transparent=True, bbox_inches='tight')
+    plt.savefig(fname, transparent=True, bbox_inches='tight')
diff --git a/pySDC/implementations/controller_classes/controller_nonMPI.py b/pySDC/implementations/controller_classes/controller_nonMPI.py
@@ -675,15 +675,15 @@ def default(local_MS_running):
             raise ControllerError('Unknown stage, got %s' % local_MS_running[0].status.stage)  # TODO
 
         # if all stages are the same (or DONE), continue, otherwise abort
-        stages = [S.status.stage for S in local_MS_active if S.status.stage is not 'DONE']
+        stages = [S.status.stage for S in local_MS_active if S.status.stage != 'DONE']
         if stages[1:] == stages[:-1]:
             stage = stages[0]
         else:
             raise ControllerError('not all stages are equal')
 
         self.logger.debug(stage)
 
-        MS_running = [S for S in local_MS_active if S.status.stage is not 'DONE']
+        MS_running = [S for S in local_MS_active if S.status.stage != 'DONE']
 
         switcher = {
             'SPREAD': spread,
diff --git a/pySDC/implementations/problem_classes/AllenCahn_Temp_MPIFFT.py b/pySDC/implementations/problem_classes/AllenCahn_Temp_MPIFFT.py
@@ -1,5 +1,4 @@
 import numpy as np
-from mpi4py import MPI
 from mpi4py_fft import PFFT
 
 from pySDC.core.Errors import ParameterError, ProblemError
diff --git a/pySDC/implementations/problem_classes/FullSolarSystem.py b/pySDC/implementations/problem_classes/FullSolarSystem.py
@@ -51,26 +51,26 @@ def u_exact(self, t):
         # initial positions and velocities taken from
         # https://www.aanda.org/articles/aa/full/2002/08/aa1405/aa1405.right.html
         me.pos[:, 0] = [0.0, 0.0, 0.0]
-        me.pos[:, 1] = [-2.503321047836E-01,	+1.873217481656E-01, +1.260230112145E-01]
-        me.pos[:, 2] = [+1.747780055994E-02,	-6.624210296743E-01, -2.991203277122E-01]
-        me.pos[:, 3] = [-9.091916173950E-01,	+3.592925969244E-01, +1.557729610506E-01]
-        me.pos[:, 4] = [+1.203018828754E+00,	+7.270712989688E-01, +3.009561427569E-01]
-        me.pos[:, 5] = [+3.733076999471E+00,	+3.052424824299E+00, +1.217426663570E+00]
+        me.pos[:, 1] = [-2.503321047836E-01, +1.873217481656E-01, +1.260230112145E-01]
+        me.pos[:, 2] = [+1.747780055994E-02, -6.624210296743E-01, -2.991203277122E-01]
+        me.pos[:, 3] = [-9.091916173950E-01, +3.592925969244E-01, +1.557729610506E-01]
+        me.pos[:, 4] = [+1.203018828754E+00, +7.270712989688E-01, +3.009561427569E-01]
+        me.pos[:, 5] = [+3.733076999471E+00, +3.052424824299E+00, +1.217426663570E+00]
         me.pos[:, 6] = [+6.164433062913E+00, +6.366775402981E+00, +2.364531109847E+00]
-        me.pos[:, 7] = [+1.457964661868E+01,	-1.236891078519E+01, -5.623617280033E+00]
-        me.pos[:, 8] = [+1.695491139909E+01,	-2.288713988623E+01, -9.789921035251E+00]
+        me.pos[:, 7] = [+1.457964661868E+01, -1.236891078519E+01, -5.623617280033E+00]
+        me.pos[:, 8] = [+1.695491139909E+01, -2.288713988623E+01, -9.789921035251E+00]
         me.pos[:, 9] = [-9.707098450131E+00, -2.804098175319E+01, -5.823808919246E+00]
 
         me.vel[:, 0] = [0.0, 0.0, 0.0]
         me.vel[:, 1] = [-2.438808424736E-02, -1.850224608274E-02, -7.353811537540E-03]
         me.vel[:, 2] = [+2.008547034175E-02, +8.365454832702E-04, -8.947888514893E-04]
-        me.vel[:, 3] = [-7.085843239142E-03,	-1.455634327653E-02, -6.310912842359E-03]
-        me.vel[:, 4] = [-7.124453943885E-03,	+1.166307407692E-02, +5.542098698449E-03]
-        me.vel[:, 5] = [-5.086540617947E-03,	+5.493643783389E-03, +2.478685100749E-03]
+        me.vel[:, 3] = [-7.085843239142E-03, -1.455634327653E-02, -6.310912842359E-03]
+        me.vel[:, 4] = [-7.124453943885E-03, +1.166307407692E-02, +5.542098698449E-03]
+        me.vel[:, 5] = [-5.086540617947E-03, +5.493643783389E-03, +2.478685100749E-03]
         me.vel[:, 6] = [-4.426823593779E-03, +3.394060157503E-03, +1.592261423092E-03]
         me.vel[:, 7] = [+2.647505630327E-03, +2.487457379099E-03, +1.052000252243E-03]
-        me.vel[:, 8] = [+2.568651772461E-03,	+1.681832388267E-03, +6.245613982833E-04]
-        me.vel[:, 9] = [+3.034112963576E-03,	-1.111317562971E-03, -1.261841468083E-03]
+        me.vel[:, 8] = [+2.568651772461E-03, +1.681832388267E-03, +6.245613982833E-04]
+        me.vel[:, 9] = [+3.034112963576E-03, -1.111317562971E-03, -1.261841468083E-03]
 
         # masses relative to the sun taken from
         # https://en.wikipedia.org/wiki/Planetary_mass#Values_from_the_DE405_ephemeris
diff --git a/pySDC/implementations/problem_classes/HeatEquation_ND_FD_forced_periodic.py b/pySDC/implementations/problem_classes/HeatEquation_ND_FD_forced_periodic.py
@@ -96,7 +96,7 @@ def __get_A(N, nu, dx, ndim, order):
             stencil = [1 / 90, -3 / 20, 3 / 2, -49 / 18, 3 / 2, -3 / 20, 1 / 90]
             zero_pos = 4
         elif order == 8:
-            stencil = [-1 / 560, 8 / 315, -1 / 5, 8 / 5, -205 / 72, 8 / 5, -1 / 5, 8 / 315,	-1 / 560]
+            stencil = [-1 / 560, 8 / 315, -1 / 5, 8 / 5, -205 / 72, 8 / 5, -1 / 5, 8 / 315, -1 / 560]
             zero_pos = 5
         else:
             raise ProblemError(f'wrong order given, has to be 2, 4, 6, or 8, got {order}')
diff --git a/pySDC/implementations/problem_classes/PenningTrap_3D.py b/pySDC/implementations/problem_classes/PenningTrap_3D.py
@@ -109,7 +109,7 @@ def u_init(self):
 
         u = self.dtype_u(((3, N), self.init[1], self.init[2]))
 
-        if u0[2][0] is not 1 or u0[3][0] is not 1:
+        if u0[2][0] != 1 or u0[3][0] != 1:
             raise ProblemError('so far only q = m = 1 is implemented')
 
         # set first particle to u0
diff --git a/pySDC/implementations/problem_classes/acoustic_helpers/standard_integrators.py b/pySDC/implementations/problem_classes/acoustic_helpers/standard_integrators.py
@@ -1,5 +1,5 @@
 import math
-from decimal import Decimal
+from decimal import Decimal, getcontext
 
 import numpy as np
 import scipy.sparse as sp
diff --git a/pySDC/implementations/problem_classes/boussinesq_helpers/standard_integrators.py b/pySDC/implementations/problem_classes/boussinesq_helpers/standard_integrators.py
@@ -1,5 +1,5 @@
 import math
-from decimal import Decimal
+from decimal import Decimal, getcontext
 
 import numpy as np
 import scipy.sparse as sp
diff --git a/pySDC/implementations/sweeper_classes/imex_1st_order.py b/pySDC/implementations/sweeper_classes/imex_1st_order.py
@@ -44,7 +44,6 @@ def integrate(self):
 
         # get current level and problem description
         L = self.level
-        P = L.prob
 
         me = []
 
diff --git a/pySDC/implementations/transfer_classes/TransferMesh.py b/pySDC/implementations/transfer_classes/TransferMesh.py
@@ -124,11 +124,11 @@ def __init__(self, fine_prob, coarse_prob, params):
                         Rspace.append(restr_factor * Pspace[-1].T)
 
                     else:
-
-                        Rspace.append(restr_factor *
-                                      th.interpolation_matrix_1d(fine_grid, coarse_grid, k=self.params.rorder,
-                                                                 periodic=self.params.periodic,
-                                                                 equidist_nested=self.params.equidist_nested).T)
+                        mat = th.interpolation_matrix_1d(fine_grid, coarse_grid,
+                                                         k=self.params.rorder,
+                                                         periodic=self.params.periodic,
+                                                         equidist_nested=self.params.equidist_nested).T
+                        Rspace.append(restr_factor * mat)
 
             # kronecker 1-d operators for n-d
             self.Pspace = Pspace[0]
diff --git a/pySDC/implementations/transfer_classes/TransferMesh_MPIFFT.py b/pySDC/implementations/transfer_classes/TransferMesh_MPIFFT.py
@@ -2,7 +2,6 @@
 from pySDC.core.SpaceTransfer import space_transfer
 from pySDC.implementations.datatype_classes.mesh import mesh, imex_mesh
 from mpi4py_fft import PFFT, newDistArray
-import numpy as np
 
 
 class fft_to_fft(space_transfer):
diff --git a/pySDC/implementations/transfer_classes/TransferMesh_NoCoarse.py b/pySDC/implementations/transfer_classes/TransferMesh_NoCoarse.py
@@ -1,6 +1,3 @@
-
-import scipy.sparse as sp
-
 from pySDC.core.Errors import TransferError
 from pySDC.core.SpaceTransfer import space_transfer
 from pySDC.implementations.datatype_classes.mesh import mesh, imex_mesh
diff --git a/pySDC/projects/AllenCahn_Bayreuth/AllenCahn_monitor.py b/pySDC/projects/AllenCahn_Bayreuth/AllenCahn_monitor.py
@@ -1,5 +1,4 @@
 import numpy as np
-import json
 from mpi4py import MPI
 from pySDC.core.Hooks import hooks
 
diff --git a/pySDC/projects/AllenCahn_Bayreuth/run_simple_forcing_verification.py b/pySDC/projects/AllenCahn_Bayreuth/run_simple_forcing_verification.py
@@ -167,7 +167,7 @@ def run_simulation(name='', spectral=None, nprocs_space=None):
         out = f'Time to solution: {timing[0][1]:.4f} sec.'
         print(out)
 
-        out = f'...Done <---------\n'
+        out = '...Done <---------\n'
         print(out)
 
         # Testing the output
diff --git a/pySDC/projects/AllenCahn_Bayreuth/run_temp_forcing_reference.py b/pySDC/projects/AllenCahn_Bayreuth/run_temp_forcing_reference.py
@@ -121,7 +121,7 @@ def run_simulation(name='', spectral=None, nprocs_space=None):
         out = f'Time to solution: {timing[0][1]:.4f} sec.'
         print(out)
 
-        out = f'...Done <---------\n'
+        out = '...Done <---------\n'
         print(out)
 
 
diff --git a/pySDC/projects/AllenCahn_Bayreuth/run_temp_forcing_verification.py b/pySDC/projects/AllenCahn_Bayreuth/run_temp_forcing_verification.py
@@ -142,7 +142,7 @@ def run_simulation(name='', spectral=None, nprocs_time=None, nprocs_space=None,
         else:
             err = abs(array - uend)
 
-        out = f'...Done <---------\n'
+        out = '...Done <---------\n'
         print(out)
 
         return err
diff --git a/pySDC/projects/AllenCahn_Bayreuth/visualize.py b/pySDC/projects/AllenCahn_Bayreuth/visualize.py
@@ -37,7 +37,7 @@ def plot_data(name=''):
         plt.colorbar()
         plt.title(f"Time: {obj['time']:6.4f}")
 
-        plt.savefig(f'data/{name}_{index}.png',  bbox_inches='tight')
+        plt.savefig(f'data/{name}_{index}.png', bbox_inches='tight')
         plt.close()
 
 
diff --git a/pySDC/projects/AllenCahn_Bayreuth/visualize_temp.py b/pySDC/projects/AllenCahn_Bayreuth/visualize_temp.py
@@ -56,7 +56,7 @@ def plot_data(path='./data', name='', output='.'):
 
         grid.cbar_axes[0].colorbar(im)
 
-        plt.savefig(f'{output}/{name}_{index}.png',  bbox_inches='tight')
+        plt.savefig(f'{output}/{name}_{index}.png', bbox_inches='tight')
         plt.close()
 
 
@@ -77,7 +77,6 @@ def make_movie(path='./data', name='', output='.'):
     data_files = sorted(glob.glob(f'{path}/{name}_*.dat'))
 
     img_list = []
-    c = 0
     for json_file, data_file in zip(json_files, data_files):
         with open(json_file, 'r') as fp:
             obj = json.load(fp)
diff --git a/pySDC/projects/AsympConv/PFASST_conv_Linf.py b/pySDC/projects/AsympConv/PFASST_conv_Linf.py
@@ -297,7 +297,7 @@ def plot_results(cwd=''):
 
         # save plot, beautify
         fname = 'data/conv_test_niter_Linf_' + type + '.png'
-        plt.savefig(fname,  bbox_inches='tight')
+        plt.savefig(fname, bbox_inches='tight')
 
         assert os.path.isfile(fname), 'ERROR: plotting did not create file'
 
diff --git a/pySDC/projects/AsympConv/PFASST_conv_tests.py b/pySDC/projects/AsympConv/PFASST_conv_tests.py
@@ -302,7 +302,7 @@ def plot_results(nsweeps):
     # plt.show()
     # save plot, beautify
     fname = 'data/conv_test_niter_NS' + str(nsweeps) + '.png'
-    plt.savefig(fname,  bbox_inches='tight')
+    plt.savefig(fname, bbox_inches='tight')
 
     assert os.path.isfile(fname), 'ERROR: plotting did not create file'
 
diff --git a/pySDC/projects/AsympConv/conv_test_to0.py b/pySDC/projects/AsympConv/conv_test_to0.py
@@ -78,7 +78,7 @@ def compute_and_plot_specrad(Nnodes, lam):
         fname = 'data/smoother_specrad_to0_L64+256_M' + str(Nnodes) + 'LU_imag.png'
     else:
         fname = 'data/smoother_specrad_to0_L64+256_M' + str(Nnodes) + 'LU_real.png'
-    plt.savefig(fname,  transparent=True, bbox_inches='tight')
+    plt.savefig(fname, transparent=True, bbox_inches='tight')
 
 
 if __name__ == "__main__":
diff --git a/pySDC/projects/AsympConv/conv_test_toinf.py b/pySDC/projects/AsympConv/conv_test_toinf.py
@@ -78,7 +78,7 @@ def compute_and_plot_specrad(Nnodes, lam):
         fname = 'data/smoother_specrad_toinf_M' + str(Nnodes) + '_LU_imag.png'
     else:
         fname = 'data/smoother_specrad_toinf_M' + str(Nnodes) + '_LU_real.png'
-    plt.savefig(fname,  transparent=True, bbox_inches='tight')
+    plt.savefig(fname, transparent=True, bbox_inches='tight')
 
 
 if __name__ == "__main__":
diff --git a/pySDC/projects/AsympConv/smoother_specrad_heatmap.py b/pySDC/projects/AsympConv/smoother_specrad_heatmap.py
@@ -147,7 +147,7 @@ def compute_and_plot_specrad():
 
         fname = 'data/heatmap_smoother_' + conv_type + '_Nsteps' + str(Nsteps) + '_M' + \
                 str(Nnodes) + '_' + qd_type + '.png'
-        plt.savefig(fname,  transparent=True, bbox_inches='tight')
+        plt.savefig(fname, transparent=True, bbox_inches='tight')
 
 
 if __name__ == "__main__":
diff --git a/pySDC/projects/FastWaveSlowWave/plot_stability.py b/pySDC/projects/FastWaveSlowWave/plot_stability.py
@@ -77,7 +77,7 @@ def compute_stability():
         for j in range(0, N_f):
             lambda_fast = lambda_f[j]
             lambda_slow = lambda_s[i]
-            if K is not 0:
+            if K != 0:
                 lambdas = [lambda_fast, lambda_slow]
                 # LHS, RHS = L.sweep.get_scalar_problems_sweeper_mats(lambdas=lambdas)
                 Mat_sweep = L.sweep.get_scalar_problems_manysweep_mat(nsweeps=K, lambdas=lambdas)
diff --git a/pySDC/projects/FastWaveSlowWave/runmultiscale_acoustic.py b/pySDC/projects/FastWaveSlowWave/runmultiscale_acoustic.py
@@ -123,9 +123,9 @@ def compute_and_plot_solutions():
     fig = plt.figure()
 
     sigma_0 = 0.1
-    k = 7.0 * 2.0 * np.pi
+    # k = 7.0 * 2.0 * np.pi
     x_0 = 0.75
-    x_1 = 0.25
+    # x_1 = 0.25
 
     print('Maximum pressure in SDC: %5.3e' % np.linalg.norm(uend[1, :], np.inf))
     print('Maximum pressure in DIRK: %5.3e' % np.linalg.norm(pnew_dirk, np.inf))
diff --git a/pySDC/projects/Hamiltonian/harmonic_oscillator.py b/pySDC/projects/Hamiltonian/harmonic_oscillator.py
@@ -119,7 +119,7 @@ def run_simulation():
     dill.dump(results, f)
     f.close()
 
-    assert os.path.isfile(fname), 'Run for %s did not create stats file' % prob
+    assert os.path.isfile(fname), 'Run did not create stats file'
 
 
 def show_results(cwd=''):
diff --git a/pySDC/projects/Hamiltonian/simple_problems.py b/pySDC/projects/Hamiltonian/simple_problems.py
@@ -133,7 +133,7 @@ def run_simulation(prob=None):
         num_procs = 100
         maxmeaniter = 5.0
     else:
-        raise NotImplemented('Problem type not implemented, got %s' % prob)
+        raise NotImplementedError('Problem type not implemented, got %s' % prob)
 
     # instantiate the controller
     controller = controller_nonMPI(num_procs=num_procs, controller_params=controller_params,
diff --git a/pySDC/projects/Hamiltonian/solar_system.py b/pySDC/projects/Hamiltonian/solar_system.py
@@ -132,7 +132,7 @@ def run_simulation(prob=None):
         num_procs = 100
         maxmeaniter = 19.0
     else:
-        raise NotImplemented('Problem type not implemented, got %s' % prob)
+        raise NotImplementedError('Problem type not implemented, got %s' % prob)
 
     f = open(prob + '_out.txt', 'w')
     out = 'Running ' + prob + ' problem with %s processors...' % num_procs
@@ -259,7 +259,7 @@ def show_results(prob=None, cwd=''):
         elif ndim == 3:
             ax.plot(pos[n, 0, :], pos[n, 1, :], pos[n, 2, :])
         else:
-            raise NotImplemented('Wrong number of dimensions for plotting, got %s' % ndim)
+            raise NotImplementedError('Wrong number of dimensions for plotting, got %s' % ndim)
 
     fname = 'data/' + prob + '_positions'
     plt_helper.savefig(fname)
diff --git a/pySDC/projects/Performance/run_simple_forcing_benchmark.py b/pySDC/projects/Performance/run_simple_forcing_benchmark.py
@@ -8,7 +8,7 @@
 from pySDC.implementations.sweeper_classes.imex_1st_order import imex_1st_order
 from pySDC.implementations.problem_classes.AllenCahn_MPIFFT import allencahn_imex, allencahn_imex_timeforcing
 from pySDC.implementations.transfer_classes.TransferMesh_MPIFFT import fft_to_fft
-from pySDC.projects.AllenCahn_Bayreuth.AllenCahn_dump import dump
+# from pySDC.projects.AllenCahn_Bayreuth.AllenCahn_dump import dump
 
 # from pySDC.projects.Performance.controller_MPI_scorep import controller_MPI
 
diff --git a/pySDC/projects/RDC/vanderpol_error_test.py b/pySDC/projects/RDC/vanderpol_error_test.py
@@ -146,7 +146,7 @@ def plot_RDC_results(cwd=''):
 
     # save plot as PNG, beautify
     fname = 'data/RDC_errors_vdp.png'
-    plt.savefig(fname,  bbox_inches='tight')
+    plt.savefig(fname, bbox_inches='tight')
 
     assert os.path.isfile(fname), 'ERROR: plot was not created'
 
diff --git a/pySDC/projects/SDC_showdown/SDC_timing_Fisher.py b/pySDC/projects/SDC_showdown/SDC_timing_Fisher.py
@@ -104,7 +104,7 @@ def run_SDC_variant(variant=None, inexact=False):
         description['problem_class'] = petsc_fisher_multiimplicit
         description['sweeper_class'] = multi_implicit
     else:
-        raise NotImplemented('Wrong variant specified, got %s' % variant)
+        raise NotImplementedError('Wrong variant specified, got %s' % variant)
 
     if inexact:
         description['problem_params']['nlsol_maxiter'] = 1
diff --git a/pySDC/projects/SDC_showdown/SDC_timing_GrayScott.py b/pySDC/projects/SDC_showdown/SDC_timing_GrayScott.py
@@ -107,7 +107,7 @@ def run_SDC_variant(variant=None, inexact=False, cwd=''):
         description['problem_class'] = petsc_grayscott_multiimplicit
         description['sweeper_class'] = multi_implicit
     else:
-        raise NotImplemented('Wrong variant specified, got %s' % variant)
+        raise NotImplementedError('Wrong variant specified, got %s' % variant)
 
     if inexact:
         description['problem_params']['lsol_maxiter'] = 2
diff --git a/pySDC/projects/TOMS/AllenCahn_contracting_circle.py b/pySDC/projects/TOMS/AllenCahn_contracting_circle.py
@@ -122,7 +122,7 @@ def run_SDC_variant(variant=None, inexact=False):
         if inexact:
             description['problem_params']['newton_maxiter'] = 1
     else:
-        raise NotImplemented('Wrong variant specified, got %s' % variant)
+        raise NotImplementedError('Wrong variant specified, got %s' % variant)
 
     if inexact:
         out = 'Working on inexact %s variant...' % variant
diff --git a/pySDC/projects/node_failure/animate_convergence.py b/pySDC/projects/node_failure/animate_convergence.py
@@ -70,7 +70,7 @@ def create_animation(cwd=''):
 
     residual = np.zeros((maxiter_full, maxstep - minstep))
     plot = plt.pcolor(residual, cmap=cmap, vmin=vmin, vmax=vmax)
-    text = plt.text(0, 0, '', horizontalalignment='center', verticalalignment='center')
+    plt.text(0, 0, '', horizontalalignment='center', verticalalignment='center')
 
     ticks = np.arange(vmin, vmax + 1, 2)
     tickpos = np.linspace(ticks[0] + 0.5, ticks[-1] - 0.5, len(ticks))
diff --git a/pySDC/projects/node_failure/boussinesq_example.py b/pySDC/projects/node_failure/boussinesq_example.py
@@ -96,7 +96,7 @@ def main(ft_strategies):
         ft.strategy = strategy
 
         # read in reference data from clean run, will provide reproducable locations for faults
-        if strategy is not 'NOFAULT':
+        if strategy != 'NOFAULT':
             reffile = np.load('data/PFASST_BOUSSINESQ_stats_hf_NOFAULT_P16.npz')
             ft.refdata = reffile['hard_stats']
 
diff --git a/pySDC/projects/node_failure/emulate_hard_faults.py b/pySDC/projects/node_failure/emulate_hard_faults.py
diff --git a/pySDC/projects/node_failure/grayscott_example.py b/pySDC/projects/node_failure/grayscott_example.py
diff --git a/pySDC/projects/node_failure/hard_faults_detail.py b/pySDC/projects/node_failure/hard_faults_detail.py
diff --git a/pySDC/projects/node_failure/hard_faults_test.py b/pySDC/projects/node_failure/hard_faults_test.py
diff --git a/pySDC/projects/node_failure/postproc_boussinesq.py b/pySDC/projects/node_failure/postproc_boussinesq.py
diff --git a/pySDC/projects/node_failure/postproc_hard_faults_detail.py b/pySDC/projects/node_failure/postproc_hard_faults_detail.py
diff --git a/pySDC/projects/node_failure/postproc_hard_faults_test.py b/pySDC/projects/node_failure/postproc_hard_faults_test.py
diff --git a/pySDC/projects/parallelSDC/AllenCahn_parallel.py b/pySDC/projects/parallelSDC/AllenCahn_parallel.py
diff --git a/pySDC/projects/parallelSDC/generic_implicit_MPI.py b/pySDC/projects/parallelSDC/generic_implicit_MPI.py
diff --git a/pySDC/projects/parallelSDC/minimization.py b/pySDC/projects/parallelSDC/minimization.py
diff --git a/pySDC/projects/soft_failure/generate_statistics.py b/pySDC/projects/soft_failure/generate_statistics.py
diff --git a/pySDC/tutorial/step_1/B_spatial_accuracy_check.py b/pySDC/tutorial/step_1/B_spatial_accuracy_check.py
diff --git a/pySDC/tutorial/step_1/D_collocation_accuracy_check.py b/pySDC/tutorial/step_1/D_collocation_accuracy_check.py
diff --git a/pySDC/tutorial/step_6/C_MPI_parallelization.py b/pySDC/tutorial/step_6/C_MPI_parallelization.py

Original file line number	Diff line number	Diff line change
`@@ -1,5 +1,4 @@`
`1`	`1`	`import numpy as np`
`2`		`-from mpi4py import MPI`
`3`	`2`	`from mpi4py_fft import PFFT`
`4`	`3`
`5`	`4`	`from pySDC.core.Errors import ParameterError, ProblemError`