Fix signs for hyperdiffusion coefficient

tum-pbs · May 24, 2024 · 6c57f93 · 6c57f93
1 parent 6b6f275
commit 6c57f93
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/apebench/scenarios/physical/_special_linear.py b/apebench/scenarios/physical/_special_linear.py
@@ -195,7 +195,8 @@ def get_ref_stepper(self) -> ex.BaseStepper:
             domain_extent=self.domain_extent,
             num_points=self.num_points,
             dt=self.dt,
-            hyper_diffusivity=self.hyp_diffusion_coef,
+            # Need minus to match the sign of the diffusion coefficient
+            hyper_diffusivity=-self.hyp_diffusion_coef,
             diffuse_on_diffuse=True,
         )
 
@@ -205,7 +206,8 @@ def get_coarse_stepper(self) -> ex.BaseStepper:
             domain_extent=self.domain_extent,
             num_points=self.num_points,
             dt=self.dt * self.coarse_proportion,
-            hyper_diffusivity=self.hyp_diffusion_coef,
+            # Need minus to match the sign of the diffusion coefficient
+            hyper_diffusivity=-self.hyp_diffusion_coef,
             diffuse_on_diffuse=True,
         )