[ET Trac] [Einstein Toolkit] #1797: allow new timelevels to be created at runtime
Einstein Toolkit
trac-noreply at einsteintoolkit.org
Fri Jul 24 10:07:29 CDT 2015
#1797: allow new timelevels to be created at runtime
--------------------------+-------------------------------------------------
Reporter: rhaas | Owner:
Type: enhancement | Status: review
Priority: optional | Milestone:
Component: Cactus | Version: development version
Resolution: | Keywords: Cactus Carpet PUGH
--------------------------+-------------------------------------------------
Comment (by rhaas):
Replying to [comment:5 hinder]:
> So timelevel cycling would then have the effect of copying tl_N into
tl_{N+1}, and leaving tl0 as it is, so that it is identical to tl1? This
sounds like a good idea to me. If we can then omit MoL's copy, I think we
can then use MoL with only a single timelevel with, e.g. RK4, which is all
it should need.
Avoiding the copying does not avoid multiple timelevels since MoL uses
tl=0 as scratch space during the computation. Eg the final
LinearCombination for RK4 is:
y_0_new = 1.0 * y_0_from_last_RK_step + scratchVar - 4/3 * y_1
What Erik proposes does not change anything from the point of the MoL time
intergrator, it just moves MoL's Initialcopy from the beginning of EVOL
(when MoL's group runs) to just before CCTK_PRESTEP (when Carpet cycles
timelevels).
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/1797#comment:7>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
More information about the Trac
mailing list