empty dump_sweep function for less verbose output

Author: Daniel Ruprecht

examples/acoustic_1d_imex/HookClass.py

@@ -22,6 +22,8 @@ def __init__(self):
         #        self.fig = plt.figure(figsize=(18,6))
         #self.fig = plt.figure(figsize=(9,9))
 
+    def dump_sweep(self,status):
+        return None
 
     def dump_step(self,status):
         """

examples/acoustic_1d_imex/runitererror.py

@@ -0,0 +1,85 @@
+
+from pySDC import CollocationClasses as collclass
+
+import numpy as np
+
+from ProblemClass_conv import acoustic_1d_imex
+from examples.acoustic_1d_imex.HookClass import plot_solution
+
+from pySDC.datatype_classes.mesh import mesh, rhs_imex_mesh
+from pySDC.sweeper_classes.imex_1st_order import imex_1st_order
+import pySDC.PFASST_stepwise as mp
+from pySDC import Log
+from pySDC.Stats import grep_stats, sort_stats
+
+from matplotlib import pyplot as plt
+
+
+if __name__ == "__main__":
+
+    # set global logger (remove this if you do not want the output at all)
+    logger = Log.setup_custom_logger('root')
+
+    num_procs = 1
+
+    # This comes as read-in for the level class
+    lparams = {}
+    lparams['restol'] = 1E-14
+
+    sparams = {}
+  
+
+    # This comes as read-in for the problem class
+    pparams = {}
+    pparams['nvars']     = [(2,250)]
+    pparams['cadv']      = 0.05
+    pparams['cs']        = 5.0
+    pparams['order_adv'] = 5
+    pparams['waveno']    = 1
+
+    # This comes as read-in for the transfer operations
+    tparams = {}
+    tparams['finter'] = True
+
+    # Fill description dictionary for easy hierarchy creation
+    description = {}
+    description['problem_class']     = acoustic_1d_imex
+    description['problem_params']    = pparams
+    description['dtype_u']           = mesh
+    description['dtype_f']           = rhs_imex_mesh
+    description['collocation_class'] = collclass.CollGaussLobatto
+    description['sweeper_class']     = imex_1st_order
+    description['level_params']      = lparams
+    description['hook_class']        = plot_solution
+    #description['transfer_class'] = mesh_to_mesh_1d
+    #description['transfer_params'] = tparams
+    
+    for order in [3]:
+    
+      # setup parameters "in time"
+      t0   = 0
+      Tend = 0.05
+    
+      if order==2:
+        description['num_nodes'] = 2
+      elif order==3:
+        description['num_nodes'] = 3
+      elif order==4:
+        description['num_nodes'] = 3
+
+      sparams['maxiter'] = 50
+
+       # quickly generate block of steps
+      MS = mp.generate_steps(num_procs,sparams,description)
+    
+      dt = Tend
+
+      # get initial values on finest level
+      P = MS[0].levels[0].prob
+      uinit = P.u_exact(t0)
+
+      # call main function to get things done...
+      uend,stats = mp.run_pfasst(MS, u0=uinit, t0=t0, dt=dt, Tend=Tend)
+      
+      # compute exact solution and compare
+      uex = P.u_exact(Tend)