Skip to content

Commit

Permalink
fix the new source fill patch
Browse files Browse the repository at this point in the history
  • Loading branch information
zingale committed Aug 31, 2019
1 parent 8a8f89c commit 719d5e9
Showing 1 changed file with 2 additions and 5 deletions.
7 changes: 2 additions & 5 deletions Source/driver/Castro_advance_sdc.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -251,11 +251,8 @@ Castro::do_advance_sdc (Real time,
// store the new solution
MultiFab::Copy(S_new, *(k_new[SDC_NODES-1]), 0, 0, S_new.nComp(), 0);


// I think this bit only needs to be done for the last iteration...

// We need to make source_old and source_new be the source terms at
// the old and new time. we never actually evaluate the sources
// the old and new time. we never actually evaluated the sources
// using the new time state (since we just constructed it). Note:
// we always use do_old_sources here, since we want the actual
// source and not a correction.
Expand All @@ -273,7 +270,7 @@ Castro::do_advance_sdc (Real time,
clean_state(S_new, cur_time, 0);
expand_state(Sborder, cur_time, Sborder.nGrow());
do_old_sources(new_source, Sborder, cur_time, dt, amr_iteration, amr_ncycle);
AmrLevel::FillPatch(*this, old_source, old_source.nGrow(), cur_time, Source_Type, 0, NUM_STATE);
AmrLevel::FillPatch(*this, new_source, new_source.nGrow(), cur_time, Source_Type, 0, NUM_STATE);
}

finalize_do_advance(time, dt, amr_iteration, amr_ncycle);
Expand Down

0 comments on commit 719d5e9

Please sign in to comment.