<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>Since Baikal does not support Llama (it does not, does it?) all you should have to do is use <code>GetBoundarySpecification</code> instead of <code>GetBoundarySizesAndTypes</code>. Lean_public which does, or at least there were plans to do so, support Llama has to do more work.</p>
<p>Ie this should work (and does for me using repos/wvuthorns/BaikalVacuum/test/BaikalVacuum_EE_O8_sgw3d.par even when increasing resolution to dx=0.0.1 so that not all MPI ranks have a boundary face all the time) :</p>
<div class="codehilite language-diff"><pre><span></span><span class="gh">diff --git a/BaikalVacuum/src/BoundaryConditions.c b/BaikalVacuum/src/BoundaryConditions.c</span>
<span class="gh">index 7240e5b..f89a9d2 100644</span>
<span class="gd">--- a/BaikalVacuum/src/BoundaryConditions.c</span>
<span class="gi">+++ b/BaikalVacuum/src/BoundaryConditions.c</span>
<span class="gu">@@ -96,11 +96,11 @@ void BaikalVacuum_BoundaryConditions_aux_gfs(CCTK_ARGUMENTS) {</span>
   CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0;

   CCTK_INT bndsize[6];
<span class="gd">-  CCTK_INT is_ghostbnd[6];</span>
<span class="gd">-  CCTK_INT is_symbnd[6];</span>
<span class="gd">-  CCTK_INT is_physbnd[6];</span>
<span class="gd">-  GetBoundarySizesAndTypes(cctkGH, 6, bndsize, is_ghostbnd, is_symbnd, is_physbnd);</span>
<span class="gi">+  CCTK_INT is_internal[6];</span>
<span class="gi">+  CCTK_INT is_staggered[6];</span>
<span class="gi">+  CCTK_INT shiftout[6];</span>

<span class="gi">+  GetBoundarySpecification(6, bndsize, is_internal, is_staggered, shiftout);</span>

   ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, bndsize[0], -1, &quot;BaikalVacuum::HGF&quot;, &quot;flat&quot;);
   if (ierr &lt; 0) CCTK_ERROR(&quot;Failed to register BC for BaikalVacuum::HGF!&quot;);
</pre></div>


<p>‌</p>
<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>