[Users] 2D configuration with Carpet

ian.hinder at aei.mpg.de ian.hinder at aei.mpg.de
Wed Feb 14 07:25:52 CST 2018



> On 13 Feb 2018, at 06:19, Miguel Zilhão <miguel.zilhao.nogueira at 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[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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20180214/0dba8ebc/attachment-0001.html 


More information about the Users mailing list