Paramètres d'une routine d'intégration IMSL

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

Salut !

Dans la documentation de QAND, regarde la définition exacte des deux paramètres ERRABS et ERRREL. A priori, la valeur de 0.01 pour ERRREL me semble raisonnable. En ce qui concerne ERRABS, ça dépend de tes données. Je te suggère de faire le petit test suivant; tu passes plusieurs fois ton programme avec ERRREL=0.01, la première avec ERRABS=0.1, la deuxième avec 0.2, la troisième avec 0.5, etc. Il ne m'étonnerait pas que ERRABS ne joue aucun rôle a condition qu'il soit assez grand (mais pas trop).

Essaie et donne nous des nouvelles.

Jean-Marc Blanc

Re-bonjour,

Merci de votre réponse.
J'étais justement en train de faire ce genre de tests. Il apparaît que:
- la valeur de ERRABS mène a une divergence des résultats dès qu'elle est en-dessous de 0.1 ;
- la valeur de ERRREL ne joue apparemment aucun rôle: pour ERRABS fixé au-dessus de 0.1, les résultats obtenus sont rigoureusement identiques.

Pour mon programme pas de souci, je sais donc maintenant quel jeu (ERRABS, ERRREL) utiliser - en espérant qu'il fonctionnera avec d'autres cas que celui que je teste - mais j'aimerais comprendre ce résultat plutôt étrange, pour la "culture" !

Salut !

Je n'ai qu'une très antique version de IMSL, où la fonction QAND n'existe pas, mais, pour chaque sous-programme était expliquée très clairement au début du fichier source. C'est là que tu dois regarder.

Jean-Marc Blanc

Je vais y jeter un œil de ce pas !

Merci de votre aide,

julien

Salut,

Il y a un forum dédié à IMSL, maintenu par l'éditeur de cette librairie (Visual Numerics): Forum IMSL.

Tu multiplieras ainsi les chances d'obtenir la solution à ton problème !

SebGR