Added Euler config

Author: Thomas Baumann

pySDC/projects/GPU/analysis_scripts/RBC3D_stability.py

@@ -7,6 +7,7 @@
     'RBC3DG4Ra1e5': [3, 1e0, 1e-1, 8e-2, 6e-2],
     'RBC3DG4RKRa1e5': [1e3, 5, 4, 1e0, 1e-1, 8e-2, 6e-2],
     'RBC3DG4R4Ra1e5': [1e0, 1e-1, 9e-2, 8e-2, 7e-2, 6e-2, 5e-2],
+    'RBC3DG4R4RKRa1e5': [1e0, 1e-1, 9e-2, 8e-2, 7e-2, 6e-2, 5e-2],
 }
 n_freefall_times = {}
 

pySDC/projects/GPU/configs/RBC3D_configs.py

@@ -490,11 +490,26 @@ def get_description(self, *args, res=-1, dt=-1, **kwargs):
         return desc
 
 
+class RBC3DverificationEuler(RBC3DverificationRK):
+
+    def get_description(self, *args, res=-1, dt=-1, **kwargs):
+        from pySDC.implementations.sweeper_classes.Runge_Kutta import IMEXEuler
+
+        desc = super().get_description(*args, res=res, dt=dt, **kwargs)
+        desc['sweeper_class'] = IMEXEuler
+        return desc
+
+
 class RBC3DverificationRKGamma4(RBC3DverificationRK):
     gamma = 4
     res_ratio = 2
 
 
+class RBC3DverificationEulerG4R4(RBC3DverificationEuler):
+    gamma = 4
+    res_ratio = 4
+
+
 class RBC3DRa1e4(RBC3Dverification):
     converged = 60
     dt = 1.0
@@ -539,7 +554,24 @@ class RBC3DG4Ra1e5(RBC3DverificationGamma4):
 
 class RBC3DG4R4Ra1e5(RBC3DverificationGamma4):
     res_ratio = 4
-    Tend = 100
+    Tend = 200
+    dt = 8e-2
+    ic_config = None
+    res = 32
+    converged = 50
+
+
+class RBC3DG4R4RKRa1e5(RBC3DverificationRKGamma4):
+    res_ratio = 4
+    Tend = 200
+    dt = 8e-2
+    ic_config = None
+    res = 32
+    converged = 50
+
+
+class RBC3DG4R4EulerRa1e5(RBC3DverificationEulerG4R4):
+    Tend = 200
     dt = 8e-2
     ic_config = None
     res = 32

pySDC/projects/GPU/configs/base_config.py

@@ -207,7 +207,7 @@ def merge_all_stats(self, controller):
                 with open(self.get_stats_path(index=i), 'rb') as file:
                     _stats = pickle.load(file)
                     stats = {**stats, **_stats}
-            except FileNotFoundError:
+            except (FileNotFoundError, EOFError):
                 print(f'Warning: No stats found at path {self.get_stats_path(index=i)}')
 
         stats = {**stats, **controller.return_stats()}