<div dir="ltr"><div dir="ltr">Thanks a lot Zach, Roland and Ian for your comments. We are currently compiling the last version of the EinsteinToolkit and we will make some runs with different values of CarpetRegrid2::min_fraction to check if the noise disappears. We will report you the results in a few days.</div><div dir="ltr"><br></div><div>Best Regards,</div><div>Toni.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El jue., 7 mar. 2019 a las 14:47, Ian Hinder (<<a href="mailto:ian.hinder@manchester.ac.uk">ian.hinder@manchester.ac.uk</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">
<div style="overflow-wrap: break-word;">
<br>
<div><br>
<blockquote type="cite">
<div>On 7 Mar 2019, at 07:40, Antoni Ramos Buades <<a href="mailto:antoniramosbuades@gmail.com" target="_blank">antoniramosbuades@gmail.com</a>> wrote:</div>
<br class="gmail-m_1493953409875899807Apple-interchange-newline">
<div>
<div dir="ltr">Hi,
<div><br>
</div>
<div>Roland you are right it is forming and destroying a common box each time the outburst of noise happens. Do you know how one could get rid of this noise, or if there are some parameters of Carpet which one could add to the parameter file to avoid
such a behaviour?</div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>Hi Antoni,</div>
<div><br>
</div>
<div>There is this parameter in CarpetRegrid2:</div>
<div><br>
</div>
<div>
<div style="background-color:rgb(255,255,254);font-family:SFMono-Medium,"SF Mono","Segoe UI Mono","Roboto Mono","Ubuntu Mono",Menlo,monospace;font-size:13px;line-height:20px;white-space:pre-wrap">
<div><span style="color:rgb(32,32,32)">CCTK_REAL</span> <span style="color:rgb(32,32,32)">
min_fraction</span> <span style="color:rgb(191,38,0)">"Minimum fraction of required refined points that need to be present in a refined region"</span>
<span style="color:rgb(32,32,32)">STEERABLE</span>=<span style="color:rgb(32,32,32)">always</span></div>
<div>{</div>
<div><span style="color:rgb(101,84,192)">0</span>:* :: <span style="color:rgb(191,38,0)">
""</span></div>
<div>} <span style="color:rgb(101,84,192)">0.9</span></div>
</div>
</div>
<div><br>
</div>
<div>When two boxes overlap, CarpetRegrid2 has to make a decision about whether to use the union of the two boxes, or to replace the two boxes with a single enclosing box. This code is in carpet/CarpetRegrid2/src/<a href="http://property.cc" target="_blank">property.cc</a> in
the function combine_regions::test_impl. The decision is made based on the number of point in the union of the two boxes, vs the number of points in a single enclosing box. It will leave the regions separate if </div>
<div><br>
</div>
</div>
<blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px">
<div>
<div>min_fraction * combined_size > regions_size</div>
<div><br>
</div>
</div>
</blockquote>
which, with the default of 0.9, means that the number of points in the single enclosing region is greater than 1.11 times the number of points in the original regions. The logic for this is that having lots of regions means lots of faces, edges and corners,
and more communication and prolongation, which can affect performance and might also be undesirable due to creating numerical error features such as reflections. On the other hand, the single enclosing region will contain more points, and hence will be more
expensive to evolve. min_fraction allows you to decide how many extra points you are willing to evolve, for the sake of having only a single box. If there would be more than 1/min_fraction times the number of points by combining, the code does not combine.
<div><br>
</div>
<div>If you set
<div>
<div>
<div><br>
</div>
</div>
</div>
<blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px">
<div>
<div>
<div>CarpetRegrid2::min_fraction = 1</div>
</div>
</div>
</blockquote>
<div>
<div>
<div><br>
</div>
<div>then Carpet will never create a single enclosing box, but will always give you a box based on the union of the points in the two boxes.</div>
<div><br>
</div>
<div>Ninja'd by Roland. Sigh...</div>
<div><br>
</div>
</div>
<div>
<div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
<div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
<div dir="auto" style="overflow-wrap: break-word;">
<div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
-- </div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
Ian Hinder</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
Research Software Engineer</div>
<div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
University of Manchester, UK</div>
</div>
</div>
</div>
</div>
<br>
</div>
</div>
</div>
</blockquote></div>