iy1
int iy0, iy1, iy2, iy3, iy4;
iy1 = ival1 & 0x007fffff;
ival1 = (iy1 + 0x20000) & 0xfffc0000;
iy1 = iy1 - ival1;
yy1 = (double) iy1 * __TBL_logf[i1 + 1];