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/Carpet/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, [email protected] wrote: > > >> On 13 Feb 2018, at 06:19, Miguel Zilhão >> <[email protected] >> <mailto:[email protected]>> 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[0] + >> boundary_width[1])' 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 > _______________________________________________ Users mailing list [email protected] http://lists.einsteintoolkit.org/mailman/listinfo/users
