std::log1p, std::log1pf, std::log1pl

[editar] Parámetros

[editar] Valor de retorno

ln(1+arg) si no se producen errores.

Si se produce un error de dominio, se devuelve un valor definido por la implementación (NaN donde se dé apoyo).

Si se produce un error de polo, 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).

[editar] Manejo de errores

Los errores se informan como se especifica en math_errhandling.

Se produce un error de dominio si arg es menor que -1.

Se puede producir un error de polo si arg es -1.

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 -1, se devuelve -∞ y se genera FE_DIVBYZERO.
• Si el argumento es menor que -1, se devuelve NaN y se genera FE_INVALID.
• Si el argumento es +∞, se devuelve +∞.
• Si el argumento es NaN, se devuelve NaN.

[editar] 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.

[editar] Ejemplo

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

log1p(0) = 0
Interés devengado en 2 días sobre $100, compuesto diariamente al 1%
 en un calendario 30/360 = 0.00555563
log(1+1e-16) = 0 log1p(1e-16) = 1e-16
log1p(-0) = -0
log1p(+Inf) = inf
log1p(-1) = -inf
    errno == ERANGE: Resultado demasiado grande
    Se generó FE_DIVBYZERO

[editar] Véase también