1263 lines
30 KiB
C
1263 lines
30 KiB
C
/*
|
|
* msvcrt.dll math functions
|
|
*
|
|
* Copyright 2000 Jon Griffiths
|
|
*
|
|
* This library is free software; you can redistribute it and/or
|
|
* modify it under the terms of the GNU Lesser General Public
|
|
* License as published by the Free Software Foundation; either
|
|
* version 2.1 of the License, or (at your option) any later version.
|
|
*
|
|
* This library is distributed in the hope that it will be useful,
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
* Lesser General Public License for more details.
|
|
*
|
|
* You should have received a copy of the GNU Lesser General Public
|
|
* License along with this library; if not, write to the Free Software
|
|
* Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA
|
|
*/
|
|
#include "config.h"
|
|
|
|
#include <stdio.h>
|
|
#define __USE_ISOC9X 1
|
|
#define __USE_ISOC99 1
|
|
#include <math.h>
|
|
#ifdef HAVE_IEEEFP_H
|
|
#include <ieeefp.h>
|
|
#endif
|
|
|
|
#include "msvcrt.h"
|
|
|
|
#include "wine/debug.h"
|
|
|
|
WINE_DEFAULT_DEBUG_CHANNEL(msvcrt);
|
|
|
|
#ifndef HAVE_FINITE
|
|
#ifndef finite /* Could be a macro */
|
|
#ifdef isfinite
|
|
#define finite(x) isfinite(x)
|
|
#else
|
|
#define finite(x) (!isnan(x)) /* At least catch some cases */
|
|
#endif
|
|
#endif
|
|
#endif
|
|
|
|
#ifndef signbit
|
|
#define signbit(x) 0
|
|
#endif
|
|
|
|
typedef int (*MSVCRT_matherr_func)(struct MSVCRT__exception *);
|
|
|
|
static MSVCRT_matherr_func MSVCRT_default_matherr_func = NULL;
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_acos (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_acos( double x )
|
|
{
|
|
if (x < -1.0 || x > 1.0 || !finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
/* glibc implements acos() as the FPU equivalent of atan2(sqrt(1 - x ^ 2), x).
|
|
* asin() uses a similar construction. This is bad because as x gets nearer to
|
|
* 1 the error in the expression "1 - x^2" can get relatively large due to
|
|
* cancellation. The sqrt() makes things worse. A safer way to calculate
|
|
* acos() is to use atan2(sqrt((1 - x) * (1 + x)), x). */
|
|
return atan2(sqrt((1 - x) * (1 + x)), x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_asin (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_asin( double x )
|
|
{
|
|
if (x < -1.0 || x > 1.0 || !finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return atan2(x, sqrt((1 - x) * (1 + x)));
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_atan (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_atan( double x )
|
|
{
|
|
if (!finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return atan(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_atan2 (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_atan2( double x, double y )
|
|
{
|
|
if (!finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return atan2(x,y);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_cos (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_cos( double x )
|
|
{
|
|
if (!finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return cos(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_cosh (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_cosh( double x )
|
|
{
|
|
if (!finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return cosh(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_exp (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_exp( double x )
|
|
{
|
|
if (!finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return exp(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_fmod (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_fmod( double x, double y )
|
|
{
|
|
if (!finite(x) || !finite(y)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return fmod(x,y);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_log (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_log( double x)
|
|
{
|
|
if (x < 0.0 || !finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
if (x == 0.0) *MSVCRT__errno() = MSVCRT_ERANGE;
|
|
return log(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_log10 (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_log10( double x )
|
|
{
|
|
if (x < 0.0 || !finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
if (x == 0.0) *MSVCRT__errno() = MSVCRT_ERANGE;
|
|
return log10(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_pow (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_pow( double x, double y )
|
|
{
|
|
/* FIXME: If x < 0 and y is not integral, set EDOM */
|
|
double z = pow(x,y);
|
|
if (!finite(z)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return z;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_sin (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_sin( double x )
|
|
{
|
|
if (!finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return sin(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_sinh (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_sinh( double x )
|
|
{
|
|
if (!finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return sinh(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_sqrt (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_sqrt( double x )
|
|
{
|
|
if (x < 0.0 || !finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return sqrt(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_tan (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_tan( double x )
|
|
{
|
|
if (!finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return tan(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* MSVCRT_tanh (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_tanh( double x )
|
|
{
|
|
if (!finite(x)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return tanh(x);
|
|
}
|
|
|
|
|
|
#if defined(__GNUC__) && defined(__i386__)
|
|
|
|
#define FPU_DOUBLE(var) double var; \
|
|
__asm__ __volatile__( "fstpl %0;fwait" : "=m" (var) : )
|
|
#define FPU_DOUBLES(var1,var2) double var1,var2; \
|
|
__asm__ __volatile__( "fstpl %0;fwait" : "=m" (var2) : ); \
|
|
__asm__ __volatile__( "fstpl %0;fwait" : "=m" (var1) : )
|
|
|
|
/*********************************************************************
|
|
* _CIacos (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIacos(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_acos(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIasin (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIasin(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_asin(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIatan (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIatan(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_atan(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIatan2 (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIatan2(void)
|
|
{
|
|
FPU_DOUBLES(x,y);
|
|
return MSVCRT_atan2(x,y);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIcos (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIcos(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_cos(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIcosh (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIcosh(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_cosh(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIexp (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIexp(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_exp(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIfmod (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIfmod(void)
|
|
{
|
|
FPU_DOUBLES(x,y);
|
|
return MSVCRT_fmod(x,y);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIlog (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIlog(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_log(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIlog10 (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIlog10(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_log10(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIpow (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIpow(void)
|
|
{
|
|
FPU_DOUBLES(x,y);
|
|
return MSVCRT_pow(x,y);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIsin (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIsin(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_sin(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIsinh (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIsinh(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_sinh(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CIsqrt (MSVCRT.@)
|
|
*/
|
|
double CDECL _CIsqrt(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_sqrt(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CItan (MSVCRT.@)
|
|
*/
|
|
double CDECL _CItan(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_tan(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _CItanh (MSVCRT.@)
|
|
*/
|
|
double CDECL _CItanh(void)
|
|
{
|
|
FPU_DOUBLE(x);
|
|
return MSVCRT_tanh(x);
|
|
}
|
|
|
|
#else /* defined(__GNUC__) && defined(__i386__) */
|
|
|
|
/* The above cannot be called on non x86 platforms, stub them for linking */
|
|
|
|
#define IX86_ONLY(func) double func(void) { return 0.0; }
|
|
|
|
IX86_ONLY(_CIacos)
|
|
IX86_ONLY(_CIasin)
|
|
IX86_ONLY(_CIatan)
|
|
IX86_ONLY(_CIatan2)
|
|
IX86_ONLY(_CIcos)
|
|
IX86_ONLY(_CIcosh)
|
|
IX86_ONLY(_CIexp)
|
|
IX86_ONLY(_CIfmod)
|
|
IX86_ONLY(_CIlog)
|
|
IX86_ONLY(_CIlog10)
|
|
IX86_ONLY(_CIpow)
|
|
IX86_ONLY(_CIsin)
|
|
IX86_ONLY(_CIsinh)
|
|
IX86_ONLY(_CIsqrt)
|
|
IX86_ONLY(_CItan)
|
|
IX86_ONLY(_CItanh)
|
|
|
|
#endif /* defined(__GNUC__) && defined(__i386__) */
|
|
|
|
/*********************************************************************
|
|
* _fpclass (MSVCRT.@)
|
|
*/
|
|
int CDECL _fpclass(double num)
|
|
{
|
|
#if defined(HAVE_FPCLASS) || defined(fpclass)
|
|
switch (fpclass( num ))
|
|
{
|
|
#ifdef FP_SNAN
|
|
case FP_SNAN: return MSVCRT__FPCLASS_SNAN;
|
|
#endif
|
|
#ifdef FP_QNAN
|
|
case FP_QNAN: return MSVCRT__FPCLASS_QNAN;
|
|
#endif
|
|
#ifdef FP_NINF
|
|
case FP_NINF: return MSVCRT__FPCLASS_NINF;
|
|
#endif
|
|
#ifdef FP_PINF
|
|
case FP_PINF: return MSVCRT__FPCLASS_PINF;
|
|
#endif
|
|
#ifdef FP_NDENORM
|
|
case FP_NDENORM: return MSVCRT__FPCLASS_ND;
|
|
#endif
|
|
#ifdef FP_PDENORM
|
|
case FP_PDENORM: return MSVCRT__FPCLASS_PD;
|
|
#endif
|
|
#ifdef FP_NZERO
|
|
case FP_NZERO: return MSVCRT__FPCLASS_NZ;
|
|
#endif
|
|
#ifdef FP_PZERO
|
|
case FP_PZERO: return MSVCRT__FPCLASS_PZ;
|
|
#endif
|
|
#ifdef FP_NNORM
|
|
case FP_NNORM: return MSVCRT__FPCLASS_NN;
|
|
#endif
|
|
#ifdef FP_PNORM
|
|
case FP_PNORM: return MSVCRT__FPCLASS_PN;
|
|
#endif
|
|
default: return MSVCRT__FPCLASS_PN;
|
|
}
|
|
#elif defined (fpclassify)
|
|
switch (fpclassify( num ))
|
|
{
|
|
case FP_NAN: return MSVCRT__FPCLASS_QNAN;
|
|
case FP_INFINITE: return signbit(num) ? MSVCRT__FPCLASS_NINF : MSVCRT__FPCLASS_PINF;
|
|
case FP_SUBNORMAL: return signbit(num) ?MSVCRT__FPCLASS_ND : MSVCRT__FPCLASS_PD;
|
|
case FP_ZERO: return signbit(num) ? MSVCRT__FPCLASS_NZ : MSVCRT__FPCLASS_PZ;
|
|
}
|
|
return signbit(num) ? MSVCRT__FPCLASS_NN : MSVCRT__FPCLASS_PN;
|
|
#else
|
|
if (!finite(num))
|
|
return MSVCRT__FPCLASS_QNAN;
|
|
return num == 0.0 ? MSVCRT__FPCLASS_PZ : (num < 0 ? MSVCRT__FPCLASS_NN : MSVCRT__FPCLASS_PN);
|
|
#endif
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _rotl (MSVCRT.@)
|
|
*/
|
|
unsigned int CDECL _rotl(unsigned int num, int shift)
|
|
{
|
|
shift &= 31;
|
|
return (num << shift) | (num >> (32-shift));
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _logb (MSVCRT.@)
|
|
*/
|
|
double CDECL _logb(double num)
|
|
{
|
|
if (!finite(num)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return logb(num);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _lrotl (MSVCRT.@)
|
|
*/
|
|
unsigned long CDECL _lrotl(unsigned long num, int shift)
|
|
{
|
|
shift &= 0x1f;
|
|
return (num << shift) | (num >> (32-shift));
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _lrotr (MSVCRT.@)
|
|
*/
|
|
unsigned long CDECL _lrotr(unsigned long num, int shift)
|
|
{
|
|
shift &= 0x1f;
|
|
return (num >> shift) | (num << (32-shift));
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _rotr (MSVCRT.@)
|
|
*/
|
|
unsigned int CDECL _rotr(unsigned int num, int shift)
|
|
{
|
|
shift &= 0x1f;
|
|
return (num >> shift) | (num << (32-shift));
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _scalb (MSVCRT.@)
|
|
*/
|
|
double CDECL _scalb(double num, long power)
|
|
{
|
|
/* Note - Can't forward directly as libc expects y as double */
|
|
double dblpower = (double)power;
|
|
if (!finite(num)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
return scalb(num, dblpower);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _hypot (MSVCRT.@)
|
|
*/
|
|
double CDECL _hypot(double x, double y)
|
|
{
|
|
/* FIXME: errno handling */
|
|
return hypot( x, y );
|
|
}
|
|
|
|
/*********************************************************************
|
|
* ceil (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_ceil( double x )
|
|
{
|
|
return ceil(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* floor (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_floor( double x )
|
|
{
|
|
return floor(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* fabs (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_fabs( double x )
|
|
{
|
|
return fabs(x);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* frexp (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_frexp( double x, int *exp )
|
|
{
|
|
return frexp( x, exp );
|
|
}
|
|
|
|
/*********************************************************************
|
|
* modf (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_modf( double x, double *iptr )
|
|
{
|
|
return modf( x, iptr );
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _matherr (MSVCRT.@)
|
|
*/
|
|
int CDECL MSVCRT__matherr(struct MSVCRT__exception *e)
|
|
{
|
|
if (e)
|
|
TRACE("(%p = %d, %s, %g %g %g)\n",e, e->type, e->name, e->arg1, e->arg2,
|
|
e->retval);
|
|
else
|
|
TRACE("(null)\n");
|
|
if (MSVCRT_default_matherr_func)
|
|
return MSVCRT_default_matherr_func(e);
|
|
ERR(":Unhandled math error!\n");
|
|
return 0;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* __setusermatherr (MSVCRT.@)
|
|
*/
|
|
void CDECL MSVCRT___setusermatherr(MSVCRT_matherr_func func)
|
|
{
|
|
MSVCRT_default_matherr_func = func;
|
|
TRACE(":new matherr handler %p\n", func);
|
|
}
|
|
|
|
/**********************************************************************
|
|
* _statusfp (MSVCRT.@)
|
|
*/
|
|
unsigned int CDECL _statusfp(void)
|
|
{
|
|
unsigned int retVal = 0;
|
|
#if defined(__GNUC__) && defined(__i386__)
|
|
unsigned int fpword;
|
|
|
|
__asm__ __volatile__( "fstsw %0" : "=m" (fpword) : );
|
|
if (fpword & 0x1) retVal |= MSVCRT__SW_INVALID;
|
|
if (fpword & 0x2) retVal |= MSVCRT__SW_DENORMAL;
|
|
if (fpword & 0x4) retVal |= MSVCRT__SW_ZERODIVIDE;
|
|
if (fpword & 0x8) retVal |= MSVCRT__SW_OVERFLOW;
|
|
if (fpword & 0x10) retVal |= MSVCRT__SW_UNDERFLOW;
|
|
if (fpword & 0x20) retVal |= MSVCRT__SW_INEXACT;
|
|
#else
|
|
FIXME(":Not implemented!\n");
|
|
#endif
|
|
return retVal;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _clearfp (MSVCRT.@)
|
|
*/
|
|
unsigned int CDECL _clearfp(void)
|
|
{
|
|
unsigned int retVal = _statusfp();
|
|
#if defined(__GNUC__) && defined(__i386__)
|
|
__asm__ __volatile__( "fnclex" );
|
|
#else
|
|
FIXME(":Not Implemented\n");
|
|
#endif
|
|
return retVal;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* __fpecode (MSVCRT.@)
|
|
*/
|
|
int * CDECL __fpecode(void)
|
|
{
|
|
return &msvcrt_get_thread_data()->fpecode;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* ldexp (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT_ldexp(double num, long exp)
|
|
{
|
|
double z = ldexp(num,exp);
|
|
|
|
if (!finite(z))
|
|
*MSVCRT__errno() = MSVCRT_ERANGE;
|
|
else if (z == 0 && signbit(z))
|
|
z = 0.0; /* Convert -0 -> +0 */
|
|
return z;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _cabs (MSVCRT.@)
|
|
*/
|
|
double CDECL MSVCRT__cabs(struct MSVCRT__complex num)
|
|
{
|
|
return sqrt(num.x * num.x + num.y * num.y);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _chgsign (MSVCRT.@)
|
|
*/
|
|
double CDECL _chgsign(double num)
|
|
{
|
|
/* FIXME: +-infinity,Nan not tested */
|
|
return -num;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _control87 (MSVCRT.@)
|
|
*/
|
|
unsigned int CDECL _control87(unsigned int newval, unsigned int mask)
|
|
{
|
|
#if defined(__GNUC__) && defined(__i386__)
|
|
unsigned int fpword = 0;
|
|
unsigned int flags = 0;
|
|
|
|
TRACE("(%08x, %08x): Called\n", newval, mask);
|
|
|
|
/* Get fp control word */
|
|
__asm__ __volatile__( "fstcw %0" : "=m" (fpword) : );
|
|
|
|
TRACE("Control word before : %08x\n", fpword);
|
|
|
|
/* Convert into mask constants */
|
|
if (fpword & 0x1) flags |= MSVCRT__EM_INVALID;
|
|
if (fpword & 0x2) flags |= MSVCRT__EM_DENORMAL;
|
|
if (fpword & 0x4) flags |= MSVCRT__EM_ZERODIVIDE;
|
|
if (fpword & 0x8) flags |= MSVCRT__EM_OVERFLOW;
|
|
if (fpword & 0x10) flags |= MSVCRT__EM_UNDERFLOW;
|
|
if (fpword & 0x20) flags |= MSVCRT__EM_INEXACT;
|
|
switch(fpword & 0xC00) {
|
|
case 0xC00: flags |= MSVCRT__RC_UP|MSVCRT__RC_DOWN; break;
|
|
case 0x800: flags |= MSVCRT__RC_UP; break;
|
|
case 0x400: flags |= MSVCRT__RC_DOWN; break;
|
|
}
|
|
switch(fpword & 0x300) {
|
|
case 0x0: flags |= MSVCRT__PC_24; break;
|
|
case 0x200: flags |= MSVCRT__PC_53; break;
|
|
case 0x300: flags |= MSVCRT__PC_64; break;
|
|
}
|
|
if (fpword & 0x1000) flags |= MSVCRT__IC_AFFINE;
|
|
|
|
/* Mask with parameters */
|
|
flags = (flags & ~mask) | (newval & mask);
|
|
|
|
/* Convert (masked) value back to fp word */
|
|
fpword = 0;
|
|
if (flags & MSVCRT__EM_INVALID) fpword |= 0x1;
|
|
if (flags & MSVCRT__EM_DENORMAL) fpword |= 0x2;
|
|
if (flags & MSVCRT__EM_ZERODIVIDE) fpword |= 0x4;
|
|
if (flags & MSVCRT__EM_OVERFLOW) fpword |= 0x8;
|
|
if (flags & MSVCRT__EM_UNDERFLOW) fpword |= 0x10;
|
|
if (flags & MSVCRT__EM_INEXACT) fpword |= 0x20;
|
|
switch(flags & (MSVCRT__RC_UP | MSVCRT__RC_DOWN)) {
|
|
case MSVCRT__RC_UP|MSVCRT__RC_DOWN: fpword |= 0xC00; break;
|
|
case MSVCRT__RC_UP: fpword |= 0x800; break;
|
|
case MSVCRT__RC_DOWN: fpword |= 0x400; break;
|
|
}
|
|
switch (flags & (MSVCRT__PC_24 | MSVCRT__PC_53)) {
|
|
case MSVCRT__PC_64: fpword |= 0x300; break;
|
|
case MSVCRT__PC_53: fpword |= 0x200; break;
|
|
case MSVCRT__PC_24: fpword |= 0x0; break;
|
|
}
|
|
if (flags & MSVCRT__IC_AFFINE) fpword |= 0x1000;
|
|
|
|
TRACE("Control word after : %08x\n", fpword);
|
|
|
|
/* Put fp control word */
|
|
__asm__ __volatile__( "fldcw %0" : : "m" (fpword) );
|
|
|
|
return flags;
|
|
#else
|
|
FIXME(":Not Implemented!\n");
|
|
return 0;
|
|
#endif
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _controlfp (MSVCRT.@)
|
|
*/
|
|
unsigned int CDECL _controlfp(unsigned int newval, unsigned int mask)
|
|
{
|
|
#ifdef __i386__
|
|
return _control87( newval, mask & ~MSVCRT__EM_DENORMAL );
|
|
#else
|
|
FIXME(":Not Implemented!\n");
|
|
return 0;
|
|
#endif
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _copysign (MSVCRT.@)
|
|
*/
|
|
double CDECL _copysign(double num, double sign)
|
|
{
|
|
/* FIXME: Behaviour for Nan/Inf? */
|
|
if (sign < 0.0)
|
|
return num < 0.0 ? num : -num;
|
|
return num < 0.0 ? -num : num;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _finite (MSVCRT.@)
|
|
*/
|
|
int CDECL _finite(double num)
|
|
{
|
|
return (finite(num)?1:0); /* See comment for _isnan() */
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _fpreset (MSVCRT.@)
|
|
*/
|
|
void CDECL _fpreset(void)
|
|
{
|
|
#if defined(__GNUC__) && defined(__i386__)
|
|
__asm__ __volatile__( "fninit" );
|
|
#else
|
|
FIXME(":Not Implemented!\n");
|
|
#endif
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _isnan (MSVCRT.@)
|
|
*/
|
|
INT CDECL _isnan(double num)
|
|
{
|
|
/* Some implementations return -1 for true(glibc), msvcrt/crtdll return 1.
|
|
* Do the same, as the result may be used in calculations
|
|
*/
|
|
return isnan(num) ? 1 : 0;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _j0 (MSVCRT.@)
|
|
*/
|
|
double CDECL _j0(double num)
|
|
{
|
|
/* FIXME: errno handling */
|
|
return j0(num);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _j1 (MSVCRT.@)
|
|
*/
|
|
double CDECL _j1(double num)
|
|
{
|
|
/* FIXME: errno handling */
|
|
return j1(num);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* jn (MSVCRT.@)
|
|
*/
|
|
double CDECL _jn(int n, double num)
|
|
{
|
|
/* FIXME: errno handling */
|
|
return jn(n, num);
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _y0 (MSVCRT.@)
|
|
*/
|
|
double CDECL _y0(double num)
|
|
{
|
|
double retval;
|
|
if (!finite(num)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
retval = y0(num);
|
|
if (_fpclass(retval) == MSVCRT__FPCLASS_NINF)
|
|
{
|
|
*MSVCRT__errno() = MSVCRT_EDOM;
|
|
retval = sqrt(-1);
|
|
}
|
|
return retval;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _y1 (MSVCRT.@)
|
|
*/
|
|
double CDECL _y1(double num)
|
|
{
|
|
double retval;
|
|
if (!finite(num)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
retval = y1(num);
|
|
if (_fpclass(retval) == MSVCRT__FPCLASS_NINF)
|
|
{
|
|
*MSVCRT__errno() = MSVCRT_EDOM;
|
|
retval = sqrt(-1);
|
|
}
|
|
return retval;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _yn (MSVCRT.@)
|
|
*/
|
|
double CDECL _yn(int order, double num)
|
|
{
|
|
double retval;
|
|
if (!finite(num)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
retval = yn(order,num);
|
|
if (_fpclass(retval) == MSVCRT__FPCLASS_NINF)
|
|
{
|
|
*MSVCRT__errno() = MSVCRT_EDOM;
|
|
retval = sqrt(-1);
|
|
}
|
|
return retval;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _nextafter (MSVCRT.@)
|
|
*/
|
|
double CDECL _nextafter(double num, double next)
|
|
{
|
|
double retval;
|
|
if (!finite(num) || !finite(next)) *MSVCRT__errno() = MSVCRT_EDOM;
|
|
retval = nextafter(num,next);
|
|
return retval;
|
|
}
|
|
|
|
/*********************************************************************
|
|
* _ecvt (MSVCRT.@)
|
|
*/
|
|
char * CDECL _ecvt( double number, int ndigits, int *decpt, int *sign )
|
|
{
|
|
thread_data_t *data = msvcrt_get_thread_data();
|
|
char *dec;
|
|
|
|
if (!data->efcvt_buffer)
|
|
data->efcvt_buffer = MSVCRT_malloc( 80 ); /* ought to be enough */
|
|
|
|
snprintf(data->efcvt_buffer, 80, "%.*e", ndigits /* FIXME wrong */, number);
|
|
*sign = (number < 0);
|
|
dec = strchr(data->efcvt_buffer, '.');
|
|
*decpt = (dec) ? dec - data->efcvt_buffer : -1;
|
|
return data->efcvt_buffer;
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _fcvt (MSVCRT.@)
|
|
*/
|
|
char * CDECL _fcvt( double number, int ndigits, int *decpt, int *sign )
|
|
{
|
|
thread_data_t *data = msvcrt_get_thread_data();
|
|
int stop, dec1, dec2;
|
|
char *ptr1, *ptr2, *first;
|
|
char buf[80]; /* ought to be enough */
|
|
|
|
if (!data->efcvt_buffer)
|
|
data->efcvt_buffer = MSVCRT_malloc( 80 ); /* ought to be enough */
|
|
|
|
if (number < 0)
|
|
{
|
|
*sign = 1;
|
|
number = -number;
|
|
} else *sign = 0;
|
|
|
|
snprintf(buf, 80, "%.*f", ndigits < 0 ? 0 : ndigits, number);
|
|
ptr1 = buf;
|
|
ptr2 = data->efcvt_buffer;
|
|
first = NULL;
|
|
dec1 = 0;
|
|
dec2 = 0;
|
|
|
|
/* For numbers below the requested resolution, work out where
|
|
the decimal point will be rather than finding it in the string */
|
|
if (number < 1.0 && number > 0.0) {
|
|
dec2 = log10(number + 1e-10);
|
|
if (-dec2 <= ndigits) dec2 = 0;
|
|
}
|
|
|
|
/* If requested digits is zero or less, we will need to truncate
|
|
* the returned string */
|
|
if (ndigits < 1) {
|
|
stop = strlen(buf) + ndigits;
|
|
} else {
|
|
stop = strlen(buf);
|
|
}
|
|
|
|
while (*ptr1 == '0') ptr1++; /* Skip leading zeroes */
|
|
while (*ptr1 != '\0' && *ptr1 != '.') {
|
|
if (!first) first = ptr2;
|
|
if ((ptr1 - buf) < stop) {
|
|
*ptr2++ = *ptr1++;
|
|
} else {
|
|
ptr1++;
|
|
}
|
|
dec1++;
|
|
}
|
|
|
|
if (ndigits > 0) {
|
|
ptr1++;
|
|
if (!first) {
|
|
while (*ptr1 == '0') { /* Process leading zeroes */
|
|
*ptr2++ = *ptr1++;
|
|
dec1--;
|
|
}
|
|
}
|
|
while (*ptr1 != '\0') {
|
|
if (!first) first = ptr2;
|
|
*ptr2++ = *ptr1++;
|
|
}
|
|
}
|
|
|
|
*ptr2 = '\0';
|
|
|
|
/* We never found a non-zero digit, then our number is either
|
|
* smaller than the requested precision, or 0.0 */
|
|
if (!first) {
|
|
if (number > 0.0) {
|
|
first = ptr2;
|
|
} else {
|
|
first = data->efcvt_buffer;
|
|
dec1 = 0;
|
|
}
|
|
}
|
|
|
|
*decpt = dec2 ? dec2 : dec1;
|
|
return first;
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _gcvt (MSVCRT.@)
|
|
*
|
|
* FIXME: uses both E and F.
|
|
*/
|
|
char * CDECL _gcvt( double number, int ndigit, char *buff )
|
|
{
|
|
sprintf(buff, "%.*E", ndigit, number);
|
|
return buff;
|
|
}
|
|
|
|
#include <stdlib.h> /* div_t, ldiv_t */
|
|
|
|
/*********************************************************************
|
|
* div (MSVCRT.@)
|
|
* VERSION
|
|
* [i386] Windows binary compatible - returns the struct in eax/edx.
|
|
*/
|
|
#ifdef __i386__
|
|
unsigned __int64 CDECL MSVCRT_div(int num, int denom)
|
|
{
|
|
div_t dt = div(num,denom);
|
|
return ((unsigned __int64)dt.rem << 32) | (unsigned int)dt.quot;
|
|
}
|
|
#else
|
|
/*********************************************************************
|
|
* div (MSVCRT.@)
|
|
* VERSION
|
|
* [!i386] Non-x86 can't run win32 apps so we don't need binary compatibility
|
|
*/
|
|
MSVCRT_div_t CDECL MSVCRT_div(int num, int denom)
|
|
{
|
|
div_t dt = div(num,denom);
|
|
MSVCRT_div_t ret;
|
|
ret.quot = dt.quot;
|
|
ret.rem = dt.rem;
|
|
|
|
return ret;
|
|
|
|
}
|
|
#endif /* ifdef __i386__ */
|
|
|
|
|
|
/*********************************************************************
|
|
* ldiv (MSVCRT.@)
|
|
* VERSION
|
|
* [i386] Windows binary compatible - returns the struct in eax/edx.
|
|
*/
|
|
#ifdef __i386__
|
|
unsigned __int64 CDECL MSVCRT_ldiv(long num, long denom)
|
|
{
|
|
ldiv_t ldt = ldiv(num,denom);
|
|
return ((unsigned __int64)ldt.rem << 32) | (unsigned long)ldt.quot;
|
|
}
|
|
#else
|
|
/*********************************************************************
|
|
* ldiv (MSVCRT.@)
|
|
* VERSION
|
|
* [!i386] Non-x86 can't run win32 apps so we don't need binary compatibility
|
|
*/
|
|
MSVCRT_ldiv_t CDECL MSVCRT_ldiv(long num, long denom)
|
|
{
|
|
ldiv_t result = ldiv(num,denom);
|
|
|
|
MSVCRT_ldiv_t ret;
|
|
ret.quot = result.quot;
|
|
ret.rem = result.rem;
|
|
|
|
return ret;
|
|
}
|
|
#endif /* ifdef __i386__ */
|
|
|
|
/***********************************************************************
|
|
* _adj_fdiv_m16i (MSVCRT.@)
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void __stdcall _adj_fdiv_m16i( short arg )
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adj_fdiv_m32 (MSVCRT.@)
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void __stdcall _adj_fdiv_m32( unsigned int arg )
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adj_fdiv_m32i (MSVCRT.@)
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void __stdcall _adj_fdiv_m32i( int arg )
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adj_fdiv_m64 (MSVCRT.@)
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void __stdcall _adj_fdiv_m64( unsigned __int64 arg )
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adj_fdiv_r (MSVCRT.@)
|
|
* FIXME
|
|
* This function is likely to have the wrong number of arguments.
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void _adj_fdiv_r(void)
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adj_fdivr_m16i (MSVCRT.@)
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void __stdcall _adj_fdivr_m16i( short arg )
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adj_fdivr_m32 (MSVCRT.@)
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void __stdcall _adj_fdivr_m32( unsigned int arg )
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adj_fdivr_m32i (MSVCRT.@)
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void __stdcall _adj_fdivr_m32i( int arg )
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adj_fdivr_m64 (MSVCRT.@)
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void __stdcall _adj_fdivr_m64( unsigned __int64 arg )
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adj_fpatan (MSVCRT.@)
|
|
* FIXME
|
|
* This function is likely to have the wrong number of arguments.
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void _adj_fpatan(void)
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adj_fprem (MSVCRT.@)
|
|
* FIXME
|
|
* This function is likely to have the wrong number of arguments.
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void _adj_fprem(void)
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adj_fprem1 (MSVCRT.@)
|
|
* FIXME
|
|
* This function is likely to have the wrong number of arguments.
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void _adj_fprem1(void)
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adj_fptan (MSVCRT.@)
|
|
* FIXME
|
|
* This function is likely to have the wrong number of arguments.
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void _adj_fptan(void)
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _adjust_fdiv (MSVCRT.@)
|
|
* FIXME
|
|
* I _think_ this function should be a variable indicating whether
|
|
* Pentium fdiv bug safe code should be used.
|
|
*/
|
|
void _adjust_fdiv(void)
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _safe_fdiv (MSVCRT.@)
|
|
* FIXME
|
|
* This function is likely to have the wrong number of arguments.
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void _safe_fdiv(void)
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _safe_fdivr (MSVCRT.@)
|
|
* FIXME
|
|
* This function is likely to have the wrong number of arguments.
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void _safe_fdivr(void)
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _safe_fprem (MSVCRT.@)
|
|
* FIXME
|
|
* This function is likely to have the wrong number of arguments.
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void _safe_fprem(void)
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|
|
|
|
/***********************************************************************
|
|
* _safe_fprem1 (MSVCRT.@)
|
|
*
|
|
* FIXME
|
|
* This function is likely to have the wrong number of arguments.
|
|
*
|
|
* NOTE
|
|
* I _think_ this function is intended to work around the Pentium
|
|
* fdiv bug.
|
|
*/
|
|
void _safe_fprem1(void)
|
|
{
|
|
TRACE("(): stub\n");
|
|
}
|