fix

Author: pancetta

examples/boussinesq_2d_imex/ProblemClass.py

@@ -104,7 +104,6 @@ def solve_system(self,rhs,factor,u0,t):
 
         b         = rhs.values.flatten()
         cb        = Callback()
-        print(self.gmres_tol)
 
         sol, info = LA.gmres( self.Id - factor*self.M, b, x0=u0.values.flatten(), tol=self.gmres_tol, restart=self.gmres_restart, maxiter=self.gmres_maxiter, callback=cb)
         # If this is a dummy call with factor==0.0, do not log because it should not be counted as a solver call

examples/boussinesq_2d_imex/rungmrescounter.py

@@ -26,7 +26,7 @@
 if __name__ == "__main__":
 
     # set global logger (remove this if you do not want the output at all)
-    logger = Log.setup_custom_logger('root')
+    # logger = Log.setup_custom_logger('root')
 
     num_procs = 1
 
@@ -37,7 +37,7 @@
     swparams = {}
     swparams['collocation_class'] = collclass.CollGaussLegendre
     swparams['num_nodes'] = 3
-    swparams['do_LU'] = True
+    swparams['do_LU'] = False
 
     sparams = {}
     sparams['maxiter'] = 4
@@ -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
@@ -64,7 +64,7 @@
     pparams['gmres_maxiter'] = [500]
     pparams['gmres_restart'] = [10]
     pparams['gmres_tol_limit'] = [1e-3]
-    pparams['gmres_tol_factor'] = [0.1]
+    pparams['gmres_tol_factor'] = [0.01]
 
     # This comes as read-in for the transfer operations
     tparams = {}