Bonjour !
Dans mon programme FORTRAN, j'utilise la librairie IMSL pour effectuer un calcul de triple intégrale (subroutine QAND). Je suis en phase de test donc les résultats obtenus sont comparés avec des résultats analytiques calculés à la main.
La routine QAND possède deux paramètres d'entrée ERRABS et ERRREL qui décrivent les erreurs absolue et relative désirées par l'utilisateur. Lorsque ERRABS=ERRREL=0.1 aucun problème, les résultats concordent avec la théorie. Par contre pour ERRABS = ERRREL = 0.01 (précision attendue plus grande), les résultats divergent. Je travaille en réels double précision.
Impossible de comprendre d'où vient ce comportement. Quelqu'un saurait-il m'aider ?
Merci d'avance !
Julien
Partager