<div dir="ltr">Miguel<div><br></div><div>Apparently you do need to keep requesting 3 ghost zones. This won'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"><<a href="mailto:miguel.zilhao.nogueira@tecnico.ulisboa.pt" target="_blank">miguel.zilhao.nogueira@tecnico.ulisboa.pt</a>></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'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>
-> 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'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'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 <<a href="mailto:miguel.zilhao.nogueira@tecnico.ulisboa.pt" target="_blank">miguel.zilhao.nogueira@tecnic<wbr>o.ulisboa.pt</a> <mailto:<a href="mailto:miguel.zilhao.nogueira@tecnico.ulisboa.pt" target="_blank">miguel.zilhao.nogueira<wbr>@tecnico.ulisboa.pt</a>>> wrote:<br>
<br>
hi Erik,<br>
<br>
thanks for your reply. i've tried to set the parameters as you described, but i'm getting a Carpet<br>
assertion failure:<br>
<br></span>
cactus_Lean_ET: Cactus/arrangements/Carpet/Car<wbr>petLib/src/gh.cc:61 <<a href="http://gh.cc:61/" rel="noreferrer" target="_blank">http://gh.cc:61/</a>>: gh::gh(const<div><div class="h5"><br>
std::vector<vect<int, 3> >&, centering, int, centering, const std::vector<std::vector<bbox<i<wbr>nt, 3> ><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
&, const i2vect&): Assertion `all(box.shape() / box.stride() >= boundary_width[0] +<br>
</blockquote>
boundary_width[1])' failed.<br>
<br>
here's how i'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'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 <<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>><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>