ReactOS  0.4.14-dev-604-gcfdd483
j0_y0.c File Reference
#include "math.h"
#include "ieee754.h"
Include dependency graph for j0_y0.c:

Go to the source code of this file.

Functions

static double pzero ()
 
static double qzero ()
 
double __ieee754_j0 (double x)
 
double __ieee754_y0 (double x)
 
static double pzero (double x)
 
static double qzero (double x)
 

Variables

static double huge = 1e300
 
static double one = 1.0
 
static double invsqrtpi = 5.64189583547756279280e-01
 
static double tpi = 6.36619772367581382433e-01
 
static double R []
 
static double S []
 
static double zero = 0.0
 
static double U []
 
static double V []
 
static double pR8 [6]
 
static double pS8 [5]
 
static double pR5 [6]
 
static double pS5 [5]
 
static double pR3 [6]
 
static double pS3 [5]
 
static double pR2 [6]
 
static double pS2 [5]
 
static double qR8 [6]
 
static double qS8 [6]
 
static double qR5 [6]
 
static double qS5 [6]
 
static double qR3 [6]
 
static double qS3 [6]
 
static double qR2 [6]
 
static double qS2 [6]
 

Function Documentation

◆ __ieee754_j0()

double __ieee754_j0 ( double  x)

Definition at line 102 of file j0_y0.c.

105 {
106  double z, s,c,ss,cc,r,u,v,r1,r2,s1,s2,z2,z4;
107  int32_t hx,ix;
108 
109  GET_HIGH_WORD(hx,x);
110  ix = hx&0x7fffffff;
111  if(ix>=0x7ff00000) return one/(x*x);
112  x = fabs(x);
113  if(ix >= 0x40000000) { /* |x| >= 2.0 */
114  __sincos (x, &s, &c);
115  ss = s-c;
116  cc = s+c;
117  if(ix<0x7fe00000) { /* make sure x+x not overflow */
118  z = -__cos(x+x);
119  if ((s*c)<zero) cc = z/ss;
120  else ss = z/cc;
121  }
122  /*
123  * j0(x) = 1/sqrt(pi) * (P(0,x)*cc - Q(0,x)*ss) / sqrt(x)
124  * y0(x) = 1/sqrt(pi) * (P(0,x)*ss + Q(0,x)*cc) / sqrt(x)
125  */
126  if(ix>0x48000000) z = (invsqrtpi*cc)/__ieee754_sqrt(x);
127  else {
128  u = pzero(x); v = qzero(x);
130  }
131  return z;
132  }
133  if(ix<0x3f200000) { /* |x| < 2**-13 */
134  if(huge+x>one) { /* raise inexact if x != 0 */
135  if(ix<0x3e400000) return one; /* |x|<2**-27 */
136  else return one - 0.25*x*x;
137  }
138  }
139  z = x*x;
140 #ifdef DO_NOT_USE_THIS
141  r = z*(R02+z*(R03+z*(R04+z*R05)));
142  s = one+z*(S01+z*(S02+z*(S03+z*S04)));
143 #else
144  r1 = z*R[2]; z2=z*z;
145  r2 = R[3]+z*R[4]; z4=z2*z2;
146  r = r1 + z2*r2 + z4*R[5];
147  s1 = one+z*S[1];
148  s2 = S[2]+z*S[3];
149  s = s1 + z2*s2 + z4*S[4];
150 #endif
151  if(ix < 0x3FF00000) { /* |x| < 1.00 */
152  return one + z*(-0.25+(r/s));
153  } else {
154  u = 0.5*x;
155  return((one+u)*(one-u)+z*(r/s));
156  }
157 }
GLsizei GLenum const GLvoid GLsizei GLenum GLbyte GLbyte GLbyte GLdouble GLdouble GLdouble GLfloat GLfloat GLfloat GLint GLint GLint GLshort GLshort GLshort GLubyte GLubyte GLubyte GLuint GLuint GLuint GLushort GLushort GLushort GLbyte GLbyte GLbyte GLbyte GLdouble GLdouble GLdouble GLdouble GLfloat GLfloat GLfloat GLfloat GLint GLint GLint GLint GLshort GLshort GLshort GLshort GLubyte GLubyte GLubyte GLubyte GLuint GLuint GLuint GLuint GLushort GLushort GLushort GLushort GLboolean const GLdouble const GLfloat const GLint const GLshort const GLbyte const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLdouble const GLfloat const GLfloat const GLint const GLint const GLshort const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort GLenum GLenum GLenum GLfloat GLenum GLint GLenum GLenum GLenum GLfloat GLenum GLenum GLint GLenum GLfloat GLenum GLint GLint GLushort GLenum GLenum GLfloat GLenum GLenum GLint GLfloat const GLubyte GLenum GLenum GLenum const GLfloat GLenum GLenum const GLint GLenum GLint GLint GLsizei GLsizei GLint GLenum GLenum const GLvoid GLenum GLenum const GLfloat GLenum GLenum const GLint GLenum GLenum const GLdouble GLenum GLenum const GLfloat GLenum GLenum const GLint GLsizei GLuint GLfloat GLuint GLbitfield GLfloat GLint GLuint GLboolean GLenum GLfloat GLenum GLbitfield GLenum GLfloat GLfloat GLint GLint const GLfloat GLenum GLfloat GLfloat GLint GLint GLfloat GLfloat GLint GLint const GLfloat GLint GLfloat GLfloat GLint GLfloat GLfloat GLint GLfloat GLfloat const GLdouble * u
Definition: glfuncs.h:240
static double huge
Definition: j0_y0.c:79
struct S2 s2
GLdouble GLdouble GLdouble r
Definition: gl.h:2055
static __inline double __ieee754_sqrt(double x)
Definition: ieee754.h:40
static DNS_RECORDW r1
Definition: record.c:37
GLint GLint GLint GLint GLint x
Definition: gl.h:1548
Definition: movable.cpp:7
static __inline double __cos(double x)
Definition: ieee754.h:42
static double one
Definition: j0_y0.c:80
#define GET_HIGH_WORD(i, d)
Definition: ieee754.h:26
static double pzero()
GLdouble GLdouble z
Definition: glext.h:5874
static double qzero()
struct S1 s1
static DNS_RECORDW r2
Definition: record.c:38
const GLubyte * c
Definition: glext.h:8905
static double zero
Definition: j0_y0.c:96
GLdouble s
Definition: gl.h:2039
static double invsqrtpi
Definition: j0_y0.c:81
INT32 int32_t
Definition: types.h:71
_Check_return_ _CRT_JIT_INTRINSIC double __cdecl fabs(_In_ double x)
uint32_t cc
Definition: isohybrid.c:75
static __inline void __sincos(double x, double *s, double *c)
Definition: ieee754.h:43
const GLdouble * v
Definition: gl.h:2040
static double R[]
Definition: j0_y0.c:84
#define c
Definition: ke_i.h:80
#define ss
Definition: i386-dis.c:432

