136 1.57079625129699707031e+00,
137 7.54978941586159635335e-08,
138 5.39030252995776476554e-15,
139 3.28200341580791294123e-22,
140 1.27065575308067607349e-29,
141 1.22933308981111328932e-36,
142 2.73370053816464559624e-44,
143 2.16741683877804819444e-51,
154 int32_t jz,jx,jv,jp,jk,carry,
n,iq[20],
i,
j,
k,
m,q0,ih;
155 double z,fw,
f[20],fq[20],
q[20];
169 jv = (e0-3)/24;
if(jv<0) jv=0;
173 j = jv-jx;
m = jx+jk;
174 for(
i=0;
i<=
m;
i++,
j++)
f[
i] = (
j<0)?
zero : (double) ipio2[
j];
180 for (
i=0;
i<=jk;
i++) {
181 for(
j=0,fw=0.0;
j<=jx;
j++) fw +=
x[
j]*
f[jx+
i-
j];
188 for(
i=0,
j=jz,
z=
q[jz];
j>0;
i++,
j--) {
194 SDL_memset(&iq[jz], 0,
sizeof (
q) - (jz *
sizeof (iq[0])));
204 i = (iq[jz-1]>>(24-q0));
n +=
i;
205 iq[jz-1] -=
i<<(24-q0);
206 ih = iq[jz-1]>>(23-q0);
208 else if(q0==0) ih = iq[jz-1]>>23;
209 else if(
z>=0.5) ih=2;
217 carry = 1; iq[
i] = 0x1000000-
j;
219 }
else iq[
i] = 0xffffff -
j;
224 iq[jz-1] &= 0x7fffff;
break;
226 iq[jz-1] &= 0x3fffff;
break;
238 for (
i=jz-1;
i>=jk;
i--)
j |= iq[
i];
240 for(
k=1;iq[jk-
k]==0;
k++);
242 for(
i=jz+1;
i<=jz+
k;
i++) {
243 f[jx+
i] = (double) ipio2[jv+
i];
244 for(
j=0,fw=0.0;
j<=jx;
j++) fw +=
x[
j]*
f[jx+
i-
j];
256 while(iq[jz]==0) { jz--;
SDL_assert(jz >= 0); q0-=24;}
279 SDL_memset(&fq[jz+1], 0,
sizeof (fq) - ((jz+1) *
sizeof (fq[0])));
286 for (
i=jz;
i>=0;
i--) fw += fq[
i];
287 y[0] = (ih==0)? fw: -fw;
292 for (
i=jz;
i>=0;
i--) fw += fq[
i];
293 y[0] = (ih==0)? fw: -fw;
295 for (
i=1;
i<=jz;
i++) fw += fq[
i];
296 y[1] = (ih==0)? fw: -fw;
309 for (fw=0.0,
i=jz;
i>=2;
i--) fw += fq[
i];
311 y[0] = fq[0];
y[1] = fq[1];
y[2] = fw;
313 y[0] = -fq[0];
y[1] = -fq[1];
y[2] = -fw;