disable conservative_fixing again

src/SRHD/rhs.jl

@@ -280,6 +280,7 @@ function rhs_mono!(cache, mesh, equation, P)
   # fd_diff_1d!(ldg_tau, tau, x, mesh)
 
 
+  # dg1d.new_broadcast_volume!(conservative_fixing, equation, mesh)
   dg1d.new_broadcast_volume!(cons2prim_kastaun, equation, mesh)
   dg1d.new_broadcast_volume!(maxspeed, equation, mesh)
 
@@ -332,9 +333,11 @@ function rhs_covariant!(cache, mesh, equation, P)
           bdry_ldg_D, bdry_ldg_S, bdry_ldg_tau,
           bdry_rhs_D, bdry_rhs_S, bdry_rhs_tau,
           bdry_smoothed_mu      = get_bdry_variables(cache)
+  @unpack x = get_static_variables(cache)
 
   # compute rhs for ∂t u + ∂x f(x) = 0
 
+  # dg1d.new_broadcast_volume!(conservative_fixing, equation, mesh)
   dg1d.new_broadcast_volume!(cons2prim_kastaun, equation, mesh)
   dg1d.new_broadcast_volume!(maxspeed, equation, mesh)
 
@@ -365,6 +368,9 @@ function rhs_covariant!(cache, mesh, equation, P)
   compute_rhs_weak_form!(ldg_D,   D,   nflx_D,   mesh)
   compute_rhs_weak_form!(ldg_S,   S,   nflx_S,   mesh)
   compute_rhs_weak_form!(ldg_tau, tau, nflx_tau, mesh)
+  # fd_diff_1d!(ldg_D,   D,   x, mesh)
+  # fd_diff_1d!(ldg_S,   S,   x, mesh)
+  # fd_diff_1d!(ldg_tau, tau, x, mesh)
 
   # need to recompute, because ldg computation overwrites nflx_{D,S,tau}
   # numerical fluxes