Implementación rápida/aproximación de la función pow() en C/C++

Question

Estoy buscando una implementación más rápida o una buena aproximación de las funciones provistas por cmath.

que han de acelerar las siguientes funciones

1. pow(x,y)
2. exp(z*pow(x,y))

donde z<0. x es de (-1.0,1.0) y y es de (0.0, 5.0)

Answer 1

Éstos son algunos approxmiations:

• Optimized pow Approximation for Java and C/C++. Esta aproximación es muy inexacta, tienes que probar por ti mismo si es lo suficientemente buena.
• Optimized Exponential Functions for Java. ¡Bastante bien! Lo uso para una red neuronal.

Si la aproximación de arriba para pow no es lo suficientemente bueno, todavía puede intentar reemplazarlo con funciones exponenciales, dependiendo de su máquina y compilador de esto podría ser más rápido:

1. x^y = e^(y*ln(x))
2. Y el resultado: e^(z * x^y) = e^(z * e^(y*ln(x)))

Otro truco es cuando algunos parámetros de la fórmula no cambian con frecuencia. Entonces, por ejemplo, xey son en su mayoría constantes, puede precalcular x^y reutilizar esto.

Answer 2

¿Cuáles son los valores posibles de x y y? Si están dentro de límites razonables, la creación de algunas tablas de búsqueda podría ayudar.

Answer 3

Recomiendo las rutinas en el libro "Math Toolkit for Real-Time Programming" de Jack W. Crenshaw.

Es posible que también desee publicar algo de su código para mostrar cómo llama a estas funciones, ya que puede haber algunas otras posibilidades de optimización de nivel superior que no se desprenden de la descripción dada hasta ahora.