[an error occurred while processing this directive]
|
быстрое приближенное вычисление 1/x (с небольшой точностью) с дальнейшем уточнением результата методом Ньютона-Рапсона: y[n+1] = y[n]*(2-x*y[n]). Каждая такая итерация увеличивает количество точных разрядов вдвое. То есть, приблизительно посчитав 1/x с точностью 8 разрядов (например по таблице), за 2 таких итерации можно получить точность 32 бита - что, как правило, достаточно. А время, занимаемое одной такой итерацией, на большинстве ДСП обычно 2-3 такта процессора.
E-mail: info@telesys.ru