diff --git a/acc/mhd_solver/stencil_process.sps b/acc/mhd_solver/stencil_process.sps index 213d879..2f9ef2b 100644 --- a/acc/mhd_solver/stencil_process.sps +++ b/acc/mhd_solver/stencil_process.sps @@ -39,9 +39,9 @@ value(in Vector uu) Scalar upwd_der6(in Vector uu, in Scalar lnrho) { - Scalar uux = value(uu).x; - Scalar uuy = value(uu).y; - Scalar uuz = value(uu).z; + Scalar uux = fabs(value(uu).x); + Scalar uuy = fabs(value(uu).y); + Scalar uuz = fabs(value(uu).z); return (Scalar){uux*der6x_upwd(lnrho) + uuy*der6y_upwd(lnrho) + uuz*der6z_upwd(lnrho)}; } #endif