#2582: WVUThorns use CCTK_GFINDEX4D when they should be using CCTK_VECTGFINDEX#D
Reporter:
Status: new
Milestone:
Version: development version
Type: bug
Priority: major
Component:
Eg in Seed_Magnetic_Field:
```
Avec[CCTK_GFINDEX4D(cctkGH,i,j,k,0)] = -(yL + 0.5*dY)*A_b*pow(MAX(Pressure_at_
```
this causes test failures when compiled with `DEBUG=yes` or when enabling grid function eidex checking via other means.
Affected thorns seem to be:
* magnetizedTOV-Baikal (from Baikal)
* GiRaFFE_tests_FFEBreak (from GiRaFFE)
* GiRaFFE_tests_FastWave (from GiRaFFE)
* GiRaFFE_tests_ThreeWave (from GiRaFFE)
* magnetizedTOV (from IllinoisGRMHD)
* magnetized_explosionTOV (from VolumeIntegrals_GRMHD)
* magnetized_explosionTOV (from VolumeIntegrals_vacuum)
* magnetized_explosionTOV (from particle_tracerET)
Errors in the log file look like this:
```
INFO (Seed_Magnetic_Fields): Seeding magnetic fields
ERROR from host ekohaes8 process 0
while executing schedule bin CCTK_INITIAL, routine Seed_Magnetic_Fields::Seed_Magnetic_Fields
in thorn Cactus, file /data/rhaas/postdoc/gr/cactus/ET_trunk/configs/sim/build/Cactus/main/DebugDefines.c:105:
-> Grid function index out of bounds. i=0 j=0 k=0 l=0 cctk_lsh=[29,29,18,0]
ERROR from host ekohaes8 process 0
while executing schedule bin CCTK_INITIAL, routine Seed_Magnetic_Fields::Seed_Magnetic_Fields
in thorn Cactus, file /data/rhaas/postdoc/gr/cactus/ET_trunk/configs/sim/build/Cactus/main/DebugDefines.c:105:
-> Grid function index out of bounds. i=0 j=0 k=0 l=0 cctk_lsh=[29,29,18,0]
cactus_sim: /data/rhaas/postdoc/gr/cactus/ET_trunk/configs/sim/build/Carpet/helpers.cc:275: int Carpet::Abort(const cGH*, int): Assertion `0' failed.
Rank 0 with PID 648211 received signal 6
Writing backtrace to magnetizedTOV-Baikal/backtrace.0.txt
```
```
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2582/wvuthorns-use-cct…
#2578: Kranc, Cactus, and MPI linking issues
Reporter: Nicholas Hugh Olsen
Status: new
Milestone:
Version: ET_2021_05
Type: bug
Priority: minor
Component: Cactus
Comment (by Ian Hinder):
@{606c600bedc14f007689694e} Is this still an issue, or did you manage to fix it? It sounds like it might be a problem with having the wrong, or multiple MPI libraries installed.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2578/kranc-cactus-and-…
#2581: Vectors test fails on AMD Ryzen 9 4900HS
Reporter: Zach Etienne
Status: new
Milestone:
Version: development version
Type: bug
Priority: major
Component: EinsteinToolkit thorn
Vectors, after correctly finding that CCTK\_REAL\_VEC\_SIZE is in fact 4 on this CPU, with flags:
```
flags : fpu vme de pse tsc msr pae mce cx8 apic sep mtrr pge mca cmov pat pse36 clflush mmx fxsr sse sse2 ht syscall nx mmxext fxsr_opt pdpe1gb rdtscp lm constant_tsc rep_good nopl nonstop_tsc cpuid extd_apicid aperfmperf pni pclmulqdq monitor ssse3 fma cx16 sse4_1 sse4_2 movbe popcnt aes xsave avx f16c rdrand lahf_lm cmp_legacy svm extapic cr8_legacy abm sse4a misalignsse 3dnowprefetch osvw ibs skinit wdt tce topoext perfctr_core perfctr_nb bpext perfctr_llc mwaitx cpb cat_l3 cdp_l3 hw_pstate ssbd mba ibrs ibpb stibp vmmcall fsgsbase bmi1 avx2 smep bmi2 cqm rdt_a rdseed adx smap clflushopt clwb sha_ni xsaveopt xsavec xgetbv1 xsaves cqm_llc cqm_occup_llc cqm_mbm_total cqm_mbm_local clzero irperf xsaveerptr rdpru wbnoinvd arat npt lbrv svm_lock nrip_save tsc_scale vmcb_clean flushbyasid decodeassists pausefilter pfthreshold avic v_vmsave_vmload vgif v_spec_ctrl umip rdpid overflow_recov succor smca sme sev sev_es
```
segfaults on the vec\_loadu\_maybe\[1\] test in Vectors.
```
INFO (Vectors): Testing vectorisation... [errors may result in segfaults]
INFO (VecINFO (Vectors): Testing vectorisation... [errors may result in segfaults]
INFO (Vectors): Test vec_set1...
INFO (Vectors): Test vec_set...
INFO (Vectors): Test vec_elt[0]...
INFO (Vectors): Test vec_elt[1]...
INFO (Vectors): Test vec_elt[2]...
INFO (Vectors): Test vec_elt[3]...
INFO (Vectors): Test vec_load...
INFO (Vectors): Test vec_loadu[0]...
INFO (Vectors): Test vec_loadu[1]...
INFO (Vectors): Test vec_loadu[2]...
INFO (Vectors): Test vec_loadu[3]...
INFO (Vectors): Test vec_loadu_maybe[0]...
INFO (Vectors): Test vec_loadu_maybe[1]...
Rank 0 with PID 218092 received signal 11
```
After commenting out this test, the test suite runs with no problems.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2581/vectors-test-fail…
#2554: init_3_timelevels produces nans on past timelevels of the initial slice
Reporter: Yosef Zlochower
Status: open
Milestone:
Version: ET_2021_05
Type: bug
Priority: minor
Component: Carpet
Comment (by Roland Haas):
I am actually not sure if this ever worked as advertised. The code in question is:
```c++
initialise_3tl_flip_timelevels(cctkGH);
initialise_3tl_evolve(cctkGH);
initialise_3tl_evolve(cctkGH);
// TODO: May want to restrict where possible (i.e. if the time
// refinement factor is one)
initialise_3tl_recycle(cctkGH);
```
which is bad because `initialise_3tl_evolve` steps all reflevel by one step, so that in step 1 when say reflevel 1 wants to do prolongation for the buffer zones there are only 2 valid timelevels on reflevel 0 \(namley tl=0 and tl=1\) and thus interpolation in time picks up nan from tl=2.
A fix is proposed in: [https://bitbucket.org/eschnett/carpet/commits/7147f842569c37719d1c97d5f6ddd… essentially moving the reflevel loop outside of `initialise_3tl_evolve`. Still has issues with the restriction involved but I don’t think those can be avoided.
Comments would be welcome.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2554/init_3_timelevels…
#2173: Test "Poisson equation" example
Reporter: Roland Haas
Status: resolved
Milestone: ET_2021_11
Version: development version
Type: task
Priority: major
Component: EinsteinToolkit website
Changes (by Miguel Gracia):
status: resolved (was open)
Comment (by Miguel Gracia):
Done for the Johnson release.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2173/test-poisson-equa…
#2174: Test "Multi Patch Scalar Wave Equation" example
Reporter: Roland Haas
Status: resolved
Milestone: ET_2021_11
Version: development version
Type: task
Priority: major
Component: EinsteinToolkit website
Changes (by William Gabella):
status: resolved (was open)
Comment (by William Gabella):
Finished for this relase, Dec 2021.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2174/test-multi-patch-…
#2175: Test "Single, stable neutron star" example
Reporter: Roland Haas
Status: resolved
Milestone: ET_2021_11
Version: development version
Type: task
Priority: major
Component: EinsteinToolkit website
Changes (by Peter Diener):
status: resolved (was open)
Comment (by Peter Diener):
Done for the November 2021 release.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2175/test-single-stabl…