Подразделение Ньютона-Рафсона для деления с плавающей точкой?

Question

Я пытаюсь реализовать алгоритм деления Ньютона-Рафсона Запись в Википедии , чтобы реализовать 32-разрядное деление IEEE-754 с плавающей запятой на процессоре, который не имеет аппаратного делителя.

В моей памяти хранятся 32-битные слова дополнения до двух, и я уже реализовал сложение, вычитание и умножение с плавающей запятой, поэтому я могу повторно использовать код для реализации алгоритма Ньютона-Рафсона. Я пытаюсь сначала реализовать все это в Matlab.

На этом этапе: X_0 = 48/17 - 32/17 * D
Как правильно сместить бит D в диапазоне от 0,5 до 1, как описано в деталях алгоритма?

Answer 1

Вы можете посмотреть библиотеку времени выполнения compiler-rt (часть LLVM), которая имеет либеральную лицензию и реализует операции с плавающей запятой для процессоров, в которых отсутствует аппаратная поддержка.

Вы также можете взглянуть на libgcc, хотя я считаю, что это GPL, что может или не может быть проблемой для вас.

На самом деле, не просто смотрите на них. Используйте один из них (или другую мягкую библиотеку). Нет необходимости заново изобретать колесо.