std::expm1, std::expm1f, std::expm1l - cppreference.com

De cppreference.com

Definido en el archivo de encabezado `<cmath>`
`float expm1 ( float arg ); float expm1f( float arg );`	(1)	(desde C++11)
`double expm1 ( double arg );`	(2)	(desde C++11)
`long double expm1 ( long double arg ); long double expm1l( long double arg );`	(3)	(desde C++11)
`double expm1 ( TipoEntero arg );`	(4)	(desde C++11)

1-3) Calcula e (el número de Euler, 2.7182818) elevado a la potencia dada arg, menos 1.0. Esta función es más precisa que la expresión std::exp(arg)-1.0 si arg está cerca de cero.

4) Un conjunto de sobrecargas o una plantilla de función que acepta un argumento de cualquier tipo entero. Equivalente a (2) (el argumento se convierte a double).

Parámetros

arg

Valor de un tipo de punto flotante o un tipo entero.

Valor de retorno

earg
-1 si no se producen errores.

Si se produce un error de rango debido a desbordamiento, se devuelve HUGE_VAL, +HUGE_VALF, o +HUGE_VALL.

Si se produce un error debido a subdesbordamiento, se devuelve el resultado correcto (después del redondeo).

Manejo de errores

Los errores se informan como se especifica en math_errhandling.

Si la implementación admite la aritmética de punto flotante IEEE (IEC 60559):

Si el argumento es +0, se devuelve +0.
Si el argumento es -0, se devuelve -0.
Si el argumento es -∞, se devuelve -1.
Si el argumento es +∞, se devuelve +∞.
Si el argumento es NaN, se devuelve NaN.

Notas

Las funciones std::expm1 y std::log1p son útiles para cálculos financieros, por ejemplo, al calcular pequeñas tasas de interés diarias: (1+x)n
-1 se puede expresar como std::expm1(n * std::log1p(x)). Estas funciones también simplifican la escritura de funciones hiperbólicas inversas precisas.

Para el tipo double compatible con IEEE, se garantiza el desbordamiento si 709.8 < arg.

Ejemplo

#include <iostream>
#include <cmath>
#include <cerrno>
#include <cstring>
#include <cfenv>
// #pragma STDC FENV_ACCESS ON
int main()
{
    std::cout << "expm1(1) = " << std::expm1(1) << '\n'
              << "Interés devengado en 2 días sobre $100, compuesto diariamente al 1%\n"
              << " en un calendario 30/360 = "
              << 100*std::expm1(2*std::log1p(0.01/360)) << '\n'
              << "exp(1e-16)-1 = " << std::exp(1e-16)-1
              << ", but expm1(1e-16) = " << std::expm1(1e-16) << '\n';
    // valores especiales
    std::cout << "expm1(-0) = " << std::expm1(-0.0) << '\n'
              << "expm1(-Inf) = " << std::expm1(-INFINITY) << '\n';
    // manejo de errores
    errno = 0;
    std::feclearexcept(FE_ALL_EXCEPT);
    std::cout << "expm1(710) = " << std::expm1(710) << '\n';
    if (errno == ERANGE)
        std::cout << "    errno == ERANGE: " << std::strerror(errno) << '\n';
    if (std::fetestexcept(FE_OVERFLOW))
        std::cout << "    Se generó FE_OVERFLOW\n";
}

Posible salida: