Skip to content

Commit

Permalink
'inverf.c' : debugging code should have been enclosed within #ifdef D…
Browse files Browse the repository at this point in the history
…EBUGGING/#endif
  • Loading branch information
Bill-Gray committed Aug 7, 2023
1 parent c765ddf commit 6d83af4
Showing 1 changed file with 2 additions and 0 deletions.
2 changes: 2 additions & 0 deletions inverf.c
Expand Up @@ -47,10 +47,12 @@ double inverf( const double y)
{
double x, diff;
const double erf_limit = .915;
#ifdef DEBUGGING
const double t = sqrt( -2. * log( (1. - y) / 2.));
const double guess = -0.70711*((2.30753+t*0.27061)/(1.+t*(0.99229+t*0.04481)) - t);

printf( "NR guess %.8f\n", guess);
#endif
if( y < -erf_limit)
return( -inverf( -y));
else if( y > erf_limit)
Expand Down

0 comments on commit 6d83af4

Please sign in to comment.