ReactOS 0.4.16-dev-2104-gb84fa49
j1_y1.c
Go to the documentation of this file.
1#include <math.h>
2#include <float.h>
3#include "ieee754/ieee754.h"
4
5int *_errno(void);
6
7/*
8 * @unimplemented
9 */
10double _j1(double num)
11{
12 if (!_finite(num)) *_errno() = EDOM;
13 return __ieee754_j1(num);
14}
15
16/*
17 * @implemented
18 */
19double _y1(double num)
20{
21 double retval;
22 int fpclass = _fpclass(num);
23
24 if (!_finite(num) || fpclass == _FPCLASS_NN ||
25 fpclass == _FPCLASS_ND || fpclass == _FPCLASS_NZ)
26 *_errno() = EDOM;
27
28 retval = __ieee754_y1(num);
29 if (_fpclass(retval) == _FPCLASS_NINF)
30 {
31 *_errno() = EDOM;
32 retval = sqrt(-1);
33 }
34 return retval;
35}
EDOM
#define EDOM
Definition: errno.h:54
_finite
_ACRTIMP int __cdecl _finite(double)
Definition: math.c:1582
_FPCLASS_ND
#define _FPCLASS_ND
Definition: float.h:111
_FPCLASS_NINF
#define _FPCLASS_NINF
Definition: float.h:109
_FPCLASS_NN
#define _FPCLASS_NN
Definition: float.h:110
_fpclass
_ACRTIMP int __cdecl _fpclass(double)
Definition: math.c:595
_FPCLASS_NZ
#define _FPCLASS_NZ
Definition: float.h:112
sqrt
_ACRTIMP double __cdecl sqrt(double)
Definition: sqrt.c:5
num
GLuint GLuint num
Definition: glext.h:9618
__ieee754_y1
double __ieee754_y1(double x)
Definition: j1_y1.c:182
__ieee754_j1
double __ieee754_j1(double x)
Definition: j1_y1.c:103
_y1
double _y1(double num)
Definition: j1_y1.c:19
_errno
int * _errno(void)
Definition: errno.c:215
_j1
double _j1(double num)
Definition: j1_y1.c:10
retval
int retval
Definition: wcstombs.cpp:91