now prints out inf norm of final pressure to clearly indicate unstable solutions

Author: Daniel Ruprecht

examples/acoustic_1d_imex/runmultiscale.py

@@ -126,6 +126,10 @@
     x_0     = 0.75
     x_1     = 0.25
 
+    print ('Maximum pressure in SDC: %5.3e' % np.linalg.norm(uend.values[1,:], np.inf))
+    print ('Maximum pressure in DIRK: %5.3e' % np.linalg.norm(pnew_dirk, np.inf))
+    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')
     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'])+')')

examples/boussinesq_2d_imex/rungmrescounter.py

@@ -52,8 +52,8 @@
 
     # This comes as read-in for the problem class
     pparams = {}
-    #pparams['nvars']    = [(4,450,30)]
-    pparams['nvars']    = [(4,300,30)]
+    pparams['nvars']    = [(4,450,30)]
+    #pparams['nvars']    = [(4,300,30)]
     pparams['u_adv']    = 0.02
     pparams['c_s']      = 0.3
     pparams['Nfreq']    = 0.01
@@ -63,7 +63,7 @@
     pparams['order_upw'] = [5]
     pparams['gmres_maxiter'] = [500]
     pparams['gmres_restart'] = [10]
-    pparams['gmres_tol']     = [1e-5]
+    pparams['gmres_tol']     = [1e-3]
 
     # This comes as read-in for the transfer operations
     tparams = {}