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: 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.00CoordBase::ymin = -48.00CoordBase::zmin = 0.00CoordBase::xmax = +48.00CoordBase::ymax = +48.00CoordBase::zmax = +0.00CoordBase::dx = 2.00CoordBase::dy = 2.00CoordBase::dz = 2.00driver::ghost_size_x = 3driver::ghost_size_y = 3driver::ghost_size_z = 0CoordBase::boundary_size_x_lower = 3CoordBase::boundary_size_y_lower = 3CoordBase::boundary_size_z_lower = 0CoordBase::boundary_size_x_upper = 0CoordBase::boundary_size_y_upper = 0CoordBase::boundary_size_z_upper = 0CoordBase::boundary_shiftout_x_lower = 0CoordBase::boundary_shiftout_y_lower = 0CoordBase::boundary_shiftout_z_lower = 0am 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):
I have not tried this recently, but it worked at one point.