improved multiscale example plots with markers

Author: Daniel Ruprecht

examples/acoustic_1d_imex/runmultiscale.py

@@ -131,23 +131,23 @@
     print ('Maximum pressure in RK-IMEX: %5.3e' % np.linalg.norm(pnew_imex, np.inf))
 
     #plt.plot(P.mesh, pnew_tp,  '-', color='c', label='Trapezoidal')
-    if np.linalg.norm(pnew_imex, np.inf)<=2:
-      plt.plot(P.mesh, pnew_imex,  '-', color='c', label='IMEX('+str(rkimex.order)+')')
-    plt.plot(P.mesh, uend.values[1,:], '--', color='b', label='SDC('+str(sparams['maxiter'])+')')
     if dirk.order==2:
-      plt.plot(P.mesh, pnew_bdf, '-', color='r', label='BDF-2')
-    plt.plot(P.mesh, pnew_dirk, color='g', label='DIRK('+str(dirk.order)+')')
+      plt.plot(P.mesh, pnew_bdf, 'd-', color='c', label='BDF-2',markevery=(50,75))
+    p_slow = np.exp(-np.square( np.mod( P.mesh-pparams['cadv']*Tend, 1.0 ) -x_0 )/(sigma_0*sigma_0))
+    plt.plot(P.mesh, p_slow, '--', color='k', markersize=fs-2, label='Slow mode', dashes=(10,2))
+    if np.linalg.norm(pnew_imex, np.inf)<=2:
+      plt.plot(P.mesh, pnew_imex,  '+-', color='r', label='IMEX('+str(rkimex.order)+')',markevery=(1,75),mew=1.0)
+    plt.plot(P.mesh, uend.values[1,:], 'o-', color='b', label='SDC('+str(sparams['maxiter'])+')',markevery=(25,75))
+    plt.plot(P.mesh, pnew_dirk, '-', color='g', label='DIRK('+str(dirk.order)+')')
     #plt.plot(P.mesh, uex.values[1,:],  '+', color='r', label='p (exact)')
     #plt.plot(P.mesh, uend.values[1,:], '-', color='b', linewidth=2.0, label='p (SDC)')
 
-    p_slow = np.exp(-np.square( np.mod( P.mesh-pparams['cadv']*Tend, 1.0 ) -x_0 )/(sigma_0*sigma_0))
-    plt.plot(P.mesh, p_slow, '+', color='k', markersize=fs-2, label='Slow mode', markevery=10)
     plt.xlabel('x', fontsize=fs, labelpad=0)
     plt.ylabel('Pressure', fontsize=fs, labelpad=0)
     fig.gca().set_xlim([0, 1.0])
     fig.gca().set_ylim([-0.5, 1.1])
     fig.gca().tick_params(axis='both', labelsize=fs)
-    plt.legend(loc='upper left', fontsize=fs, prop={'size':fs})
+    plt.legend(loc='upper left', fontsize=fs, prop={'size':fs}, handlelength=3)
     fig.gca().grid()
     #plt.show()
     filename = 'sdc-fwsw-multiscale-K'+str(sparams['maxiter'])+'-M'+str(description['num_nodes'])+'.pdf'