Miguel Apparently you do need to keep requesting 3 ghost zones. This won't extend the grid since you still have 0 boundary zones.
-erik On Wed, Feb 14, 2018 at 10:21 AM, Miguel Zilhão < miguel.zilhao.nogue...@tecnico.ulisboa.pt> wrote: > hi Ian, > > many thanks for the parameter file. this allowed me to go a little further > (i was missing CoordBase::boundary_shiftout_z_lower = 1 and > CoordBase::boundary_shiftout_z_upper = 1). > however, i'm still getting the following error: > > ERROR from host meurglysIII process 0 > while executing schedule bin (none), routine (no thorn)::(no routine) > in thorn Carpet, file ./Cactus/arrangements/Carpet/C > arpet/src/SetupGH.cc:2512: > -> 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. > > i'm setting: > > Carpet::prolongation_order_space = 5 > Carpet::prolongation_order_time = 2 > > driver::ghost_size_x = 3 > driver::ghost_size_y = 3 > driver::ghost_size_z = 0 > > 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? > > thanks, > Miguel > > On 14/02/18 22:25, ian.hin...@aei.mpg.de wrote: > >> >> >> On 13 Feb 2018, at 06:19, Miguel Zilhão <miguel.zilhao.nogueira@tecnic >>> o.ulisboa.pt <mailto:miguel.zilhao.nogue...@tecnico.ulisboa.pt>> wrote: >>> >>> hi Erik, >>> >>> thanks for your reply. i've tried to set the parameters as you >>> described, but i'm getting a Carpet >>> assertion failure: >>> >>> cactus_Lean_ET: Cactus/arrangements/Carpet/CarpetLib/src/gh.cc:61 < >>> http://gh.cc:61/>: gh::gh(const >>> >>> std::vector<vect<int, 3> >&, centering, int, centering, const >>> std::vector<std::vector<bbox<int, 3> > >>> >>>> &, const i2vect&): Assertion `all(box.shape() / box.stride() >= >>>> boundary_width + >>>> >>> boundary_width)' failed. >>> >>> here's how i'm specifying my grid: >>> >>> CoordBase::xmin = -48.00 >>> CoordBase::ymin = -48.00 >>> CoordBase::zmin = 0.00 >>> CoordBase::xmax = +48.00 >>> CoordBase::ymax = +48.00 >>> CoordBase::zmax = +0.00 >>> CoordBase::dx = 2.00 >>> CoordBase::dy = 2.00 >>> CoordBase::dz = 2.00 >>> >>> driver::ghost_size_x = 3 >>> driver::ghost_size_y = 3 >>> driver::ghost_size_z = 0 >>> >>> CoordBase::boundary_size_x_lower = 3 >>> CoordBase::boundary_size_y_lower = 3 >>> CoordBase::boundary_size_z_lower = 0 >>> CoordBase::boundary_size_x_upper = 0 >>> CoordBase::boundary_size_y_upper = 0 >>> CoordBase::boundary_size_z_upper = 0 >>> >>> CoordBase::boundary_shiftout_x_lower = 0 >>> CoordBase::boundary_shiftout_y_lower = 0 >>> CoordBase::boundary_shiftout_z_lower = 0 >>> >>> am i missing something? >>> >> >> Hi Miguel, >> >> 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): >> >> https://github.com/ianhinder/Kranc/blob/master/Examples/laplace.par >> >> I have not tried this recently, but it worked at one point. >> >> -- >> Ian Hinder >> http://members.aei.mpg.de/ianhin >> >> -- Erik Schnetter <schnet...@cct.lsu.edu> http://www.perimeterinstitute.ca/personal/eschnetter/
_______________________________________________ Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users