<html>#2288: Lean's boundary condition setting for BSSN constraints only fills in a single ghostzone
<table style='border-spacing: 1ex 0pt; '>
<tr><td style='text-align:right'> Reporter:</td><td>Zach Etienne</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></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>EinsteinToolkit thorn</td></tr>
</table>
<p>Comment (by Roland Haas):</p>
<p>The function to use is actually <code>GetBoundarySpecification</code> not <code>GetBoundarySizesAndTypes</code>. Easiest is most likely to copy the code from ML_ADMConstraint’s <code>Kranc.cc</code> and <code>ML_ADMConstraints_evaluate.cc</code> files, ie:</p>
<div class="codehilite language-c++"><pre><span></span><span class="cm">/*********************************************************************</span>
<span class="cm"> * GetBoundaryWidths</span>
<span class="cm"> *********************************************************************/</span>
<span class="kt">void</span> <span class="nf">GetBoundaryWidths</span><span class="p">(</span><span class="n">cGH</span> <span class="k">const</span> <span class="o">*</span> <span class="k">restrict</span> <span class="k">const</span> <span class="n">cctkGH</span><span class="p">,</span> <span class="n">CCTK_INT</span> <span class="n">nboundaryzones</span><span class="p">[</span><span class="mi">6</span><span class="p">])</span>
<span class="p">{</span>
<span class="n">CCTK_INT</span> <span class="n">is_internal</span><span class="p">[</span><span class="mi">6</span><span class="p">];</span>
<span class="n">CCTK_INT</span> <span class="n">is_staggered</span><span class="p">[</span><span class="mi">6</span><span class="p">];</span>
<span class="n">CCTK_INT</span> <span class="n">shiftout</span><span class="p">[</span><span class="mi">6</span><span class="p">];</span>
<span class="kt">int</span> <span class="n">ierr</span> <span class="o">=</span> <span class="o">-</span><span class="mi">1</span><span class="p">;</span>
<span class="k">if</span> <span class="p">(</span><span class="n">CCTK_IsFunctionAliased</span> <span class="p">(</span><span class="s">"MultiPatch_GetBoundarySpecification"</span><span class="p">))</span> <span class="p">{</span>
<span class="kt">int</span> <span class="k">const</span> <span class="n">map</span> <span class="o">=</span> <span class="n">MultiPatch_GetMap</span> <span class="p">(</span><span class="n">cctkGH</spa
<p>--<br/>
Ticket URL: <a href='https://bitbucket.org/einsteintoolkit/tickets/issues/2288/leans-boundary-condition-setting-for-bssn'>https://bitbucket.org/einsteintoolkit/tickets/issues/2288/leans-boundary-condition-setting-for-bssn</a></p>
</html>