Потеря точности при преобразовании типов?

Question

[isvirin]

Вот такая функция:

void fp()
{
    int64_t tmpn=(((1<<30)-1)/3*2)*1000LL;
    double tmpd=tmpn/1000.0;
    int64_t tmpm=(int64_t)(tmpd*1000.0);
    printf("%lld %.3f %lld\n", tmpn, tmpd, tmpm);
    fflush(stdout);
}

Вызываю ее в начале и в конце функции main(). Выхлоп отличается. В начале:
715827882000 715827882.000 715827882000

В конце:
715827882000 715827904.000 715827904512


Если внимательно посмотреть на число 715827904.000, можно увидеть, что в мантиссе осталось 24 бита, остальные — нули. Т.е. как будто целое число преобразовывали не в double, а во float.


Проявляется это только на винде, в 32-битной mingw. В 64-битных линухах всё ок.


А теперь два вопроса, которые поднимал еще Чернышевский:

Кто виноват?

Что делать?

Answer 1

[grechnik]

Похоже, что кто-то где-то внутри функций, вызываемых из main, вызвал функцию _controlfp (или машинную инструкцию fldcw) и не вернул состояние назад. Как решение, можно вызвать _controlfp(_PC_64, _MCW_PC).

Comments

[isvirin]

Спасибо. Это, действительно помогло! Правда непонятно, как нам это знание теперь использовать:) У нас достаточно большой проект, использующий целый ряд внешних библиотек, в т.ч. кодеков (на них и грешим). Мы сами явно _controlfp не вызываем, а значит это делает кто-то. Непонятно, в какой момент времени нам делать вызов _controlfp(_PC_64, _MCW_PC), чтобы не слишком рано и не слишком поздно:)

[grechnik]

Может, создать отдельный поток для выполнения части кода? Точность FPU, если я не ошибаюсь, индивидуальна для каждого потока.

[isvirin]

У нас и так под сотню потоков. Двигаясь по этому пути остается только изолировать каждую внешнюю библиотеку в отдельном потоке, но это весьма непростой вариант. Мне кажется, что это некорректное поведение библиотеки «вызвать» и не вернуть обратно. Правильный путь — найти, какая библиотека «гадит» и избавиться от нее (заменить).

[grechnik]

Найти, наверное, проще всего экспериментально, вставляя проверки состояния в стратегически важных местах. Обнаружить изменение (без созерцания выхлопов printf) можно той же _controlfp, (_controlfp(0, 0) & _MCW_PC) == _PC_64 или _PC_53.