msvcrt: Import erff implementation from musl.

Signed-off-by: Piotr Caban <piotr@codeweavers.com>
Signed-off-by: Alexandre Julliard <julliard@winehq.org>
This commit is contained in:
Piotr Caban 2021-05-19 15:26:29 +02:00 committed by Alexandre Julliard
parent de11fe6dcd
commit 98525b04bc
6 changed files with 41 additions and 20 deletions

1
configure vendored
View File

@ -19624,7 +19624,6 @@ for ac_func in \
atanhf \
erf \
erfc \
erff \
exp2 \
exp2f \
expm1 \

View File

@ -2664,7 +2664,6 @@ AC_CHECK_FUNCS(\
atanhf \
erf \
erfc \
erff \
exp2 \
exp2f \
expm1 \

View File

@ -5306,10 +5306,50 @@ static float erfc2f(UINT32 ix, float x)
/*********************************************************************
* erff (MSVCR120.@)
*
* Copied from musl: src/math/erff.c
*/
float CDECL erff(float x)
{
return unix_funcs->erff( x );
static const float efx8 = 1.0270333290e+00,
pp0 = 1.2837916613e-01,
pp1 = -3.2504209876e-01,
pp2 = -2.8481749818e-02,
pp3 = -5.7702702470e-03,
pp4 = -2.3763017452e-05,
qq1 = 3.9791721106e-01,
qq2 = 6.5022252500e-02,
qq3 = 5.0813062117e-03,
qq4 = 1.3249473704e-04,
qq5 = -3.9602282413e-06;
float r, s, z, y;
UINT32 ix;
int sign;
ix = *(UINT32*)&x;
sign = ix >> 31;
ix &= 0x7fffffff;
if (ix >= 0x7f800000) {
/* erf(nan)=nan, erf(+-inf)=+-1 */
return 1 - 2 * sign + 1 / x;
}
if (ix < 0x3f580000) { /* |x| < 0.84375 */
if (ix < 0x31800000) { /* |x| < 2**-28 */
/*avoid underflow */
return 0.125f * (8 * x + efx8 * x);
}
z = x * x;
r = pp0 + z * (pp1 + z * (pp2 + z * (pp3 + z * pp4)));
s = 1 + z * (qq1 + z * (qq2 + z * (qq3 + z * (qq4 + z * qq5))));
y = r / s;
return x + x * y;
}
if (ix < 0x40c00000) /* |x| < 6 */
y = 1 - erfc2f(ix, x);
else
y = 1 - FLT_MIN;
return sign ? -y : y;
}
/*********************************************************************

View File

@ -173,18 +173,6 @@ static double CDECL unix_erf(double x)
#endif
}
/*********************************************************************
* erff
*/
static float CDECL unix_erff(float x)
{
#ifdef HAVE_ERFF
return erff(x);
#else
return unix_erf(x);
#endif
}
/*********************************************************************
* erfc
*/
@ -567,7 +555,6 @@ static const struct unix_funcs funcs =
unix_coshf,
unix_erf,
unix_erfc,
unix_erff,
unix_exp,
unix_expf,
unix_exp2,

View File

@ -35,7 +35,6 @@ struct unix_funcs
float (CDECL *coshf)(float x);
double (CDECL *erf)(double x);
double (CDECL *erfc)(double x);
float (CDECL *erff)(float x);
double (CDECL *exp)(double x);
float (CDECL *expf)(float x);
double (CDECL *exp2)(double x);

View File

@ -116,9 +116,6 @@
/* Define to 1 if you have the `erfc' function. */
#undef HAVE_ERFC
/* Define to 1 if you have the `erff' function. */
#undef HAVE_ERFF
/* Define to 1 if you have the `exp2' function. */
#undef HAVE_EXP2