[ET Trac] #2601: Add four more parameters for carpet

Liwei Ji trac-noreply at einsteintoolkit.org
Mon Jun 20 07:16:47 CDT 2022


#2601: Add four more parameters for carpet

 Reporter: Liwei Ji
   Status: new
Milestone: 
  Version: development version
     Type: enhancement
 Priority: major
Component: Carpet

Comment (by Liwei Ji):

Hi @{557058:59e031ba-9bb5-4298-a472-7b99d0ae6f22} , here is a test case

```
ActiveThorns = "
        Boundary
        Carpet
        CarpetIOASCII
        CarpetLib
        # CarpetRegrid2
        CartGrid3D
        CoordBase
        IOUtil
        InitBase
        MoL
        SymBase
"

$ghost_size = 1
$boundary_size = 1
$num_levels = 1

IOASCII::out1D_vars = "
        grid::x
"

IO::out_dir       = $parfile
IO::out_fileinfo  = "none"
IO::parfile_write = no

IOASCII::out1D_every = 1
IOASCII::out1D_d     = no
IOASCII::out1D_y     = no
IOASCII::out1D_z     = no

Carpet::max_refinement_levels    = 1
Carpet::use_buffer_zones         = yes
Carpet::prolongation_order_space =  1
Carpet::prolongation_order_time  =  1

Carpet::init_fill_timelevels = yes

CartGrid3D::type         = "CoordBase"
CartGrid3D::domain       = "full"
CartGrid3D::avoid_origin = no

CoordBase::xmin = -1.
CoordBase::xmax = +1 - 0.125

CoordBase::ymin = -0.125
CoordBase::ymax = +0.125
CoordBase::zmin = -0.125
CoordBase::zmax = +0.125
CoordBase::dx   =  0.125
CoordBase::dy   =  0.125
CoordBase::dz   =  0.125

driver::ghost_size  = $ghost_size
Cactus::cctk_itlast = 0

Carpet::domain_from_CoordBase = yes
#Carpet::split_interior_points_x = yes
Carpet::processor_topology = "manual"
Carpet::processor_topology_3d_x = 4
Carpet::processor_topology_3d_y = 1
Carpet::processor_topology_3d_z = 1

CoordBase::boundary_size_x_lower = $boundary_size
CoordBase::boundary_size_y_lower = $boundary_size
CoordBase::boundary_size_z_lower = $boundary_size
CoordBase::boundary_size_x_upper = $boundary_size
CoordBase::boundary_size_y_upper = $boundary_size
CoordBase::boundary_size_z_upper = $boundary_size
CoordBase::boundary_shiftout_x_lower  = 1
CoordBase::boundary_shiftout_x_upper  = 1
#CoordBase::boundary_shiftout_y_lower  = 1
#CoordBase::boundary_shiftout_y_upper  = 1
#CoordBase::boundary_shiftout_z_lower  = 1
#CoordBase::boundary_shiftout_z_upper  = 1

#CarpetRegrid2::regrid_every   = 0
#CarpetRegrid2::verbose        = yes

#CarpetRegrid2::num_centres  = 1
#CarpetRegrid2::num_levels_1 = $num_levels
#CarpetRegrid2::position_x_1 = 0.0
#CarpetRegrid2::radius_1[1]  = 0.2
```

I made the following changes:

1. remove `CarpetRegrid2`, it seems that the new pars are not consistent with it,
2. set `Carpet::max_refinement_levels = 1` the new pars only can work with no refinement case for now,
3. add `CoordBase::boundary_shiftout_x_lower  = 1` and `CoordBase::boundary_shiftout_x_upper  = 1`, for some reason only when I turn on these. The boundary ghost points are counted as global points. And these option are on in `SphericalNR`.

After these changes, when `Carpet::split_interior_points_x = yes`, I got

```
ipos = [0,0,0], nprocs_dir = [4,1,1], cub = [5,3,3], rub0 = [18,3,3]
ipos = [1,0,0], nprocs_dir = [4,1,1], cub = [9,3,3], rub0 = [18,3,3]
ipos = [2,0,0], nprocs_dir = [4,1,1], cub = [13,3,3], rub0 = [18,3,3]
ipos = [3,0,0], nprocs_dir = [4,1,1], cub = [18,3,3], rub0 = [18,3,3]
```

when `Carpet::split_interior_points_x = no`, I got

```
ipos = [0,0,0], nprocs_dir = [4,1,1], cub = [5,3,3], rub0 = [18,3,3]
ipos = [1,0,0], nprocs_dir = [4,1,1], cub = [10,3,3], rub0 = [18,3,3]
ipos = [2,0,0], nprocs_dir = [4,1,1], cub = [14,3,3], rub0 = [18,3,3]
ipos = [3,0,0], nprocs_dir = [4,1,1], cub = [18,3,3], rub0 = [18,3,3]
```

basally, the distribution of `(clb, cub)` change in `x` change from `[(0,5), (5,9), (9,13), (13,18)]` to `[(0,5), (5,10), (10,14), (14,18)]`

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2601/add-four-more-parameters-for-carpet
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/trac/attachments/20220620/df1f6cdf/attachment.html 


More information about the Trac mailing list