reactos/sdk/lib/crt/math/j1_y1.c

35 lines
576 B
C

#include <math.h>
#include <float.h>
#include "ieee754/ieee754.h"
int *_errno(void);
/*
* @unimplemented
*/
double _j1(double num)
{
if (!_finite(num)) *_errno() = EDOM;
return __ieee754_j1(num);
}
/*
* @implemented
*/
double _y1(double num)
{
double retval;
int fpclass = _fpclass(num);
if (!_finite(num) || fpclass == _FPCLASS_NN ||
fpclass == _FPCLASS_ND || fpclass == _FPCLASS_NZ)
*_errno() = EDOM;
retval = __ieee754_y1(num);
if (_fpclass(retval) == _FPCLASS_NINF)
{
*_errno() = EDOM;
retval = sqrt(-1);
}
return retval;
}