Referenced by __ieee754_jn(), __ieee754_y0(), and _j0().

◆ __ieee754_y0()

double __ieee754_y0 ( double  x)

Definition at line 179 of file j0_y0.c.

182 {
183  double z, s,c,ss,cc,u,v,z2,z4,z6,u1,u2,u3,v1,v2;
184  int32_t hx,ix,lx;
185 
186  EXTRACT_WORDS(hx,lx,x);
187  ix = 0x7fffffff&hx;
188  /* Y0(NaN) is NaN, y0(-inf) is Nan, y0(inf) is 0, y0(0) is -inf. */
189  if(ix>=0x7ff00000) return one/(x+x*x);
190  if((ix|lx)==0) return -HUGE_VAL+x; /* -inf and overflow exception. */
191  if(hx<0) return zero/(zero*x);
192  if(ix >= 0x40000000) { /* |x| >= 2.0 */
193  /* y0(x) = sqrt(2/(pi*x))*(p0(x)*sin(x0)+q0(x)*cos(x0))
194  * where x0 = x-pi/4
195  * Better formula:
196  * cos(x0) = cos(x)cos(pi/4)+sin(x)sin(pi/4)
197  * = 1/sqrt(2) * (sin(x) + cos(x))
198  * sin(x0) = sin(x)cos(3pi/4)-cos(x)sin(3pi/4)
199  * = 1/sqrt(2) * (sin(x) - cos(x))
200  * To avoid cancellation, use
201  * sin(x) +- cos(x) = -cos(2x)/(sin(x) -+ cos(x))
202  * to compute the worse one.
203  */
204  __sincos (x, &s, &c);
205  ss = s-c;
206  cc = s+c;
207  /*
208  * j0(x) = 1/sqrt(pi) * (P(0,x)*cc - Q(0,x)*ss) / sqrt(x)
209  * y0(x) = 1/sqrt(pi) * (P(0,x)*ss + Q(0,x)*cc) / sqrt(x)
210  */
211  if(ix<0x7fe00000) { /* make sure x+x not overflow */
212  z = -__cos(x+x);
213  if ((s*c)<zero) cc = z/ss;
214  else ss = z/cc;
215  }
216  if(ix>0x48000000) z = (invsqrtpi*ss)/__ieee754_sqrt(x);
217  else {
218  u = pzero(x); v = qzero(x);
220  }
221  return z;
222  }
223  if(ix<=0x3e400000) { /* x < 2**-27 */
224  return(U[0] + tpi*__ieee754_log(x));
225  }
226  z = x*x;
227 #ifdef DO_NOT_USE_THIS
228  u = u00+z*(u01+z*(u02+z*(u03+z*(u04+z*(u05+z*u06)))));
229  v = one+z*(v01+z*(v02+z*(v03+z*v04)));
230 #else
231  u1 = U[0]+z*U[1]; z2=z*z;
232  u2 = U[2]+z*U[3]; z4=z2*z2;
233  u3 = U[4]+z*U[5]; z6=z4*z2;
234  u = u1 + z2*u2 + z4*u3 + z6*U[6];
235  v1 = one+z*V[0];
236  v2 = V[1]+z*V[2];
237  v = v1 + z2*v2 + z4*V[3];
238 #endif
239  return(u/v + tpi*(__ieee754_j0(x)*__ieee754_log(x)));
240 }
GLsizei GLenum const GLvoid GLsizei GLenum GLbyte GLbyte GLbyte GLdouble GLdouble GLdouble GLfloat GLfloat GLfloat GLint GLint GLint GLshort GLshort GLshort GLubyte GLubyte GLubyte GLuint GLuint GLuint GLushort GLushort GLushort GLbyte GLbyte GLbyte GLbyte GLdouble GLdouble GLdouble GLdouble GLfloat GLfloat GLfloat GLfloat GLint GLint GLint GLint GLshort GLshort GLshort GLshort GLubyte GLubyte GLubyte GLubyte GLuint GLuint GLuint GLuint GLushort GLushort GLushort GLushort GLboolean const GLdouble const GLfloat const GLint const GLshort const GLbyte const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLdouble const GLfloat const GLfloat const GLint const GLint const GLshort const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort const GLdouble const GLfloat const GLint const GLshort GLenum GLenum GLenum GLfloat GLenum GLint GLenum GLenum GLenum GLfloat GLenum GLenum GLint GLenum GLfloat GLenum GLint GLint GLushort GLenum GLenum GLfloat GLenum GLenum GLint GLfloat const GLubyte GLenum GLenum GLenum const GLfloat GLenum GLenum const GLint GLenum GLint GLint GLsizei GLsizei GLint GLenum GLenum const GLvoid GLenum GLenum const GLfloat GLenum GLenum const GLint GLenum GLenum const GLdouble GLenum GLenum const GLfloat GLenum GLenum const GLint GLsizei GLuint GLfloat GLuint GLbitfield GLfloat GLint GLuint GLboolean GLenum GLfloat GLenum GLbitfield GLenum GLfloat GLfloat GLint GLint const GLfloat GLenum GLfloat GLfloat GLint GLint GLfloat GLfloat GLint GLint const GLfloat GLint GLfloat GLfloat GLint GLfloat GLfloat GLint GLfloat GLfloat const GLdouble * u
Definition: glfuncs.h:240
GLdouble GLdouble u2
Definition: glext.h:8308
static __inline double __ieee754_sqrt(double x)
Definition: ieee754.h:40
GLint GLint GLint GLint GLint x
Definition: gl.h:1548
double __ieee754_j0(double x)
Definition: j0_y0.c:102
GLdouble u1
Definition: glext.h:8308
static __inline double __cos(double x)
Definition: ieee754.h:42
static double one
Definition: j0_y0.c:80
static double pzero()
GLdouble GLdouble z
Definition: glext.h:5874
static double qzero()
#define EXTRACT_WORDS(ix0, ix1, d)
Definition: ieee754.h:16
const GLubyte * c
Definition: glext.h:8905
static double zero
Definition: j0_y0.c:96
GLdouble s
Definition: gl.h:2039
static double invsqrtpi
Definition: j0_y0.c:81
INT32 int32_t
Definition: types.h:71
GLfloat GLfloat GLfloat v2
Definition: glext.h:6063
uint32_t cc
Definition: isohybrid.c:75
static __inline void __sincos(double x, double *s, double *c)
Definition: ieee754.h:43
static double tpi
Definition: j0_y0.c:82
const GLdouble * v
Definition: gl.h:2040
#define HUGE_VAL
Definition: math.h:51
#define c
Definition: ke_i.h:80
static __inline double __ieee754_log(double x)
Definition: ieee754.h:41
static BYTE u3[]
Definition: msg.c:580
#define ss
Definition: i386-dis.c:432
GLfloat GLfloat v1
Definition: glext.h:6062
static double U[]
Definition: j0_y0.c:164

