log1p, log1pf, log1pl

Параметры

Возвращаемое значение

Если ошибок не возникает, ln(1 + arg) возвращается.

Если происходит ошибка области определения, возвращается значение, определяемое реализацией (NaN, если поддерживается).

Если возникает ошибка полюса, - HUGE_VAL , -HUGE_VALF , или -HUGE_VALL возвращается.

Если происходит ошибка диапазона из-за потери значимости (underflow), возвращается корректный результат (после округления).

Обработка ошибок

Ошибки сообщаются, как указано в math_errhandling .

Ошибка области определения возникает, если arg меньше -1 .

Ошибка полюса может возникнуть, если arg равен -1 .

Если реализация поддерживает арифметику с плавающей запятой IEEE (IEC 60559),

• Если аргумент равен ±0, он возвращается без изменений.
• Если аргумент равен -1, возвращается -∞ и FE_DIVBYZERO устанавливается.
• Если аргумент меньше -1, возвращается NaN и FE_INVALID устанавливается.
• Если аргумент равен +∞, возвращается +∞.
• Если аргумент равен NaN, возвращается NaN.

Примечания

Функции expm1 и log1p полезны для финансовых расчетов, например, при вычислении небольших ежедневных процентных ставок: (1+x) n
-1 может быть выражено как expm1 ( n * log1p ( x ) ) . Эти функции также упрощают написание точных обратных гиперболических функций.

Пример

#include <errno.h>
#include <fenv.h>
#include <float.h>
#include <math.h>
#include <stdio.h>
// #pragma STDC FENV_ACCESS ON
int main(void)
{
    printf("log1p(0) = %f\n", log1p(0));
    printf("Процентный доход за 2 дня на $100, начисляемый ежедневно по 1%%\n"
           " по календарю 30/360 = %f\n",
           100*expm1(2*log1p(0.01/360)));
    printf("log(1+1e-16) = %g, но log1p(1e-16) = %g\n",
           log(1+1e-16), log1p(1e-16));
    // специальные значения
    printf("log1p(-0) = %f\n", log1p(-0.0));
    printf("log1p(+Inf) = %f\n", log1p(INFINITY));
    // обработка ошибок
    errno = 0; feclearexcept(FE_ALL_EXCEPT);
    printf("log1p(-1) = %f\n", log1p(-1));
    if (errno == ERANGE)
        perror("    errno == ERANGE");
    if (fetestexcept(FE_DIVBYZERO))
        puts("    FE_DIVBYZERO raised");
}

log1p(0) = 0.000000
Процентный доход за 2 дня на $100, начисляемый ежедневно по 1%
 по календарю 30/360 = 0.005556
log(1+1e-16) = 0, но log1p(1e-16) = 1e-16
log1p(-0) = -0.000000
log1p(+Inf) = Inf
log1p(-1) = -Inf
    errno == ERANGE: Result too large
    FE_DIVBYZERO raised

Ссылки

Смотрите также