<html>#2639: CarpetLib internal error with PreSync and analysis thorns
<table style='border-spacing: 1ex 0pt; '>
<tr><td style='text-align:right'> Reporter:</td><td>Samuel Cupp</td></tr>
<tr><td style='text-align:right'> Status:</td><td>new</td></tr>
<tr><td style='text-align:right'>Milestone:</td><td></td></tr>
<tr><td style='text-align:right'> Version:</td><td>development version</td></tr>
<tr><td style='text-align:right'> Type:</td><td>bug</td></tr>
<tr><td style='text-align:right'> Priority:</td><td>major</td></tr>
<tr><td style='text-align:right'>Component:</td><td>Carpet</td></tr>
</table>
<p>Comment (by Samuel Cupp):</p>
<p>So, the heart of the problem comes from (I think) the function</p>
<div class="codehilite"><pre><span></span><span class="c">th::regrid()</span>
</pre></div>
<p>in repos/carpet/CarpetLib/src/th.cc:51. When this runs, it sets the</p>
<div class="codehilite"><pre><span></span><span>deltas.AT(ml).AT(rl)</span>
</pre></div>
<p>object to 1.0 (possibly divided by the timereffacts), which then causes an incorrect cctk_time. This seems to be polluting some functions which call get_time. This happens without PreSync, but by miracle (design?) it gets reset to the right value before anything tries to use it. This does not occur in time with PreSync.</p>
<p>--<br/>
Ticket URL: <a href='https://bitbucket.org/einsteintoolkit/tickets/issues/2639/carpetlib-internal-error-with-presync-and'>https://bitbucket.org/einsteintoolkit/tickets/issues/2639/carpetlib-internal-error-with-presync-and</a></p>
</html>