Merge pull request #72 from danielru/fix/playground_bouss

Fix playground in Boussinesq example

Author: pancetta

examples/boussinesq_2d_imex/playground.py

@@ -24,7 +24,7 @@
     # set global logger (remove this if you do not want the output at all)
     logger = Log.setup_custom_logger('root')
 
-    num_procs = 8
+    num_procs = 1
 
     # This comes as read-in for the level class
     lparams = {}
@@ -33,7 +33,7 @@
     swparams = {}
     swparams['collocation_class'] = collclass.CollGaussLobatto
     swparams['num_nodes'] = 3
-    swparams['do_LU'] = True
+    swparams['do_LU'] = False
 
     sparams = {}
     sparams['maxiter'] = 12
@@ -48,17 +48,18 @@
 
     # This comes as read-in for the problem class
     pparams = {}
-    pparams['nvars']    = [(4,450,30), (4,450,30)]
+    pparams['nvars']    = [(4,450,30)]
     pparams['u_adv']    = 0.02
     pparams['c_s']      = 0.3
     pparams['Nfreq']    = 0.01
     pparams['x_bounds'] = [(-150.0, 150.0)]
     pparams['z_bounds'] = [(   0.0,  10.0)]
-    pparams['order']    = [4, 2] # [fine_level, coarse_level]
-    pparams['order_upw'] = [5, 1]
-    pparams['gmres_maxiter'] = [50, 50]
-    pparams['gmres_restart'] = [20, 20]
-    pparams['gmres_tol']     = [1e-8, 1e-8]
+    pparams['order']    = [4] # [fine_level, coarse_level]
+    pparams['order_upw'] = [5]
+    pparams['gmres_maxiter'] = [50]
+    pparams['gmres_restart'] = [20]
+    pparams['gmres_tol_limit'] = [1e-5]
+    pparams['gmres_tol_factor'] = [0.1]
 
     # This comes as read-in for the transfer operations
     tparams = {}
@@ -74,8 +75,8 @@
     description['sweeper_class']     = imex_1st_order
     description['level_params']      = lparams
     description['hook_class']        = plot_solution
-    description['transfer_class']    = mesh_to_mesh_2d
-    description['transfer_params']   = tparams
+  #  description['transfer_class']    = mesh_to_mesh_2d
+  #  description['transfer_params']   = tparams
 
     # quickly generate block of steps
     MS = mp.generate_steps(num_procs,sparams,description)