Как произвести быстрое вычисление функции exp(x)?

Question

[Pi-314]

Cразу оговорюсь - не программист, но нужда в виде необходимости численного решения системы уравнений мат.физики столкнула с Free Pascal / Lazarus.

Код написал, численная схема работает, но задача ускорения выполнения программы осложнена медленной работой встроенной в FPC функции для вычисления exp(x).

Вопрос гуглил, для Pascal решений не нашел, просьба помочь.

Answer 1

[Pi-314]

Отвечаю сам себе:
Обнаружил статью по теме (www.rsdn.ru/?article/alg/fastpow.xml), изучаю вопрос,скорее всего, нужно расширить предложенное там решение на отрицательные значения аргумента, понятно как это сделать на Pasal, но может кто может доработать сразу ASM ? )) Ау, боги ассемблера!

Из статьи следует, что если Free Pascal в той или иной степени базировался на библиотеках Borland - проблема exp (x) - давняя болезнь Pascal/Delphi...

Answer 2

[Иван Старков]

Бывает в уравнениях мат физики точность некоторых операций на большей части итераций не так уж и важна, так что может поместить в массив предрасчитанные значения e^x на достаточно мелких интервалах, и использовать их в алгоритме до какого то момента пока error_value > N

И только на последней части итераций пользовать "истинное" значение

(я к тому что неточное значение e^x будет как то увеличивать ошибку но это увеличение по сравнению с ошибкой на начальных итерациях вычисления незначительно)

Answer 3

[Иван Старков]

Вобщем это не совсем шутка что в военное время Pi может равняться и трем и даже четырем, для некоторых задач такой точности достаточно