log1p, log1pf, log1pl

Paramètres

Valeur de retour

Si aucune erreur ne se produit ln(1 + arg) est retourné.

Si une erreur de domaine se produit, une valeur définie par l'implémentation est retournée (NaN là où supporté).

Si une erreur de pôle se produit, - HUGE_VAL , -HUGE_VALF , ou -HUGE_VALL est retourné.

Si une erreur de plage se produit en raison d'un dépassement inférieur, le résultat correct (après arrondi) est retourné.

Gestion des erreurs

Les erreurs sont signalées comme spécifié dans math_errhandling .

Une erreur de domaine se produit si arg est inférieur à -1 .

Une erreur de pôle peut se produire si arg est -1 .

Si l'implémentation prend en charge l'arithmétique à virgule flottante IEEE (IEC 60559),

• Si l'argument est ±0, il est retourné inchangé.
• Si l'argument est -1, -∞ est retourné et FE_DIVBYZERO est déclenché.
• Si l'argument est inférieur à -1, NaN est retourné et FE_INVALID est déclenché.
• Si l'argument est +∞, +∞ est retourné.
• Si l'argument est NaN, NaN est retourné.

Notes

Les fonctions expm1 et log1p sont utiles pour les calculs financiers, par exemple, lors du calcul de petits taux d'intérêt quotidiens : (1+x) n
-1 peut être exprimé comme expm1 ( n * log1p ( x ) ) . Ces fonctions simplifient également l'écriture de fonctions hyperboliques inverses précises.

Exemple

#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("Intérêts gagnés en 2 jours sur 100$, composés quotidiennement à 1%%\n"
           " sur un calendrier 30/360 = %f\n",
           100*expm1(2*log1p(0.01/360)));
    printf("log(1+1e-16) = %g, mais log1p(1e-16) = %g\n",
           log(1+1e-16), log1p(1e-16));
    // valeurs spéciales
    printf("log1p(-0) = %f\n", log1p(-0.0));
    printf("log1p(+Inf) = %f\n", log1p(INFINITY));
    // gestion d'erreur
    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
Intérêts gagnés en 2 jours sur 100$, composés quotidiennement à 1%
 sur un calendrier 30/360 = 0.005556
log(1+1e-16) = 0, mais log1p(1e-16) = 1e-16
log1p(-0) = -0.000000
log1p(+Inf) = Inf
log1p(-1) = -Inf
    errno == ERANGE: Résultat trop grand
    FE_DIVBYZERO raised

Références

Voir aussi