[ET Trac] [Einstein Toolkit] #610: IDAnalyticBH: Smooth out singularities

[Einstein Toolkit]

#610: IDAnalyticBH: Smooth out singularities
-----------------------------------+----------------------------------------
 Reporter:  eschnett               |       Owner:     
     Type:  defect                 |      Status:  new
 Priority:  major                  |   Milestone:     
Component:  EinsteinToolkit thorn  |     Version:     
 Keywords:                         |  
-----------------------------------+----------------------------------------
 IDAnalyticBH has an "epsilon" parameter to smooth out singularites. This
 parameter should also be used in Brill-Lindquist initial data. Otherwise,
 moving puncture evolution of these data is not possible without an
 explicit smoothing step.

 {{{
 Index: src/BrillLindquist.c
 ===================================================================
 --- src/BrillLindquist.c        (revision 180)
 +++ src/BrillLindquist.c        (working copy)
 @@ -143,7 +143,7 @@
        tmp1 = sqrt(SQR(xval+hole_x0[n])
                   +SQR(yval+hole_y0[n])
                   +SQR(zval+hole_z0[n])
 -                 +1.0e-20);
 +                 +SQR(epsilon));

        psi[i] += hole_mass[n]/tmp1*0.5;

 }}}

-- 
Ticket URL: <https://trac.einsteintoolkit.org/ticket/610>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit