Отправлено TAP 30 января 2004 г. 11:18 В ответ на: А неуж-то AD на double разродились? По IEEE? отправлено SМ 30 января 2004 г. 10:51

у них появилось ETSI (незнаю что ет такое) библиотеки и наверное ieee
сьест не сьест :)
float powf (float x, float y) {
int sign, *p, xind; float v;

if (x == 0.0)
return 0.0;
sign = 0; /* positive sign */
p = (int *) &x; /* point to x */
if (x < 0.0) {
sign = (int) y; /* last bit for sign */
if ( sign != y) {
errno=EDOM;
return 0.0;
}
p[0] &= 0x7fff; /* x is positive now */
}
xind = p[0] >> 7;
xind -= 0x7f; /* convert ieee exponent */
p[0] &= 0x007f;
p[0] |= 0x3f80; /* x is in [1.0, 2.0) */

v = x - 1.0; /* v in [0.0, 1.0) */
v = _log2_1 (v); /* v in [0.0, 1.0) */
v += (double) xind; /* v = log2 (x); */
v *= y; /* v = y * log2 (x) */

y = floor (v); /* y is integer */
xind = (int) y; /* save exponent */
if ( xind < -126 )
return 0.0;
else
if (xind > 127 )
if ( sign&0x1 )
return FLT_MIN;
else
return FLT_MAX;
xind += 0x7f;
xind <<= 7; /* ready for ieee exponent */
v = v - y; /* v is in [0.0, 1.0) */
v = _exp2 (v); /* v is in [1.0, 2.0) */

p = (int *) &v;
p[0] &= 0x007f;
p[0] |= xind; /* v *= 2^exponent */
if (sign & 0x01)
p[0] |= 0x8000; /* negate v if sign is true */
return v;
}

• Я имею в виду 64-битную аппаратно у проца. В С-то можно и не такое придумать — SМ   (30.01.2004 11:21, пустое)
  • про апаратную незнаю Поскольку 218х — TAP   (30.01.2004 13:21, пустое)
    • Ой. Никогда не думал, что кто-то юзает плавучку на не предназначенных для этого камнях... — SМ   (30.01.2004 13:33, пустое)
      • Да ещё и не такое делаем:))вот на 54 в плавучке в узких местах — Бывает   (30.01.2004 16:44, пустое)
      • Если есть процессорное время почему бы не ... — TAP   (30.01.2004 13:51, 60 байт)