<div dir="ltr">Hi Roland,<div><br></div><div>thanks for the quick response.</div><div><br></div><div>Great points. I wasn't really sure how to use the backtrace, but now following the addresses there I found that the problem was in the Outflow thorn. In fact, paying closer attention to the parfile you shared a while back, I see that you don't use the Outflow thorn. Is it possible to use it within Llama maybe tweaking the "coord_system" parameter or the interpolator? Or should I consider alternatives?</div><div><br></div><div>Thanks again.</div><div><br></div><div>Cheers.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El jue, 21 oct 2021 a las 13:03, Roland Haas (<<a href="mailto:rhaas@illinois.edu">rhaas@illinois.edu</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello Luciano,<br>
<br>
> I'm following a simple parfile in LlamaToy that uses Thornburg04nc and one<br>
> that Roland shared a few months ago but I'm having seemingly basic issues.<br>
> In particular, I'm getting the following error:<br>
<br>
> In the code, it's mentioned that this could be a problem from CarpetInterp?<br>
> but not sure what's going on, honestly. I attach the parfile and the full<br>
> .out, .err.<br>
<br>
Based on the error and the comment in the code around the assert,<br>
namely:<br>
<br>
// This does not hold if the caller requests the same<br>
// interpolation to be done into different output arrays.<br>
// This may happen e.g. when CarpetInterp needs to<br>
// differentiate in time. This is arguably a performance bug<br>
// in CarpetInterp. (See whether this goes away now.) (It<br>
// should!)<br>
<br>
something seems to request data at a time not commensurate with the<br>
coarsest timestep.<br>
<br>
Since you have no mesh refinement that is somewhat strange.<br>
<br>
Could you provide one of the backtrace.txt files? So that one can get<br>
some idea which part of the code called the interpolator?<br>
<br>
Looking at your parfile (thank you for including it) You should be able<br>
to use:<br>
<br>
Carpet::prolongation_order_time = 0<br>
<br>
since you only have 1 refinement level and also remove<br>
<br>
CarpetLib::restriction_order_space = 3<br>
<br>
which is not used for unigrid (or vertex centered refinement where<br>
restriction is always an exact copy so has infinite order).<br>
<br>
Neither of these two should really do anything to your run though. The<br>
only effect may be that setting <br>
<br>
Carpet::prolongation_order_time = 0<br>
<br>
may give you a more useful error message from a routine earlier in the<br>
call stack.<br>
<br>
Yours,<br>
Roland<br>
<br>
-- <br>
My email is as private as my paper mail. I therefore support encrypting<br>
and signing email messages. Get my PGP key from <a href="http://pgp.mit.edu" rel="noreferrer" target="_blank">http://pgp.mit.edu</a> .<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="trebuchet ms, sans-serif"><b>Luc Combi</b></font></div><div dir="ltr"><div style="font-family:sans-serif"><div style="color:rgb(117,117,117)">Instituto Argentino de Radioastronomía </div><div style="color:rgb(117,117,117)">CONICET, Buenos Aires, Argentina</div><div style="color:rgb(117,117,117)"><a href="http://sites.google.com/view/lucianocombi" target="_blank">web page.</a><br></div></div><div style="color:rgb(33,33,33);font-family:sans-serif"><br><br></div></div></div></div></div></div></div></div></div></div></div></div>