[Users] 2D configuration with Carpet

Miguel Zilhão miguel.zilhao.nogueira at tecnico.ulisboa.pt
Wed Feb 14 09:21:25 CST 2018


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, ian.hinder at aei.mpg.de wrote:
> 
> 
>> On 13 Feb 2018, at 06:19, Miguel Zilhão <miguel.zilhao.nogueira at tecnico.ulisboa.pt 
>> <mailto: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
> 


More information about the Users mailing list