[ET Trac] #2534: Evolution of grid arrays with Mol no longer works

Roland Haas trac-noreply at einsteintoolkit.org
Fri Aug 6 16:11:38 CDT 2021


#2534: Evolution of grid arrays with Mol no longer works

 Reporter: 
   Status: new
Milestone: 
  Version: development version
     Type: bug
 Priority: major
Component: Carpet

Comment (by Roland Haas):

Writing my own parfile:

```
ActiveThorns = "Time MoL"
ActiveThorns = "Coordbase CartGrid3d Boundary StaticConformal"
ActiveThorns = "SymBase InitBase"
ActiveThorns = "IOUtil"
ActiveThorns = "Carpet CarpetLib CarpetRegrid2"

driver::ghost_size = 1
CoordBase::boundary_size_x_lower = 1
CoordBase::boundary_size_y_lower = 1
CoordBase::boundary_size_z_lower = 1
CoordBase::boundary_size_x_upper = 1
CoordBase::boundary_size_y_upper = 1
CoordBase::boundary_size_z_upper = 1

CartGrid3D::type      = "coordbase"
CoordBase::domainsize = "minmax"
CoordBase::xmin = -32.0
CoordBase::ymin = -32.0
CoordBase::zmin = -32.0
CoordBase::xmax =  32.0
CoordBase::ymax =  32.0
CoordBase::zmax =  32.0
CoordBase::dx   =   8.0
CoordBase::dy   =   8.0
CoordBase::dz   =   8.0

cactus::cctk_itlast = 1
Cactus::terminate           = "iteration"

Carpet::domain_from_coordbase = "yes"
Carpet::use_buffer_zones         = "yes"


Carpet::init_fill_timelevels     = "yes"

CarpetLib::poison_new_memory  = "yes"
CarpetLib::poison_value       = 114
Carpet::poison_new_timelevels = "yes"

Carpet::max_refinement_levels    = 2

CarpetRegrid2::regrid_every = 1
CarpetRegrid2::num_centres  = 1
CarpetRegrid2::num_levels_1 = 2
CarpetRegrid2::radius_1[1]  =16.0

time::dtfac = 0.25
MoL::ODE_Method             = "rk4"
MoL::MoL_Intermediate_Steps = 4
MoL::MoL_Num_Scratch_Levels = 1

ActiveThorns = EvolveArray
```

I can get a error similar to what is being described in the ticket.  
My guess is that this branch rhaas/poisonarray \([https://bitbucket.org/eschnett/carpet/branch/rhaas/poisonarray](https://bitbucket.org/eschnett/carpet/branch/rhaas/poisonarray)\) will fix things since it applies the same logic as is used otherwhere in PreSync.cc to [Poison.cc](http://Poison.cc) when it come to handling non-grid functions.

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2534/evolution-of-grid-arrays-with-mol-no
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/trac/attachments/20210806/aa0b60ca/attachment.html 


More information about the Trac mailing list