<div dir="ltr">Ian<div><br></div><div>The operating system will zero all pages before the application sees them. Malloc and friends will reuse memory from the operating system without zeroing it. In practice, the first malloc calls will see all zeros, later ones (especially for smaller regions) might see non-zeros. Carpet does not set anything to zero, but it can set things to nan if you ask it to. MoL will also set the RHS to zero by default before every iteration.</div><div><br></div><div>-erik</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Oct 7, 2016 at 9:31 AM, Roland Haas <span dir="ltr"><<a href="mailto:rhaas@illinois.edu" target="_blank">rhaas@illinois.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello Ian,<br>
<br>
ah, these are all tricky questions.<br>
<span class=""><br>
> which is a bit surprising. This might come from the old issue of the<br>
> intersection between outer boundaries and interpatch boundaries not<br>
> being correctly initialised by the combination of Interpolate2 and<br>
> McLachlan. Last time I checked though, this may have been fixed. If<br>
> you don't use poison, does Carpet zero the variables, or leave them<br>
> uninitialised? In which case, does the OS zero them? Presumably new<br>
> memory from the OS is zeroed (for security/privacy reasons), but if<br>
> the memory has been malloc/freed already by the cactus process, then<br>
> maybe it is random.<br>
</span>Yes, the OS would zero, malloc does not have to. new() on the other<br>
hand will also zero (you'd have to check if CarpetLib uses new() or<br>
malloc() or posix_memalign() to get new memory, I believe it is<br>
actually malloc).<br>
<span class=""><br>
> Now, another related question. Is it possible to change the number<br>
> of CarpetRegrid2 centres during evolution? The num_centres parameter<br>
> is non-steerable, and indeed CarpetRegrid2 only looks at the centre<br>
> parameters initially (otherwise every recovery would overwrite the<br>
> dynamical centre variables with the parameter values). Ideally, I<br>
> would have set num_centres to a larger value when I started the run,<br>
> but I didn't.<br>
><br>
> We just need to make sure that CarpetRegrid2 runs INIT_CENTRE (from<br>
> initialise.cc) on any newly added centres on recovery. This is done<br>
> in WRAGH, and then similar to the previous issue, overwritten by<br>
> recovery. So this should just work?<br>
</span>Yes that would work. You just have to make sure you can identify new<br>
levels, eg. by their number of levels being -42 or similar. The<br>
alternative is to use the Trigger thorn and set the grid scalars (see<br>
the NSNS gallery parfile and look for<br>
<br>
Trigger::Trigger_Steered_<wbr>Scalar[1] = "CarpetRegrid2::num_levels[2]"<br>
<br>
Since you can change the grid scalars and "activate" more levels.<br>
<div class="HOEnZb"><div class="h5"><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://keys.gnupg.net" rel="noreferrer" target="_blank">http://keys.gnupg.net</a>.<br>
</div></div><br>______________________________<wbr>_________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.<wbr>org/mailman/listinfo/users</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>><br><a href="http://www.perimeterinstitute.ca/personal/eschnetter/" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a></div>
</div>