diff --git a/src_mpi/claw.cc b/src_mpi/claw.cc index 726c8c6..c1c6e95 100644 --- a/src_mpi/claw.cc +++ b/src_mpi/claw.cc @@ -740,9 +740,12 @@ void ConservationLaw::iterate_explicit (IntegratorExplicit& integrator // Forward euler step in case of explicit scheme // In case of implicit scheme, this is the update current_solution += newton_update; + current_solution.compress (VectorOperation::insert); // current_solution = ark*old_solution + (1-ark)*current_solution current_solution.sadd (1.0-ark[rk], ark[rk], old_solution); + current_solution.compress (VectorOperation::insert); + compute_cell_average (); compute_shock_indicator ();