added a couple of useful notes to hook class

Author: Daniel Ruprecht

examples/boussinesq_2d_imex/HookClass.py

@@ -66,4 +66,8 @@ def dump_step(self,status):
           self.counter = self.counter + 1
           plt.pause(0.001)
 
+          # NOTE: Can use ffmpeg to collate images into movie.
+          # HELPFUL WEBSITE: http://hamelot.co.uk/visualization/using-ffmpeg-to-convert-a-set-of-images-into-a-video/
+          # USEFUL COMMAND: ffmpeg -r 25 -i boussinesq-%04d.jpeg -vcodec libx264 -crf 25 test1800.avi
+          # WEBSITE TO CONVERT TO MP4: http://video.online-convert.com/convert-to-mp4
         return None

examples/boussinesq_2d_imex/playground.py

@@ -28,7 +28,7 @@
 
     # This comes as read-in for the level class
     lparams = {}
-    lparams['restol'] = 1E-6
+    lparams['restol'] = 1E-14
 
     swparams = {}
     swparams['collocation_class'] = collclass.CollGaussLobatto
@@ -41,7 +41,9 @@
     # setup parameters "in time"
     t0     = 0
     Tend   = 3000
-    Nsteps =  250
+    Nsteps =  500
+    Tend   = 30
+    Nsteps =  5
     dt = Tend/float(Nsteps)
 
     # This comes as read-in for the problem class
@@ -71,8 +73,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)
@@ -97,10 +99,10 @@
     print('error at time %s: %9.5e' %(Tend,np.linalg.norm(uex.values[2,:,:].flatten()-uend.values[2,:,:].flatten(),np.inf)/np.linalg.norm(
         uex.values.flatten(),np.inf)))
 
-    #P.report_log()
+    P.report_log()
 
     #plt.show()
 
-    # extract_stats = grep_stats(stats,iter=-1,type='residual')
-    # sortedlist_stats = sort_stats(extract_stats,sortby='step')
-    # print(extract_stats,sortedlist_stats)
+    #extract_stats = grep_stats(stats,iter=-1,type='residual')
+    #sortedlist_stats = sort_stats(extract_stats,sortby='step')
+    #print(extract_stats,sortedlist_stats)