[ET Trac] [Einstein Toolkit] #905: NaNs in ADMBase::gxx using Development Version of ET

Einstein Toolkit trac-noreply at einsteintoolkit.org
Wed May 16 16:09:11 CDT 2012


#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


More information about the Trac mailing list