now reports on total number of gmres iterations of sdc and comparison run with dirk-4; plots horizontal cut through end solution for validation

Author: Daniel Ruprecht

examples/boussinesq_2d_imex/HookClass.py

@@ -20,8 +20,11 @@ def __init__(self):
 
         # add figure object for further use
         #        self.fig = plt.figure(figsize=(18,6))
-        self.fig = plt.figure(figsize=(15,5))
-        self.counter = 0
+        #self.fig = plt.figure(figsize=(15,5))
+        #self.counter = 0
+
+    def dump_sweep(self,status):
+        pass
 
     def dump_step(self,status):
         """

examples/boussinesq_2d_imex/ProblemClass.py

@@ -106,7 +106,7 @@ def solve_system(self,rhs,factor,u0,t):
         sol, info = LA.gmres( self.Id - factor*self.M, b, x0=u0.values.flatten(), tol=self.gmres_tol, restart=self.gmres_restart, maxiter=self.gmres_maxiter, callback=cb)
         # If this is a dummy call with factor==0.0, do not log because it should not be counted as a solver call
         if factor!=0.0:
-          print "Number of GMRES iterations: %3i --- Final residual: %6.3e" % ( cb.getcounter(), cb.getresidual() )
+          print "SDC: Number of GMRES iterations: %3i --- Final residual: %6.3e" % ( cb.getcounter(), cb.getresidual() )
           self.logger.add(cb.getcounter())
         me        = mesh(self.nvars)
         me.values = unflatten(sol, 4, self.N[0], self.N[1])
@@ -197,8 +197,3 @@ def u_exact(self,t):
         me.values[2,:,:] = dtheta*np.sin( np.pi*self.zz/H )/( 1.0 + np.square(self.xx - x_c)/(a*a))
         me.values[3,:,:] = 0.0*self.xx
         return me
-        
-    def report_log(self):
-      print "Number of calls to implicit solver: %5i" % self.logger.solver_calls
-      print "Total number of iterations: %5i" % self.logger.iterations
-      print "Average number of iterations per call: %6.3f" % (float(self.logger.iterations)/float(self.logger.solver_calls))