Hello Luciano,
I'm following a simple parfile in LlamaToy that uses Thornburg04nc and one that Roland shared a few months ago but I'm having seemingly basic issues. In particular, I'm getting the following error:
In the code, it's mentioned that this could be a problem from CarpetInterp? but not sure what's going on, honestly. I attach the parfile and the full .out, .err.
Based on the error and the comment in the code around the assert, namely:
// This does not hold if the caller requests the same // interpolation to be done into different output arrays. // This may happen e.g. when CarpetInterp needs to // differentiate in time. This is arguably a performance bug // in CarpetInterp. (See whether this goes away now.) (It // should!)
something seems to request data at a time not commensurate with the coarsest timestep.
Since you have no mesh refinement that is somewhat strange.
Could you provide one of the backtrace.txt files? So that one can get some idea which part of the code called the interpolator?
Looking at your parfile (thank you for including it) You should be able to use:
Carpet::prolongation_order_time = 0
since you only have 1 refinement level and also remove
CarpetLib::restriction_order_space = 3
which is not used for unigrid (or vertex centered refinement where restriction is always an exact copy so has infinite order).
Neither of these two should really do anything to your run though. The only effect may be that setting
Carpet::prolongation_order_time = 0
may give you a more useful error message from a routine earlier in the call stack.
Yours, Roland