#905: NaNs in ADMBase::gxx using Development Version of ET
---------------------------------+------------------------------------------
Reporter: yosef@… | Type: defect
Status: new | Priority: major
Milestone: | Component: Other
Version: | Keywords:
---------------------------------+------------------------------------------
While testing the CCE in the ET, I stumbled across an issue
interpolating the admbase::metric at CCTK_POSTSTEP in GLOBAL mode
at particular spatial points.
The "bug" (or parfile/scheduling error) seems to be sensitive to
the exact grid setup. Basically, when I use the
qc0-mclachlan-CCE_Cauchy.par parfile, I get nans
for gxx. I made a small test thorn that only interpolates
gxx at a particular point. Nans show up at iteration 1.
attached is the parfile I used and the test thorn.
Here is the output from the above test thorn:
gxx(3.966899,0.000000,-17.630662) = 1.1160523626974641
at iter 0 on proc 1
gxx(3.966899,0.000000,-17.630662) = 1.1160523626974641
at iter 0 on proc 0
gxx(3.966899,0.000000,-17.630662) = -nan
at iter 8 on proc 1
gxx(3.966899,0.000000,-17.630662) = -nan
at iter 8 on proc 0
gxx(3.966899,0.000000,-17.630662) = -nan
at iter 16 on proc 1
gxx(3.966899,0.000000,-17.630662) = -nan
at iter 16 on proc 0
etc....
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/905>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#1376: Cactus loops endlessly if given the same name for aliased function and
actual implementation
--------------------+-------------------------------------------------------
Reporter: knarf | Owner:
Type: defect | Status: new
Priority: minor | Milestone:
Component: Cactus | Version: development version
Keywords: |
--------------------+-------------------------------------------------------
Currently, if an aliased function is declared in interface.ccl with
{{{
PROVIDES FUNCTION fun WITH fun LANGUAGE C
}}}
Cactus goes into an infinite loop when calling that function. It does so
because it creates a function 'fun' (the first) itself, which then calls
(fun) the second.
Instead, Cactus should already produce an error for above declaration. The
two 'fun' must be different for this to work, and Cactus should catch
this.
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/1376>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
#1375: Slab does not fille in outmost layer of points in mesh refinement levels
-----------------------------------+----------------------------------------
Reporter: rhaas | Owner:
Type: defect | Status: new
Priority: minor | Milestone:
Component: EinsteinToolkit thorn | Version: development version
Keywords: Slab |
-----------------------------------+----------------------------------------
This happens only for mesh refined runs and on reflevel>0 ie in the buffer
zones. This happens only if the "useghosts" option is not set (it is unset
by default). This behaviour is technically allowed by the docs since Slab
only claims to fill in ghost zones and the outermost points are not ghosts
but buffer zones (though Carpet disagrees somewhat on that in that it does
not count them in buffer_width).
Anyhow, this does not to my knowledge actually affect RotatingSymmetryXX
which are the most prominent users, at least their poison_boundaries
option does not seem to pick this up.
To test, run the attached modified slabtest.par test file whcih fails in
the identity test for level 1 unless useghosts=1 is used in slabtest.c
(and fails int the skip-every-second point test if it is used).
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/1375>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit