fix scaling issue for 1D advection (#46)

pdebench · Aug 24, 2023 · 524848f · 524848f
1 parent f9fd131
commit 524848f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/pdebench/data_gen/data_gen_NLE/AdvectionEq/advection_multi_solution_Hydra.py b/pdebench/data_gen/data_gen_NLE/AdvectionEq/advection_multi_solution_Hydra.py
@@ -261,7 +261,7 @@ def flux(u):
         fL = uL * beta
         fR = uR * beta
         # upwind advection scheme
-        f_upwd = (fR[1:cfg.multi.nx+2] + fL[2:cfg.multi.nx+3]
+        f_upwd = 0.5 * (fR[1:cfg.multi.nx+2] + fL[2:cfg.multi.nx+3]
                - jnp.abs(beta)*(uL[2:cfg.multi.nx+3] - uR[1:cfg.multi.nx+2]))
         return f_upwd