Fixes

Author: Thomas Baumann

pySDC/projects/GPU/analysis_scripts/RBC3D_order.py

@@ -9,8 +9,12 @@
 
 step_sizes = {
     'RBC3DG4R4Ra1e5': [8e-2, 4e-2, 2e-2, 1e-2, 5e-3],
-    'RBC3DG4R4SDC23Ra1e5': [8e-2, 4e-2, 2e-2, 1e-2, 5e-3, 2.5e-3],
-    'RBC3DG4R4SDC34Ra1e5': [8e-2, 4e-2, 2e-2, 1e-2, 5e-3],
+    # 'RBC3DG4R4SDC23Ra1e5': [8e-2, 4e-2, 2e-2, 1e-2, 5e-3, 2.5e-3],
+    # 'RBC3DG4R4SDC34Ra1e5': [8e-2, 4e-2, 2e-2, 1e-2, 5e-3],
+    'RBC3DG4R4SDC23Ra1e5': [5e-3 * 2**i for i in range(5)],
+    'RBC3DG4R4SDC44Ra1e5': [5e-3 * 2**i for i in range(5)],
+    # 'RBC3DG4R4SDC44Ra1e5': [1e-3 * 2**i for i in range(8)],
+    'RBC3DG4R4SDC34Ra1e5': [1e-3 * 2**i for i in range(10)],
     'RBC3DG4R4RKRa1e5': [8e-2, 4e-2, 2e-2, 1e-2, 5e-3, 2.5e-3],
     # 'RBC3DG4R4EulerRa1e5': [8e-2, 4e-2, 2e-2, 1e-2, 5e-3],
     'RBC3DG4R4EulerRa1e5': [1e-3 * 2**i for i in range(5)],
@@ -30,7 +34,7 @@ def get_path(args):
 
 def compute_errors(args, dts, Tend):
     errors = {'u': [], 'v': [], 'w': [], 'T': [], 'p': [], 'dt': []}
-    prob = RayleighBenard3D(nx=4, ny=4, nz=4)
+    prob = RayleighBenard3D(nx=4, ny=4, nz=4, comm=MPI.COMM_SELF)
 
     dts = np.sort(dts)[::-1]
     ref = run(args, dts[-1], Tend)
@@ -104,7 +108,7 @@ def run(args, dt, Tend):
     config.Tend = n_freefall_times.get(type(config).__name__, 3)
 
     ic_config_name = type(config).__name__
-    for name in ['RK', 'Euler', 'O3', 'O4', 'SDC23', 'SDC34']:
+    for name in ['RK', 'Euler', 'O3', 'O4', 'SDC23', 'SDC34', 'SDC44']:
         ic_config_name = ic_config_name.replace(name, '')
     ic_config = get_config({**args, 'config': ic_config_name})
     config.ic_config['config'] = type(ic_config)

pySDC/projects/GPU/configs/RBC3D_configs.py

@@ -221,6 +221,15 @@ def get_description(self, *args, **kwargs):
         return desc
 
 
+class RBC3DM4K4(RBC3Dverification):
+
+    def get_description(self, *args, **kwargs):
+        desc = super().get_description(*args, **kwargs)
+        desc['level_params']['nsweeps'] = 4
+        desc['sweeper_params']['num_nodes'] = 4
+        return desc
+
+
 class RBC3DverificationRK(RBC3Dverification):
 
     def get_description(self, *args, res=-1, dt=-1, **kwargs):
@@ -268,6 +277,13 @@ class RBC3DG4R4SDC34Ra1e5(RBC3DM3K4):
     converged = 50
 
 
+class RBC3DG4R4SDC44Ra1e5(RBC3DM4K4):
+    Tend = 200
+    dt = 6e-2
+    res = 32
+    converged = 50
+
+
 class RBC3DG4R4RKRa1e5(RBC3DverificationRK):
     Tend = 200
     dt = 8e-2