\r
/* compute q[0],q[1],...q[jk] */\r
for (i=0;i<=jk;i++) {\r
- for(j=0,fw=0.0;j<=jx;j++) fw += x[j]*f[jx+i-j]; q[i] = fw;\r
+ for(j=0,fw=0.0;j<=jx;j++) fw += x[j]*f[jx+i-j];\r
+ q[i] = fw;\r
}\r
\r
jz = jk;\r
#include "math.h"\r
#include "math_private.h"\r
\r
+#ifndef _IEEE_LIBM\r
static const double\r
o_threshold= 7.09782712893383973096e+02, /* 0x40862E42, 0xFEFA39EF */\r
u_threshold= -7.45133219101941108420e+02; /* 0xc0874910, 0xD52D3051 */\r
+#endif\r
\r
double\r
exp(double x) /* wrapper exp */\r
for (i = 0; i < sp->typecnt; ++i)\r
seen[i] = FALSE;\r
nseen = 0;\r
- for (i = sp->timecnt - 1; i >= 0; --i)\r
- if (!seen[sp->types[i]]) {\r
- seen[sp->types[i]] = TRUE;\r
- types[nseen++] = sp->types[i];\r
- }\r
+ for (i = sp->timecnt - 1; i >= 0; --i)\r
+ if (!seen[sp->types[i]]) {\r
+ seen[sp->types[i]] = TRUE;\r
+ types[nseen++] = sp->types[i];\r
+ }\r
for (sameind = 0; sameind < nseen; ++sameind) {\r
samei = types[sameind];\r
if (sp->ttis[samei].tt_isdst != tmp->tm_isdst)\r