Referenced by __ieee754_yn(), and _y0().

◆ pzero() [1/2]

static double pzero ( )
static

Referenced by __ieee754_j0(), and __ieee754_y0().

◆ pzero() [2/2]

static double pzero ( double  x)
static

Definition at line 350 of file j0_y0.c.

353 {
354 #ifdef __STDC__
355  const double *p = 0,*q = 0;
356 #else
357  double *p = 0,*q = 0;
358 #endif
359  double z,r,s,z2,z4,r1,r2,r3,s1,s2,s3;
360  int32_t ix;
361  GET_HIGH_WORD(ix,x);
362  ix &= 0x7fffffff;
363  if(ix>=0x40200000) {p = pR8; q= pS8;}
364  else if(ix>=0x40122E8B){p = pR5; q= pS5;}
365  else if(ix>=0x4006DB6D){p = pR3; q= pS3;}
366  else if(ix>=0x40000000){p = pR2; q= pS2;}
367  z = one/(x*x);
368 #ifdef DO_NOT_USE_THIS
369  r = p[0]+z*(p[1]+z*(p[2]+z*(p[3]+z*(p[4]+z*p[5]))));
370  s = one+z*(q[0]+z*(q[1]+z*(q[2]+z*(q[3]+z*q[4]))));
371 #else
372  r1 = p[0]+z*p[1]; z2=z*z;
373  r2 = p[2]+z*p[3]; z4=z2*z2;
374  r3 = p[4]+z*p[5];
375  r = r1 + z2*r2 + z4*r3;
376  s1 = one+z*q[0];
377  s2 = q[1]+z*q[2];
378  s3 = q[3]+z*q[4];
379  s = s1 + z2*s2 + z4*s3;
380 #endif
381  return one+ r/s;
382 }
struct S2 s2
static DNS_RECORDW r3
Definition: record.c:39
static double pR3[6]
Definition: j0_y0.c:302
GLdouble GLdouble GLdouble r
Definition: gl.h:2055
static DNS_RECORDW r1
Definition: record.c:37
static double pR2[6]
Definition: j0_y0.c:326
GLint GLint GLint GLint GLint x
Definition: gl.h:1548
static double pR5[6]
Definition: j0_y0.c:278
static double one
Definition: j0_y0.c:80
#define GET_HIGH_WORD(i, d)
Definition: ieee754.h:26
static double pS3[5]
Definition: j0_y0.c:314
GLdouble GLdouble z
Definition: glext.h:5874
struct S1 s1
static DNS_RECORDW r2
Definition: record.c:38
GLdouble GLdouble GLdouble GLdouble q
Definition: gl.h:2063
static double pS5[5]
Definition: j0_y0.c:290
GLdouble s
Definition: gl.h:2039
INT32 int32_t
Definition: types.h:71
static double pR8[6]
Definition: j0_y0.c:254
static double pS2[5]
Definition: j0_y0.c:338
static double pS8[5]
Definition: j0_y0.c:266
GLfloat GLfloat p
Definition: glext.h:8902

