added output of CFL numbers

Author: Daniel Ruprecht

examples/acoustic_1d_imex/runconvergence.py

@@ -54,7 +54,7 @@
     #description['transfer_class'] = mesh_to_mesh_1d
     #description['transfer_params'] = tparams
 
-    Nsteps = [20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80]
+    Nsteps = [15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80]
 
     for order in [2, 3, 4]:
 
@@ -87,6 +87,10 @@
         # get initial values on finest level
         P = MS[0].levels[0].prob
         uinit = P.u_exact(t0)
+        if ii==0:
+          print "Time step: %4.2f" % dt
+          print "Fast CFL number: %4.2f" % (pparams['cs']*dt/P.dx) 
+          print "Slow CFL number: %4.2f" % (pparams['cadv']*dt/P.dx) 
 
         # call main function to get things done...
         uend,stats = mp.run_pfasst(MS, u0=uinit, t0=t0, dt=dt, Tend=Tend)

examples/acoustic_1d_imex/runitererror.py

@@ -83,6 +83,9 @@
         # get initial values on finest level
         P = MS[0].levels[0].prob
         uinit = P.u_exact(t0)
+        
+        print "Fast CFL number: %4.2f" % (pparams['cs']*dt/P.dx) 
+        print "Slow CFL number: %4.2f" % (pparams['cadv']*dt/P.dx) 
 
         # call main function to get things done...
         uend,stats = mp.run_pfasst(MS, u0=uinit, t0=t0, dt=dt, Tend=Tend)
@@ -113,12 +116,12 @@
     for ii in range(0,np.size(cs_v)):
       x = np.arange(1,lastiter[ii,0])
       y = convrate[ii, 0, 0:lastiter[ii,0]-1]
-      plt.plot(x, y, shape[ii], markersize=fs-2, color=color[ii], label=r'$c_{s}$=%4.2f' % cs_v[ii])       
+      plt.plot(x, y, shape[ii], markersize=fs-2, color=color[ii], label=r'$C_{\rm fast}$=%4.2f' % (cs_v[ii]*dt/P.dx))       
       #plt.plot(x, 0.0*y+avg_convrate[ii,0], '--', color=color[ii])
 
     plt.legend(loc='upper right', fontsize=fs, prop={'size':fs})
     plt.xlabel('Iteration', fontsize=fs)
-    plt.ylabel('Convergence rate', fontsize=fs, labelpad=2)
+    plt.ylabel(r'$|| r^{k+1} ||_{\infty}/|| r^k ||_{\infty}$', fontsize=fs, labelpad=2)
     plt.xlim([0, sparams['maxiter']])
     plt.ylim([0, 0.8])
     plt.yticks(fontsize=fs)