some comments added

Author: Daniel Ruprecht

examples/SWFW/plot_stability.py

@@ -16,10 +16,12 @@
     pparams['u0'] = 1.0
     swparams = {}
     swparams['collocation_class'] = collclass.CollGaussLobatto
-    swparams['num_nodes'] = 3
+    swparams['num_nodes'] = 9
     K = 1
 
-
+    #
+    # ...this is functionality copied from test_imexsweeper. Ideally, it should be available in one place.
+    #
     step = stepclass.step(params={})
     L = lvl.level(problem_class=swfw_scalar, problem_params=pparams, dtype_u=mesh, dtype_f=rhs_imex_mesh, sweeper_class=imex, sweeper_params=swparams, level_params={}, hook_class=hookclass.hooks, id="stability")
     step.register_level(L)

examples/SWFW/plot_stifflimit_specrad.py

@@ -64,6 +64,7 @@
     plt.plot(nodes_v, specrad[1,:], 'bo-', markersize=fs-2, label=r'$\lambda_{\rm fast} = %2.0f i$' % problem.lambda_f[0].imag)
     plt.xlabel(r'Number of nodes $M$', fontsize=fs)
     plt.ylabel(r'Spectral radius  $\sigma\left( \mathbf{E} \right)$', fontsize=fs, labelpad=2)
+    #plt.title(r'$\Delta t \left| \lambda_{\rm slow} \right|$ = %2.1f' % step.status.dt*abs(problem.lambda_s[0]), fontsize=fs)
     plt.legend(loc='lower right', fontsize=fs, prop={'size':fs})
     plt.xlim([np.min(nodes_v), np.max(nodes_v)])
     plt.ylim([0, 1.0])