◆ qzero() [1/2]

static double qzero ( )

Referenced by __ieee754_j0(), and __ieee754_y0().

◆ qzero() [2/2]

static double qzero ( double  x)
static

Definition at line 497 of file j0_y0.c.

500 {
501 #ifdef __STDC__
502  const double *p = 0,*q = 0;
503 #else
504  double *p = 0,*q = 0;
505 #endif
506  double s,r,z,z2,z4,z6,r1,r2,r3,s1,s2,s3;
507  int32_t ix;
508  GET_HIGH_WORD(ix,x);
509  ix &= 0x7fffffff;
510  if(ix>=0x40200000) {p = qR8; q= qS8;}
511  else if(ix>=0x40122E8B){p = qR5; q= qS5;}
512  else if(ix>=0x4006DB6D){p = qR3; q= qS3;}
513  else if(ix>=0x40000000){p = qR2; q= qS2;}
514  z = one/(x*x);
515 #ifdef DO_NOT_USE_THIS
516  r = p[0]+z*(p[1]+z*(p[2]+z*(p[3]+z*(p[4]+z*p[5]))));
517  s = one+z*(q[0]+z*(q[1]+z*(q[2]+z*(q[3]+z*(q[4]+z*q[5])))));
518 #else
519  r1 = p[0]+z*p[1]; z2=z*z;
520  r2 = p[2]+z*p[3]; z4=z2*z2;
521  r3 = p[4]+z*p[5]; z6=z4*z2;
522  r= r1 + z2*r2 + z4*r3;
523  s1 = one+z*q[0];
524  s2 = q[1]+z*q[2];
525  s3 = q[3]+z*q[4];
526  s = s1 + z2*s2 + z4*s3 +z6*q[5];
527 #endif
528  return (-.125 + r/s)/x;
529 }
struct S2 s2
static DNS_RECORDW r3
Definition: record.c:39
GLdouble GLdouble GLdouble r
Definition: gl.h:2055
static double qS5[6]
Definition: j0_y0.c:434
static double qR5[6]
Definition: j0_y0.c:422
static DNS_RECORDW r1
Definition: record.c:37
GLint GLint GLint GLint GLint x
Definition: gl.h:1548
static double one
Definition: j0_y0.c:80
#define GET_HIGH_WORD(i, d)
Definition: ieee754.h:26
GLdouble GLdouble z
Definition: glext.h:5874
static double qS2[6]
Definition: j0_y0.c:484
static double qR8[6]
Definition: j0_y0.c:397
struct S1 s1
static DNS_RECORDW r2
Definition: record.c:38
static double qR3[6]
Definition: j0_y0.c:447
GLdouble GLdouble GLdouble GLdouble q
Definition: gl.h:2063
GLdouble s
Definition: gl.h:2039
static double qS3[6]
Definition: j0_y0.c:459
INT32 int32_t
Definition: types.h:71
static double qS8[6]
Definition: j0_y0.c:409
static double qR2[6]
Definition: j0_y0.c:472
GLfloat GLfloat p
Definition: glext.h:8902

