<html><head><meta http-equiv="Content-Type" content="text/html charset=iso-8859-1"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><br><div><div>On 26 Oct 2017, at 23:23, Miguel Zilhão &lt;<a href="mailto:miguel.zilhao.nogueira@tecnico.ulisboa.pt">miguel.zilhao.nogueira@tecnico.ulisboa.pt</a>&gt; wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">hi all,<br><br>for the purposes of trying a particular gauge condition, we'd like to have a grid function (let's <br>call it "shift_t0") that's set to be equal to the shift at t=0. we don't want this function to <br>evolve at all; we just want this grid function to be accessible at all times (with exactly the same <br>values).<br><br>is there something in particular that we need to watch out for, when doing this? we ask because <br>we've noticed that when doing things in a "naive"* way, sometimes NaNs start developing at some <br>point in the evolution (seemingly near refinement boundaries and when regriding). we noticed that <br>exactly the same thing happens for the "puncture_u" grid function in TwoPunctures: this function is <br>only used at t=0; however, if one activates storage for it and asks for its output, at some point <br>NaNs start developing in the output, even though the function is never used again. so, what would be <br>the standard way of defining a grid function that is constant in time?<br><br>thanks,<br>Miguel &amp; Helvi<br><br>* by "naive", we mean essentially the following<br><br>in schedule.ccl we have<br><br>STORAGE: shift_t0[1]<br><br>in interface.ccl we have<br><br>CCTK_REAL shift_t0 type=gf timelevels=1 tags='tensortypealias="U" tags='prolongation="none"'<br>{<br> &nbsp;&nbsp;betax_t0 betay_t0 betaz_t0<br>} "shift vector at t=0"<br><br>and we set its value (on the whole grid) at a function call at CCTK_INITIAL after ADMBase_PostInitial<br></blockquote><div><br></div><div>Does this happen only after regridding? &nbsp;By not prolongating the refinement boundaries, it's not clear how the "new" bits of the grid are supposed to be filled when the fine grids move. &nbsp;Could this be the problem? &nbsp;If you give the gridfunction 3 timelevels and remove "prolongation=none", do the NaNs go away?</div><div><br></div></div><div apple-content-edited="true">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div>--&nbsp;</div><div>Ian Hinder</div><div><a href="http://members.aei.mpg.de/ianhin">http://members.aei.mpg.de/ianhin</a></div></div></div></div></div>
</div>
<br></body></html>