Hi All,
I am working with some piece of code based on a multipatch coordinate system (llama, http://dx.doi.org/10.1103/PhysRevD.83.044045). Currently, I'm interested in doing a test in a cubical kind of domain, that is, without the multipatch system. Closely following the qc0-mclachlan.par parameter file (ET_2014_05), I converted my multipatch parameter file to use the following: Carpet::domain_from_multipatch = "no" Carpet::domain_from_coordbase = yes CartGrid3D::type = "coordbase"
Unlike when using a multipatch sysytem, I've had to add the ReflectionSymmetry thorn in my thornlist for running on a cubical grid. The simulation is able to build successfully. But when running, it breaks with the following error from Carpet: cactus/arrangements/Carpet/Carpet/src/SetupGH.cc:1744: void Carpet::get_domain_specification(const cGH*, int, const ivect&, rvect&, rvect&, rvect&): Assertion `not CCTK_IsFunctionAliased ("MultiPatch_GetDomainSpecification")' failed.
Looking at the function that fails and a few ideas I got from the user mailing list, my naive thought is that this might possibly be an issue of cctk_lsh vs cctk_ash in Carpet. (e.g. http://lists.einsteintoolkit.org/pipermail/users/2014-October/003849.html). But I'm yet to test this idea. The code I'm working with uses cctk_lsh.
I would gladly appreciate help from anyone who could have an idea about this.
I'm using ET_2014_05 release of the EinsteinToolkit/Cactus. I have attached my test parameter file for viewing in case it will be useful in diagnosing the problem(s).
Thank you in advance for your kind assistance.
Regards, Dumsani