Variable Documentation

◆ huge

double huge = 1e300
static

Definition at line 79 of file j0_y0.c.

Referenced by __ieee754_j0().

◆ invsqrtpi

double invsqrtpi = 5.64189583547756279280e-01
static

Definition at line 81 of file j0_y0.c.

Referenced by __ieee754_j0(), and __ieee754_y0().

◆ one

◆ pR2

double pR2[6]
static
Initial value:
= {
-8.87534333032526411254e-08,
-7.03030995483624743247e-02,
-1.45073846780952986357e+00,
-7.63569613823527770791e+00,
-1.11931668860356747786e+01,
-3.23364579351335335033e+00,
}

Definition at line 326 of file j0_y0.c.

Referenced by pzero().

◆ pR3

double pR3[6]
static
Initial value:
= {
-2.54704601771951915620e-09,
-7.03119616381481654654e-02,
-2.40903221549529611423e+00,
-2.19659774734883086467e+01,
-5.80791704701737572236e+01,
-3.14479470594888503854e+01,
}

Definition at line 302 of file j0_y0.c.

Referenced by pzero().

◆ pR5

double pR5[6]
static
Initial value:
= {
-1.14125464691894502584e-11,
-7.03124940873599280078e-02,
-4.15961064470587782438e+00,
-6.76747652265167261021e+01,
-3.31231299649172967747e+02,
-3.46433388365604912451e+02,
}

Definition at line 278 of file j0_y0.c.

Referenced by pzero().

◆ pR8

