<div dir="ltr">Miguel<div><br></div><div>Apparently you do need to keep requesting 3 ghost zones. This won&#39;t extend the grid since you still have 0 boundary zones.</div><div><br></div><div>-erik</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 14, 2018 at 10:21 AM, Miguel Zilhão <span dir="ltr">&lt;<a href="mailto:miguel.zilhao.nogueira@tecnico.ulisboa.pt" target="_blank">miguel.zilhao.nogueira@tecnico.ulisboa.pt</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">hi Ian,<br>
<br>
many thanks for the parameter file. this allowed me to go a little further (i was missing CoordBase::boundary_shiftout_z<wbr>_lower = 1 and CoordBase::boundary_shiftout_z<wbr>_upper = 1).<br>
however, i&#39;m still getting the following error:<br>
<br>
ERROR from host meurglysIII process 0<br>
  while executing schedule bin (none), routine (no thorn)::(no routine)<br>
  in thorn Carpet, file ./Cactus/arrangements/Carpet/C<wbr>arpet/src/SetupGH.cc:2512:<br>
  -&gt; There are not enough ghost zones for the desired spatial prolongation order on map 0, refinement level 0.  With a spatial prolongation order of 5, you need at least 3 ghost zones.<br>
<br>
i&#39;m setting:<br>
<br>
Carpet::prolongation_order_spa<wbr>ce        = 5<br>
Carpet::prolongation_order_tim<wbr>e         = 2<span class=""><br>
<br>
driver::ghost_size_x                 = 3<br>
driver::ghost_size_y                 = 3<br>
driver::ghost_size_z                 = 0<br>
<br></span>
i&#39;m guessing that in the example you provided things worked because you only had one grid? is there any way of doing this with more inner levels?<br>
<br>
thanks,<br>
Miguel<br>
<br>
On 14/02/18 22:25, <a href="mailto:ian.hinder@aei.mpg.de" target="_blank">ian.hinder@aei.mpg.de</a> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">
On 13 Feb 2018, at 06:19, Miguel Zilhão &lt;<a href="mailto:miguel.zilhao.nogueira@tecnico.ulisboa.pt" target="_blank">miguel.zilhao.nogueira@tecnic<wbr>o.ulisboa.pt</a> &lt;mailto:<a href="mailto:miguel.zilhao.nogueira@tecnico.ulisboa.pt" target="_blank">miguel.zilhao.nogueira<wbr>@tecnico.ulisboa.pt</a>&gt;&gt; wrote:<br>
<br>
hi Erik,<br>
<br>
thanks for your reply. i&#39;ve tried to set the parameters as you described, but i&#39;m getting a Carpet<br>
assertion failure:<br>
<br></span>
cactus_Lean_ET: Cactus/arrangements/Carpet/Car<wbr>petLib/src/gh.cc:61 &lt;<a href="http://gh.cc:61/" rel="noreferrer" target="_blank">http://gh.cc:61/</a>&gt;: gh::gh(const<div><div class="h5"><br>
std::vector&lt;vect&lt;int, 3&gt; &gt;&amp;, centering, int, centering, const std::vector&lt;std::vector&lt;bbox&lt;i<wbr>nt, 3&gt; &gt;<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
&amp;, const i2vect&amp;): Assertion `all(box.shape() / box.stride() &gt;= boundary_width[0] +<br>
</blockquote>
boundary_width[1])&#39; failed.<br>
<br>
here&#39;s how i&#39;m specifying my grid:<br>
<br>
CoordBase::xmin                      = -48.00<br>
CoordBase::ymin                      = -48.00<br>
CoordBase::zmin                      =   0.00<br>
CoordBase::xmax                      = +48.00<br>
CoordBase::ymax                      = +48.00<br>
CoordBase::zmax                      =  +0.00<br>
CoordBase::dx                        =   2.00<br>
CoordBase::dy                        =   2.00<br>
CoordBase::dz                        =   2.00<br>
<br>
driver::ghost_size_x                 = 3<br>
driver::ghost_size_y                 = 3<br>
driver::ghost_size_z                 = 0<br>
<br>
CoordBase::boundary_size_x_low<wbr>er     = 3<br>
CoordBase::boundary_size_y_low<wbr>er     = 3<br>
CoordBase::boundary_size_z_low<wbr>er     = 0<br>
CoordBase::boundary_size_x_upp<wbr>er     = 0<br>
CoordBase::boundary_size_y_upp<wbr>er     = 0<br>
CoordBase::boundary_size_z_upp<wbr>er     = 0<br>
<br>
CoordBase::boundary_shiftout_x<wbr>_lower = 0<br>
CoordBase::boundary_shiftout_y<wbr>_lower = 0<br>
CoordBase::boundary_shiftout_z<wbr>_lower = 0<br>
<br>
am i missing something?<br>
</div></div></blockquote><div><div class="h5">
<br>
Hi Miguel,<br>
<br>
There is an example parameter file in Kranc, for a true 2D Laplace equation (i.e. it doesn&#39;t take derivatives in the z direction):<br>
<br>
<a href="https://github.com/ianhinder/Kranc/blob/master/Examples/laplace.par" rel="noreferrer" target="_blank">https://github.com/ianhinder/K<wbr>ranc/blob/master/Examples/lapl<wbr>ace.par</a><br>
<br>
I have not tried this recently, but it worked at one point.<br>
<br>
-- <br>
Ian Hinder<br>
<a href="http://members.aei.mpg.de/ianhin" rel="noreferrer" target="_blank">http://members.aei.mpg.de/ianh<wbr>in</a><br>
<br>
</div></div></blockquote>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>Erik Schnetter &lt;<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>&gt;<br><a href="http://www.perimeterinstitute.ca/personal/eschnetter/" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a></div><div><br></div></div></div>
</div>