#2644: NaNs in ADMBASE metric components with RNSID
Reporter: Davide Guerra
Status: new
Milestone: ET_2022_11
Version: ET_2022_11
Type: bug
Priority: critical
Component: EinsteinToolkit thorn
I have found a problem related to a parfile that I am attaching here. I used the same parfile with the version of ET 2019\_10 “Mayer” having no problems at all, and the latest version of ET having the problem with the presence of NaNs in the ADMBASE:gij where ij are all the components of g. It happens after creating the initial data with HYDRO\_RNSID which I checked to be the same between this latest simulation and the previous one with the old version of ET.
Therefore, I attach here the parfile, the .out and the .err files that I used. They refers to the model D2 written in Table I studied by Luca Baiotti et al in the paper “[Accurate simulations of the dynamical bar-mode instability in full General Relativity](https://arxiv.org/pdf/astro-ph/0609473.pdf)”
Thank you in advance.
attachment: D2_test.out (https://api.bitbucket.org/2.0/repositories/einsteintoolkit/tickets/issues/2…)
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2644/nans-in-admbase-m…
#2642: Toolkit Fails to Compile on Sawtooth (Idaho National Lab)
Reporter: Leonardo Werneck
Status: new
Milestone:
Version: ET_2022_11
Type: bug
Priority: major
Component:
Comment (by Roland Haas):
This is the sympton caused by using an outdated stdc\+\+ template [library. In](http://library.In) you fix you are basically re-introducing the compatibility hacks \(more or less\) that we have just removed from Cactus. A better fix for the compilation failure is to make sure that icpc uses a new enough g\+\+ \(and its STL library\) by passing a `-gxx-name` option and the path to a new enough g\+\+.
If on the other hand you really want an isnan that persists even when -Ofast is used \(though you are defeating the Ofast purpose in that case and with Ofast there just may not be any NaN produce in a consitent way to begin with\) you could use your change. But it addresses a different issue than the compilation failure I think.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2642/toolkit-fails-to-…
#2642: Toolkit Fails to Compile on Sawtooth (Idaho National Lab)
Reporter: Leonardo Werneck
Status: new
Milestone:
Version: ET_2022_11
Type: bug
Priority: major
Component:
Comment (by Leonardo Werneck):
@{557058:8bc23f2a-45c0-477d-8ac4-a5a16c734278} @{557058:59e031ba-9bb5-4298-a472-7b99d0ae6f22} The fix to `IllinoisGRMHD` is proposed here: [https://bitbucket.org/zach\_etienne/wvuthorns/pull-requests/9](https://bitb…. Please review the request at your earliest convenience. This fix was already part of the version of IllinoisGRMHD that supports tabulated equations of state.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2642/toolkit-fails-to-…
#2642: Toolkit Fails to Compile on Sawtooth (Idaho National Lab)
Reporter: Leonardo Werneck
Status: new
Milestone:
Version: ET_2022_11
Type: bug
Priority: major
Component: Carpet
Compilation of the Riemann release of the toolkit fails with the Intel-19.0.5 compiler on Sawtooth \(Idaho National Lab\). Attached is the full configuration used for the compilation \(sawtooth-intel19.cfg\). The compilation errors are the following:
```
/home/rosaleon/sawtooth/test_compile/Cactus/arrangements/Carpet/CarpetMask/src/mask_surface.cc(144): error: more than one instance of overloaded function "CarpetMask::isnan" matches the argument list:
function "isnan(double)"
function "std::isnan(double)"
argument types are: (CCTK_REAL8)
assert(not isnan(theta));
^
/home/rosaleon/sawtooth/test_compile/Cactus/arrangements/Carpet/CarpetMask/src/mask_surface.cc(160): error: more than one instance of overloaded function "CarpetMask::isnan" matches the argument list:
function "isnan(double)"
function "std::isnan(double)"
argument types are: (CCTK_REAL8)
assert(not isnan(phi));
^
Checking status of thorn IllinoisGRMHD
compilation aborted for /home/rosaleon/sawtooth/test_compile/Cactus/configs/intel19/build/CarpetMask/mask_surface.cc (code 2)
make[4]: *** [mask_surface.cc.o] Error 2
make[3]: *** [make.checked] Error 2
make[2]: *** [/home/rosaleon/sawtooth/test_compile/Cactus/configs/intel19/lib/libthorn_CarpetMask.a] Error 2
make[2]: *** Waiting for unfinished jobs....
COMPILING WVUThorns/IllinoisGRMHD/src/driver_conserv_to_prims.C
/home/rosaleon/sawtooth/test_compile/Cactus/arrangements/WVUThorns/IllinoisGRMHD/src/driver_conserv_to_prims.C(204): error: more than one instance of overloaded function "isnan" matches the argument list:
function "isnan(double)"
function "std::isnan(double)"
argument types are: (double)
if(isnan(CONSERVS[RHOSTAR]*CONSERVS[STILDEX]*CONSERVS[STILDEY]*CONSERVS[STILDEZ]*CONSERVS[TAUENERGY]*PRIMS[BX_CENTER]*PRIMS[BY_CENTER]*PRIMS[BZ_CENTER])) {
^
compilation aborted for /home/rosaleon/sawtooth/test_compile/Cactus/configs/intel19/build/IllinoisGRMHD/driver_conserv_to_prims.C (code 2)
make[4]: *** [driver_conserv_to_prims.C.o] Error 2
make[3]: *** [make.checked] Error 2
make[2]: *** [/home/rosaleon/sawtooth/test_compile/Cactus/configs/intel19/lib/libthorn_IllinoisGRMHD.a] Error 2
make[1]: *** [intel19] Error 2
make: *** [default-target] Error 2
Using the only available configuration: intel19.
^Cmake[1]: *** [intel19] Interrupt
make: *** [default-target] Interrupt
```
I am currently working on a pull request to fix the issue in `IllinoisGRMHD`.
Cheers,
Leo
attachment: sawtooth-intel19.cfg (https://api.bitbucket.org/2.0/repositories/einsteintoolkit/tickets/issues/2…)
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2642/toolkit-fails-to-…
#2641: ML_BSSN macro usage bypasses personalized read/write macros
Reporter: Samuel Cupp
Status: new
Milestone:
Version: development version
Type: bug
Priority: trivial
Component: Kranc
The new macros are designed to only give access to the variables listed by read/write declarations, but ML\_BSSN bypasses this sanity check with its \*\_Body internal functions. As an example, in ML\_BSSN/src/ML\_BSSN\_ADMBaseInterior.cc is the scheduled function
```
extern "C" void ML_BSSN_ADMBaseInterior(CCTK_ARGUMENTS)
{
#ifdef DECLARE_CCTK_ARGUMENTS_ML_BSSN_ADMBaseInterior
DECLARE_CCTK_ARGUMENTS_CHECKED(ML_BSSN_ADMBaseInterior);
#else
DECLARE_CCTK_ARGUMENTS;
#endif
DECLARE_CCTK_PARAMETERS;
```
which then calls
```
static void ML_BSSN_ADMBaseInterior_Body(const cGH* restrict const cctkGH, const int dir, const int face, const CCTK_REAL normal[3], const CCTK_REAL tangentA[3], const CCTK_REAL tangentB[3], const int imin[3], const int imax[3], const int n_subblock_gfs, CCTK_REAL* restrict const subblock_gfs[])
{
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
```
The various \*\_Body should \(ideally\) either have their required variables passed explicitly or use the macro of the parent function. While the two macros don’t produce different numerical results, the personalized macros help validate the code’s behavior. It is admittedly a pretty low-priority task, but using the new macro is preferred when possible.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2641/ml_bssn-macro-usa…
#2640: NRPyEllipticET fails to compile on Stampede2
Reporter: Roland Haas
Status: resolved
Milestone:
Version:
Type: bug
Priority: blocker
Component:
Changes (by Roland Haas):
status: resolved (was new)
Comment (by Roland Haas):
this commit seems to fix the issue (at least on stampede2-skx). Thank you for the quick bugfix.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2640/nrpyellipticet-fa…
#2639: CarpetLib internal error with PreSync and analysis thorns
Reporter: Samuel Cupp
Status: new
Milestone:
Version: development version
Type: bug
Priority: major
Component: Carpet
Comment (by Steven R. Brandt):
One thing I sometimes try when nothing is making sense is to run the code through valgrind.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2639/carpetlib-interna…