double pR8[6]
static
Initial value:
= {
0.00000000000000000000e+00,
-7.03124999999900357484e-02,
-8.08167041275349795626e+00,
-2.57063105679704847262e+02,
-2.48521641009428822144e+03,
-5.25304380490729545272e+03,
}

Definition at line 254 of file j0_y0.c.

Referenced by pzero().

◆ pS2

double pS2[5]
static
Initial value:
= {
2.22202997532088808441e+01,
1.36206794218215208048e+02,
2.70470278658083486789e+02,
1.53875394208320329881e+02,
1.46576176948256193810e+01,
}

Definition at line 338 of file j0_y0.c.

Referenced by pzero().

◆ pS3

double pS3[5]
static
Initial value:
= {
3.58560338055209726349e+01,
3.61513983050303863820e+02,
1.19360783792111533330e+03,
1.12799679856907414432e+03,
1.73580930813335754692e+02,
}

Definition at line 314 of file j0_y0.c.

Referenced by pzero().

◆ pS5

double pS5[5]
static
Initial value:
= {
6.07539382692300335975e+01,
1.05125230595704579173e+03,
5.97897094333855784498e+03,
9.62544514357774460223e+03,
2.40605815922939109441e+03,
}

Definition at line 290 of file j0_y0.c.

Referenced by pzero().

◆ pS8

double pS8[5]
static
Initial value:
= {
1.16534364619668181717e+02,
3.83374475364121826715e+03,
4.05978572648472545552e+04,
1.16752972564375915681e+05,
4.76277284146730962675e+04,
}

Definition at line 266 of file j0_y0.c.

Referenced by pzero().

◆ qR2

double qR2[6]
static
Initial value:
= {
1.50444444886983272379e-07,
7.32234265963079278272e-02,
1.99819174093815998816e+00,
1.44956029347885735348e+01,
3.16662317504781540833e+01,
1.62527075710929267416e+01,
}

Definition at line 472 of file j0_y0.c.

Referenced by qzero().

◆ qR3

double qR3[6]
static
Initial value:
= {
4.37741014089738620906e-09,
7.32411180042911447163e-02,
3.34423137516170720929e+00,
4.26218440745412650017e+01,
1.70808091340565596283e+02,
1.66733948696651168575e+02,
}

Definition at line 447 of file j0_y0.c.

Referenced by qzero().

◆ qR5

double qR5[6]
static
Initial value:
= {
1.84085963594515531381e-11,
7.32421766612684765896e-02,
5.83563508962056953777e+00,
1.35111577286449829671e+02,
1.02724376596164097464e+03,
1.98997785864605384631e+03,
}

Definition at line 422 of file j0_y0.c.

Referenced by qzero().

◆ qR8

double qR8[6]
static
Initial value:
= {
0.00000000000000000000e+00,
7.32421874999935051953e-02,
1.17682064682252693899e+01,
5.57673380256401856059e+02,
8.85919720756468632317e+03,
3.70146267776887834771e+04,
}

Definition at line 397 of file j0_y0.c.

Referenced by qzero().

◆ qS2

double qS2[6]
static
Initial value:
= {
3.03655848355219184498e+01,
2.69348118608049844624e+02,
8.44783757595320139444e+02,
8.82935845112488550512e+02,
2.12666388511798828631e+02,
-5.31095493882666946917e+00,
}

Definition at line 484 of file j0_y0.c.

Referenced by qzero().

◆ qS3

double qS3[6]
static
Initial value:
= {
4.87588729724587182091e+01,
7.09689221056606015736e+02,
3.70414822620111362994e+03,
6.46042516752568917582e+03,
2.51633368920368957333e+03,
-1.49247451836156386662e+02,
}

Definition at line 459 of file j0_y0.c.

Referenced by qzero().

◆ qS5

double qS5[6]
static
Initial value:
= {
8.27766102236537761883e+01,
2.07781416421392987104e+03,
1.88472887785718085070e+04,
5.67511122894947329769e+04,
3.59767538425114471465e+04,
-5.35434275601944773371e+03,
}

Definition at line 434 of file j0_y0.c.

Referenced by qzero().

◆ qS8

double qS8[6]
static
Initial value:
= {
1.63776026895689824414e+02,
8.09834494656449805916e+03,
1.42538291419120476348e+05,
8.03309257119514397345e+05,
8.40501579819060512818e+05,
-3.43899293537866615225e+05,
}

