<div dir="ltr">Miguel<div><br></div><div>Yes, Carpet prolongates with the same order in each direction. However, in your case, prolongation in the z direction is exact and does not require any ghost zones, since nothing needs to happen in the z direction. Carpet should handle this correctly. Can you post a complete, simplified parameter file to reproduce the problem?</div><div><br></div><div>-erik</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 14, 2018 at 10:36 PM, 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 all,<br>
<br>
many thanks for the parameter file, Steve. it works, but as soon as i set<br>
Carpet::prolongation_order_<wbr>space = 5, i get inconsistent grid structure type of errors:<br>
<br>
WARNING level 1 from host meurglysIII process 1<br>
<span class=""> while executing schedule bin (none), routine (no thorn)::(no routine)<br>
</span> in thorn CarpetLib, file /home/mzilhao/dev/ET/Cactus/<wbr>arrangements/Carpet/CarpetLib/<wbr>src/dh.cc:161:<br>
-><br>
/home/mzilhao/dev/ET/Cactus/<wbr>arrangements/Carpet/CarpetLib/<wbr>src/dh.cc:813:<br>
[ml=0 rl=1 c=1] The following grid structure consistency check failed:<br>
Refinement prolongation: Send region must be contained in exterior<br>
send <= obox.exterior<br>
<br>
WARNING level 0 from host meurglysIII process 0<br>
<span class=""> while executing schedule bin (none), routine (no thorn)::(no routine)<br>
</span> in thorn CarpetLib, file /home/mzilhao/dev/ET/Cactus/<wbr>arrangements/Carpet/CarpetLib/<wbr>src/dh.cc:2118:<br>
-> The grid structure is inconsistent. It is impossible to continue.<br>
cactus_Lean_ET_master: /home/mzilhao/dev/ET/Cactus/<wbr>arrangements/Carpet/Carpet/<wbr>src/helpers.cc:275:<br>
int Carpet::Abort(const cGH*, int): Assertion `0' failed.<br>
<br>
so i guess it's as Ian said, that Carpet tries to prolongate with the same order in every direction?<br>
i find that in order for this to work i either have to decrease the prolongation to order 1, or put<br>
the boundary along the z direction back to a non-zero value. but i guess in this latter case i no<br>
longer have a truly 2D configuration, right?<br>
<br>
thanks,<br>
Miguel<br>
<div class="HOEnZb"><div class="h5"><br>
On 15/02/18 03:00, Steven R. Brandt wrote:<br>
> I have a parameter file which I have used fairly recently that does this. A very short running test<br>
> case is attached.<br>
><br>
> --Steve<br>
><br>
><br>
> On 02/14/2018 09:21 AM, Miguel Zilhão wrote:<br>
>> hi Ian,<br>
>><br>
>> many thanks for the parameter file. this allowed me to go a little further (i was missing<br>
>> CoordBase::boundary_shiftout_<wbr>z_lower = 1 and CoordBase::boundary_shiftout_<wbr>z_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/<wbr>Carpet/src/SetupGH.cc:2512:<br>
>> -> There are not enough ghost zones for the desired spatial prolongation order on map 0,<br>
>> 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_<wbr>space = 5<br>
>> Carpet::prolongation_order_<wbr>time = 2<br>
>><br>
>> driver::ghost_size_x <wbr> = 3<br>
>> driver::ghost_size_y <wbr> = 3<br>
>> driver::ghost_size_z <wbr> = 0<br>
>><br>
>> i'm guessing that in the example you provided things worked because you only had one grid? is there<br>
>> 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">ian.hinder@aei.mpg.de</a> wrote:<br>
>>><br>
>>>> On 13 Feb 2018, at 06:19, Miguel Zilhão <<a href="mailto:miguel.zilhao.nogueira@tecnico.ulisboa.pt">miguel.zilhao.nogueira@<wbr>tecnico.ulisboa.pt</a><br>
>>>> <mailto:<a href="mailto:miguel.zilhao.nogueira@tecnico.ulisboa.pt">miguel.zilhao.<wbr>nogueira@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>
>>>> cactus_Lean_ET: Cactus/arrangements/Carpet/<wbr>CarpetLib/src/gh.cc:61 <<a href="http://gh.cc:61/" rel="noreferrer" target="_blank">http://gh.cc:61/</a>>: gh::gh(const<br>
>>>> std::vector<vect<int, 3> >&, centering, int, centering, const std::vector<std::vector<bbox<<wbr>int,<br>
>>>> 3> ><br>
>>>>> &, const i2vect&): Assertion `all(box.shape() / box.stride() >= boundary_width[0] +<br>
>>>> boundary_width[1])' failed.<br>
>>>><br>
>>>> here's how i'm specifying my grid:<br>
>>>><br>
>>>> CoordBase::xmin <wbr> = -48.00<br>
>>>> CoordBase::ymin <wbr> = -48.00<br>
>>>> CoordBase::zmin <wbr> = 0.00<br>
>>>> CoordBase::xmax <wbr> = +48.00<br>
>>>> CoordBase::ymax <wbr> = +48.00<br>
>>>> CoordBase::zmax <wbr> = +0.00<br>
>>>> CoordBase::dx <wbr> = 2.00<br>
>>>> CoordBase::dy <wbr> = 2.00<br>
>>>> CoordBase::dz <wbr> = 2.00<br>
>>>><br>
>>>> driver::ghost_size_x <wbr> = 3<br>
>>>> driver::ghost_size_y <wbr> = 3<br>
>>>> driver::ghost_size_z <wbr> = 0<br>
>>>><br>
>>>> CoordBase::boundary_size_x_<wbr>lower = 3<br>
>>>> CoordBase::boundary_size_y_<wbr>lower = 3<br>
>>>> CoordBase::boundary_size_z_<wbr>lower = 0<br>
>>>> CoordBase::boundary_size_x_<wbr>upper = 0<br>
>>>> CoordBase::boundary_size_y_<wbr>upper = 0<br>
>>>> CoordBase::boundary_size_z_<wbr>upper = 0<br>
>>>><br>
>>>> CoordBase::boundary_shiftout_<wbr>x_lower = 0<br>
>>>> CoordBase::boundary_shiftout_<wbr>y_lower = 0<br>
>>>> CoordBase::boundary_shiftout_<wbr>z_lower = 0<br>
>>>><br>
>>>> am i missing something?<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<br>
>>> 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/<wbr>Kranc/blob/master/Examples/<wbr>laplace.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/<wbr>ianhin</a><br>
>>><br>
>> ______________________________<wbr>_________________<br>
>> Users mailing list<br>
>> <a href="mailto:Users@einsteintoolkit.org">Users@einsteintoolkit.org</a><br>
>> <a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.<wbr>org/mailman/listinfo/users</a><br>
><br>
><br>
><br>
> ______________________________<wbr>_________________<br>
> Users mailing list<br>
> <a href="mailto:Users@einsteintoolkit.org">Users@einsteintoolkit.org</a><br>
> <a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.<wbr>org/mailman/listinfo/users</a><br>
><br>
______________________________<wbr>_________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.<wbr>org/mailman/listinfo/users</a><br>
</div></div></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>