|
|
@ -237,6 +237,7 @@ double TMF1 , TMF2 , TMF3 , TMF4 ; |
|
|
* smoothZero: flooring to zero. |
|
|
* smoothZero: flooring to zero. |
|
|
* y = 0.5 ( x + sqrt( x^2 + 4 delta^2 ) ) |
|
|
* y = 0.5 ( x + sqrt( x^2 + 4 delta^2 ) ) |
|
|
*-----------------*/ |
|
|
*-----------------*/ |
|
|
|
|
|
#if 0 |
|
|
static double smoothZero |
|
|
static double smoothZero |
|
|
( |
|
|
( |
|
|
double x, |
|
|
double x, |
|
|
@ -248,6 +249,7 @@ static double smoothZero |
|
|
if (dx) *dx = 0.5 * ( 1.0 + x / sqr ) ; |
|
|
if (dx) *dx = 0.5 * ( 1.0 + x / sqr ) ; |
|
|
return 0.5 * ( x + sqr ) ; |
|
|
return 0.5 * ( x + sqr ) ; |
|
|
} |
|
|
} |
|
|
|
|
|
#endif |
|
|
/*---------------------------------------------------* |
|
|
/*---------------------------------------------------* |
|
|
* CeilingPow: ceiling for positive x, flooring for negative x. |
|
|
* CeilingPow: ceiling for positive x, flooring for negative x. |
|
|
* y = x * xmax / ( x^{2m} + xmax^{2m} )^{1/(2m)} |
|
|
* y = x * xmax / ( x^{2m} + xmax^{2m} )^{1/(2m)} |
|
|
@ -256,6 +258,7 @@ static double smoothZero |
|
|
* - -xmax < y < xmax. |
|
|
* - -xmax < y < xmax. |
|
|
* - dy/dx|_{x=0} = 1. |
|
|
* - dy/dx|_{x=0} = 1. |
|
|
*-----------------*/ |
|
|
*-----------------*/ |
|
|
|
|
|
#if 0 |
|
|
static double CeilingPow |
|
|
static double CeilingPow |
|
|
( |
|
|
( |
|
|
double x, |
|
|
double x, |
|
|
@ -298,6 +301,7 @@ static double CeilingPow |
|
|
(*dx) = xmax * xmp * dnm / arg ; |
|
|
(*dx) = xmax * xmp * dnm / arg ; |
|
|
return result ; |
|
|
return result ; |
|
|
} |
|
|
} |
|
|
|
|
|
#endif |
|
|
|
|
|
|
|
|
/*---------------------------------------------------* |
|
|
/*---------------------------------------------------* |
|
|
* CeilingPow: ceiling for positive x, flooring for negative x. |
|
|
* CeilingPow: ceiling for positive x, flooring for negative x. |
|
|
|