[Users] Multiple time levels
Steven R. Brandt
sbrandt at cct.lsu.edu
Tue Jan 5 15:11:14 CST 2016
OK, so copying data in CCTK_PRESTEP seems to work.
My next question is, I think, trickier... My par file looks something
Carpet::time_refinement_factors = "[1,1]"
Time::timestep_method = "given"
Time::timestep = 0.05
Carpet::prolongation_order_time = 0
Unfortunately, when I try to run the code I get this:
WARNING level 0 from host 18-5e-f-12-9b-2a.wlan.lsu.edu process 0
while executing schedule bin (none), routine (no thorn)::(no routine)
in thorn CarpetLib, file
-> Internal error: extrapolation in time. variable=FUNWAVE::eta
int Carpet::Abort(const cGH*, int): Assertion `0' failed.
It's trying to get data from the same time level, but it seems to only
have the wrong one.
My theory is that prolongation_order_time = 1 is going to give me what I
want, because it's going to prolong to one of the endpoints (the code
runs, at least). Does this mean that prolongation_order_time = 0 is broken?
On 01/05/2016 01:06 PM, Erik Schnetter wrote:
> If a grid variable has multiple time levels, then it will be cycled at
> the beginning of every time step. You need to copy the past to the
> current time level after time level cycling to preserve the previous
> On Tue, Jan 5, 2016 at 1:46 PM, Steven R. Brandt <sbrandt at cct.lsu.edu> wrote:
>> In the code I'm working on, I have a set of variables living at a single
>> timestep. For various reasons, I want to shift to having them at two
>> timelevels, but simply making the change in interface and schedule
>> results in failure to run old test suites. Is there a recipe I should be
>> following to make this work?
>> Users mailing list
>> Users at einsteintoolkit.org
More information about the Users