WIP: added missing factor of 0.5 in surface flux

paranumal · Sep 13, 2023 · 99cbb2a · 99cbb2a
1 parent 2c97279
commit 99cbb2a
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/solvers/maxwell/okl/maxwellSurfaceTet3D.okl b/solvers/maxwell/okl/maxwellSurfaceTet3D.okl
@@ -176,7 +176,7 @@ void upwind(const dfloat nx,
             }
 
             // evaluate "flux" terms: (sJ/J)*(A*nx+B*ny)*(q^* - q^-)
-            const dfloat sc = invJ*sJ;
+            const dfloat sc = p_half*invJ*sJ;
 
             dfloat fluxHx, fluxHy, fluxHz;
 	    dfloat fluxEx, fluxEy, fluxEz;