tweaked parameter for convergence plots; use two sine modes as initial data instead of only one

Author: Daniel Ruprecht

examples/acoustic_1d_imex/ProblemClass_conv.py

@@ -20,7 +20,7 @@
 from buildWave1DMatrix import getWave1DMatrix, getWave1DAdvectionMatrix
 
 def u_initial(x, k):
-    return np.sin(k*2.0*np.pi*x)
+    return np.sin(k*2.0*np.pi*x) + np.sin(2.0*np.pi*x)
 
 class acoustic_1d_imex(ptype):
     """

examples/acoustic_1d_imex/plotconvdata.py

@@ -45,17 +45,17 @@
   plt.loglog(nsteps_plot[ii,:], error_plot[ii,:], shape[ii], markersize=fs, color=color[ii], label='p='+str(int(order_plot[ii])))
 
 
-plt.legend(loc='upper right', fontsize=fs, prop={'size':fs})
+plt.legend(loc='lower left', fontsize=fs, prop={'size':fs})
 plt.xlabel('Number of time steps', fontsize=fs)
 plt.ylabel('Relative error', fontsize=fs, labelpad=2)
 plt.xlim([0.9*np.min(nsteps_plot), 1.1*np.max(nsteps_plot)])
-plt.ylim([1e-7, 1e1])
-plt.yticks([1e-7, 1e-6, 1e-5, 1e-4, 1e-3, 1e-2, 1e-1, 1e0, 1e1],fontsize=fs)
-plt.xticks([25, 50, 100], fontsize=fs)
+plt.ylim([1e-5, 1e0])
+plt.yticks([1e-5, 1e-4, 1e-3, 1e-2, 1e-1, 1e0],fontsize=fs)
+plt.xticks([30, 40, 60, 80], fontsize=fs)
 plt.gca().get_xaxis().get_major_formatter().labelOnlyBase = False
 plt.gca().get_xaxis().set_major_formatter(ScalarFormatter())
 plt.show()
-filename = 'sdc_fwsw_convergence.pdf'
+filename = 'sdc-fwsw-convergence.pdf'
 fig.savefig(filename,bbox_inches='tight')
 call(["pdfcrop", filename, filename])
 

examples/acoustic_1d_imex/runconvergence.py

@@ -31,15 +31,11 @@
 
     # This comes as read-in for the problem class
     pparams = {}
-    pparams['nvars']     = [(2,400)]
-    pparams['cadv']      = 0.05
+    pparams['nvars']     = [(2,300)]
+    pparams['cadv']      = 0.1
     pparams['cs']        = 1.00
     pparams['order_adv'] = 5
-    pparams['waveno']    = 2
-
-    # This comes as read-in for the transfer operations
-    tparams = {}
-    tparams['finter'] = True
+    pparams['waveno']    = 5
 
     # Fill description dictionary for easy hierarchy creation
     description = {}
@@ -53,7 +49,7 @@
     description['level_params']      = lparams
     description['hook_class']        = plot_solution
 
-    Nsteps = [15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80]
+    Nsteps = [30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80]
 
     for order in [2, 3, 4]:
 

examples/acoustic_1d_imex/runitererror.py

@@ -46,11 +46,11 @@
 
       # This comes as read-in for the problem class
       pparams = {}
-      pparams['nvars']     = [(2,250)]
-      pparams['cadv']      = 0.05
+      pparams['nvars']     = [(2,300)]
+      pparams['cadv']      = 0.1
       pparams['cs']        = cs_v[cs_ind]
       pparams['order_adv'] = 5
-      pparams['waveno']    = 1
+      pparams['waveno']    = 5
 
       # This comes as read-in for the transfer operations
       tparams = {}
@@ -72,7 +72,7 @@
       for nodes_ind in np.arange(np.size(nodes_v)):
         # setup parameters "in time"
         t0   = 0
-        Tend = 0.05
+        Tend = 0.025
         description['num_nodes'] = nodes_v[nodes_ind]
 
          # quickly generate block of steps
@@ -119,15 +119,15 @@
       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.legend(loc='upper right', fontsize=fs, prop={'size':fs-2})
     plt.xlabel('Iteration', fontsize=fs)
     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.ylim([0, 1.0])
     plt.yticks(fontsize=fs)
     plt.xticks(fontsize=fs)
     plt.show()
-    filename = 'sdc_fwsw_iteration.pdf'
+    filename = 'sdc-fwsw-iteration.pdf'
     fig.savefig(filename,bbox_inches='tight')
     call(["pdfcrop", filename, filename])