[ET Trac] [Einstein Toolkit] #1429: Assertion error when using "eval" & new UIUC speedup in TwoPunctures

Einstein Toolkit trac-noreply at einsteintoolkit.org
Wed Sep 4 23:01:36 CDT 2013

#1429: Assertion error when using "eval" & new UIUC speedup in TwoPunctures
  Reporter:  bernard.j.kelly@…         |       Owner:                      
      Type:  defect                    |      Status:  new                 
  Priority:  major                     |   Milestone:  ET_2013_11          
 Component:  EinsteinToolkit thorn     |     Version:  development version 
Resolution:                            |    Keywords:  TwoPunctures, malloc

Comment (by bernard.j.kelly@…):

 Hi Roland. Apologies for the CCTK_REAL/CCTK_REAL* typeof; I'd corrected
 that myself locally, or so I thought.

 Anyway, I can confirm that fixing the one obvious logic bug in line 115 of
 the original code does *not* fix the actual issue; I'd tried that first in
 my own run. But it should be fixed, anyway, if the Numerical Recipes code
 is to persist.

 I suspect Erik's observation --- "It seems the problem is that
 retval[nrl][0] is never set." --- is the real problem in the original
 code. Nevertheless, I think the Numerical Recipes Fortran-masquerading-
 as-C code should be junked anyway. I was originally going to submit my
 simpler d3tensor as a replacement for what's there, but then thought that
 this would be more minimalist. That still leaves "dvector" etc ...

Ticket URL: <https://trac.einsteintoolkit.org/ticket/1429#comment:5>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit

More information about the Trac mailing list