Definition at line 409 of file j0_y0.c.

Referenced by qzero().

◆ R

double R[]
static
Initial value:
= {0.0, 0.0, 1.56249999999999947958e-02,
-1.89979294238854721751e-04,
1.82954049532700665670e-06,
-4.61832688532103189199e-09}

Definition at line 84 of file j0_y0.c.

Referenced by __ieee754_j0().

◆ S

double S[]
static
Initial value:
= {0.0, 1.56191029464890010492e-02,
1.16926784663337450260e-04,
5.13546550207318111446e-07,
1.16614003333790000205e-09}

Definition at line 88 of file j0_y0.c.

◆ tpi

double tpi = 6.36619772367581382433e-01
static

Definition at line 82 of file j0_y0.c.

Referenced by __ieee754_y0().

◆ U

double U[]
static
Initial value:
= {-7.38042951086872317523e-02,
1.76666452509181115538e-01,
-1.38185671945596898896e-02,
3.47453432093683650238e-04,
-3.81407053724364161125e-06,
1.95590137035022920206e-08,
-3.98205194132103398453e-11}

Definition at line 164 of file j0_y0.c.

Referenced by __ieee754_y0().

◆ V

double V[]
static
Initial value:
= {1.27304834834123699328e-02,
7.60068627350353253702e-05,
2.59150851840457805467e-07,
4.41110311332675467403e-10}

Definition at line 171 of file j0_y0.c.

◆ zero

double zero = 0.0
static

Definition at line 96 of file j0_y0.c.

Referenced by __ieee754_j0(), __ieee754_y0(), ACTION_FileVersionMatches(), add_userkey(), apartment_release(), AtapiRegCheckParameterValue(), bsc_thread(), CdRomCheckRegistryForMediaChangeValue(), ClassMoniker_Load(), ClassMoniker_Save(), ClasspSetMediaChangeStateEx(), ClassQueryTimeOutRegistryValue(), clear_first_megabyte(), ComponentFactory_CreateMetadataReaderFromContainer(), cowait_unmarshal_thread(), create_body_offset_list(), create_decoder(), CryptSIPRetrieveSubjectGuid(), D3DXQuaternionSquadSetup(), do_file_copyW(), encode_backchannel_sec_parms(), encode_op_bind_conn_to_session(), encode_op_close(), encode_op_delegpurge(), encode_op_exchange_id(), encode_op_layoutcommit(), encode_op_layoutreturn(), encode_op_reclaim_complete(), EngAllocMem(), Ext2SetReparsePoint(), ftruncate_growable(), g_malloc(), get_block_bh(), get_block_bh_mdl(), get_block_bh_pin(), get_stream_size(), HTTP_ParseDate(), InputProcessorProfiles_AddLanguageProfile(), interp_step(), IPADDRESS_UpdateText(), IStream_Reset(), marshal_stgmed(), marshal_WdtpInterfacePointer(), MCIAVI_mciInfo(), METAFILE_CreateCompressedImageStream(), MIDI_mciInfo(), MimeMessage_Load(), msi_init_string_table(), nfs41_send_compound(), nfstime_abs(), open_file2(), pshader_hw_texdepth(), pshader_hw_texm3x2depth(), ReadBytes(), RegisterDragDrop(), ScsiClassQueryTimeOutRegistryValue(), shader_hw_mov(), shader_hw_nrm(), shader_hw_pow(), SHIsEmptyStream(), StdGlobalInterfaceTable_RegisterInterfaceInGlobal(), StorageBaseImpl_DeleteStorageTree(), stream_get_chunk(), test_bsc_marshaling(), test_decode_msg_get_param(), test_ExitCode(), test_external_connection(), test_GetAddrInfoW(), test_GetDataHere_Proxy(), test_IStream_invalid_operations(), test_marshal_dispinterface(), test_MessageOptions(), test_send(), TransactedSnapshotImpl_Commit(), TransactedSnapshotImpl_DestroyTemporaryCopy(), TransactedSnapshotImpl_Revert(), TransactedSnapshotImpl_StreamSetSize(), trio_nzero(), unmarshal_stgmed(), VfdLoadLink(), vshader_add_footer(), WAVE_mciInfo(), WCachePreReadPacket__(), WCacheUpdatePacket(), xdr_channel_attrs4(), and xslprocessor_transform().