Some fixes

apps/showcases/PhaseFieldAllenCahn/GPU/multiphase_codegen.py

@@ -84,7 +84,7 @@ with CodeGeneration() as ctx:
     method_phase = create_lb_method(lbm_config=lbm_config_phase)
 
     lbm_config_hydro = LBMConfig(stencil=stencil_hydro, method=Method.CENTRAL_MOMENT, compressible=False,
-                                 weighted=True, relaxation_rates=[omega, ] * stencil_hydro.Q,
+                                 weighted=True, relaxation_rate=omega,
                                  force=sp.symbols(f"F_:{stencil_hydro.D}"),
                                  output={'velocity': u})
     method_hydro = create_lb_method(lbm_config=lbm_config_hydro)

apps/showcases/PhaseFieldAllenCahn/GPU/taylor_bubble.py

@@ -58,23 +58,23 @@ class Scenario:
     def __init__(self, cuda_enabled_mpi=False):
         self.density_liquid = 1.0
         self.reference_time = 16000
-        self.reference_length = 64
+        self.reference_length = 128
         d1 = 0.0254  # (0.0508, 0.0381, 0.0254)
         d2 = 0.0127  # (0.0254, 0.0127, 0.0127)
-        self.interface_width = 4
+        self.interface_width = 5
         self.mobility = 0.05
 
         num_processes = wlb.mpi.numProcesses()
 
         # output frequencies
-        self.vtkWriteFrequency = 100 # self.reference_time // 25
+        self.vtkWriteFrequency = 500 # self.reference_time // 25
         self.dbWriteFrequency = 0 # self.reference_time // 25
         self.meshWriteFrequency = self.reference_time
         self.pngWriteFrequency = self.reference_time
 
         # simulation parameters
         self.diameter = self.reference_length
-        self.timesteps = 1000 # self.reference_time #  * 15 + 1
+        self.timesteps = 5001 # self.reference_time #  * 15 + 1
         self.cells = (self.diameter, (self.reference_length * 15) // num_processes, self.diameter)
         self.blocks = (1, num_processes, 1)
         self.periodic = (0, 0, 0)