Здравствуйте, alpha21264, Вы писали:

A>А для 64-битных double... на целых числах!!!
A>Исходники лежат здесь: (корень)/sysdeps/ieee754/dbl-64
A>Соответсвенно скорость падает в СТО раз.

1. Стоит почитать коммит логи. Возможно там есть объяснение.
Например не исключено что в другом варианте будет теряться точность.

Вообще вычисления математических функций с высокой точностью очень сложная задача с кучей не очевидных моментов.