Stuff

#!!!!!! WARNING: FLAKEHEAVEN FAILED !!!!!!: 

#:

Author: Thomas Baumann

pySDC/implementations/problem_classes/RayleighBenard3D.py

@@ -141,6 +141,10 @@ def __init__(
         U01 = self.get_basis_change_matrix(p_in=0, p_out=1)
         U12 = self.get_basis_change_matrix(p_in=1, p_out=2)
         U02 = self.get_basis_change_matrix(p_in=0, p_out=2)
+        self.eliminate_zeros(S1)
+        self.eliminate_zeros(S2)
+        self.eliminate_zeros(Dz)
+        self.eliminate_zeros(Dzz)
 
         self.Dx = Dx
         self.Dxx = Dxx
@@ -158,13 +162,12 @@ def __init__(
 
         # construct operators
         _D = U02 @ (Dxx + Dyy) + Dzz
-        L_lhs = {
-            'p': {'u': U01 @ Dx, 'v': U01 @ Dy, 'w': Dz},  # divergence free constraint
-            'u': {'p': U02 @ Dx, 'u': -self.nu * _D},
-            'v': {'p': U02 @ Dy, 'v': -self.nu * _D},
-            'w': {'p': U12 @ Dz, 'w': -self.nu * _D, 'T': -U02 @ Id},
-            'T': {'T': -self.kappa * _D},
-        }
+        L_lhs = {}
+        L_lhs['p'] = {'u': U01 @ Dx, 'v': U01 @ Dy, 'w': Dz}  # divergence free constraint
+        L_lhs['u'] = {'p': U02 @ Dx, 'u': -self.nu * _D}
+        L_lhs['v'] = {'p': U02 @ Dy, 'v': -self.nu * _D}
+        L_lhs['w'] = {'p': U12 @ Dz, 'w': -self.nu * _D, 'T': -U02 @ Id}
+        L_lhs['T'] = {'T': -self.kappa * _D}
         self.setup_L(L_lhs)
 
         # mass matrix

pySDC/projects/GPU/analysis_scripts/RBC3d.py

@@ -22,7 +22,15 @@
 config = get_config(args)
 
 desc = config.get_description(**args)
-P = desc['problem_class'](**{**desc['problem_params'], 'spectral_space': False, 'comm': comm})
+P = desc['problem_class'](
+    **{
+        **desc['problem_params'],
+        'spectral_space': False,
+        'comm': comm,
+        'Dirichlet_recombination': False,
+        'left_preconditioner': False,
+    }
+)
 P.setUpFieldsIO()
 xp = P.xp
 

pySDC/projects/GPU/configs/RBC3D_configs.py

@@ -16,6 +16,7 @@ class RayleighBenard3DRegular(Config):
     Tend = 50
     gamma = 1
     res_ratio = 1
+    dealiasing = 3.0 / 2.0
 
     def get_file_name(self):
         res = self.args['res']
@@ -78,6 +79,7 @@ def get_description(self, *args, MPIsweeper=False, res=-1, **kwargs):
         desc['problem_params']['Ly'] = self.gamma
         desc['problem_params']['Lz'] = 1
         desc['problem_params']['heterogeneous'] = True
+        desc['problem_params']['dealiasing'] = self.dealiasing
 
         desc['step_params']['maxiter'] = 3
 
@@ -535,6 +537,15 @@ class RBC3DG4Ra1e5(RBC3DverificationGamma4):
     converged = 50
 
 
+class RBC3DG4R4Ra1e5(RBC3DverificationGamma4):
+    res_ratio = 4
+    Tend = 100
+    dt = 8e-2
+    ic_config = None
+    res = 32
+    converged = 50
+
+
 class RBC3DG4RKRa1e5(RBC3DverificationRKGamma4):
     Tend = 200
     dt = 7e-2
@@ -595,7 +606,34 @@ class RBC3DG4R4Ra1e7(RBC3DverificationGamma4):
     dt = 5e-3
     ic_config = RBC3DG4Ra1e6
     res = 64
-    converged = 11
+    # converged = 11
+
+
+class RBC3DG4R4D2Ra1e7(RBC3DverificationGamma4):
+    res_ratio = 4
+    Tend = 40
+    dt = 5e-3
+    ic_config = RBC3DG4Ra1e6
+    res = 64
+    dealiasing = 1
+
+
+class RBC3DG4R4D4Ra1e7(RBC3DverificationGamma4):
+    res_ratio = 4
+    Tend = 40
+    dt = 5e-3
+    ic_config = RBC3DG4Ra1e6
+    res = 64
+    dealiasing = 2
+
+
+class RBC3DG4R4D42Ra1e7(RBC3DverificationGamma4):
+    res_ratio = 4
+    Tend = 40
+    dt = 5e-3
+    ic_config = RBC3DG4R4Ra1e6
+    res = 64
+    dealiasing = 2
 
 
 class RBC3DG4RKRa1e7(RBC3DverificationRKGamma4):
@@ -606,6 +644,14 @@ class RBC3DG4RKRa1e7(RBC3DverificationRKGamma4):
     res = 96
 
 
+class RBC3DG4R4RKRa1e7(RBC3DverificationRKGamma4):
+    res_ratio = 4
+    Tend = 100
+    dt = 2e-2
+    ic_config = RBC3DG4Ra1e6
+    res = 96
+
+
 class RBC3DG4Ra1e8(RBC3DverificationGamma4):
     # 8e-3 crashes at t=2
     Tend = 100
@@ -624,7 +670,7 @@ class RBC3DG4RKRa1e8(RBC3DverificationRKGamma4):
     converged = 20
 
 
-class RBC3DG4R4Ra1e8(RBC3DverificationGamma4):
+class RBC3DG4Ra1e8(RBC3DverificationGamma4):
     # nz=128: dt=5e-3 crash at t=0.9
     Tend = 60
     dt = 6e-3
@@ -633,13 +679,31 @@ class RBC3DG4R4Ra1e8(RBC3DverificationGamma4):
     # converged=20
 
 
+class RBC3DG4R42Ra1e8(RBC3DverificationGamma4):
+    Tend = 60
+    res_ratio = 4
+    dt = 6e-3
+    ic_config = None
+    res = 96
+    # converged=20
+
+
 class RBC3DG4R4Ra1e8(RBC3DverificationGamma4):
     res_ratio = 4
     Tend = 60
-    dt = 6e-3
+    dt = 5e-3
     ic_config = RBC3DG4R4Ra1e7
     res = 96
-    # converged=20
+    # converged=14
+
+
+class RBC3DG4R4RKRa1e8(RBC3DverificationRKGamma4):
+    res_ratio = 4
+    Tend = 40
+    dt = 3e-3  # limit
+    ic_config = RBC3DG4R4Ra1e7
+    res = 256
+    # converged = 20
 
 
 # class RBC3DG4Ra1e8(RBC3DverificationGamma4):