Télécharger la plaquette Cast3M


fr en

Télécharger valpr1.eso

Retour à la liste

Numérotation des lignes : 

valpr1

  1. C VALPR1    SOURCE    CHAT      05/01/13    04:00:58     5004
  2.       SUBROUTINE VALPR1 (IPDIAG,IPSOUS,NBITER,PRECI1,PRECI2,IPVALP)
  3. ************************************************************************
  4. *
  5. *                             V A L P R 1
  6. *                             -----------
  7. *
  8. * FONCTION:
  9. * ---------
  10. *
  11. *     VOIR LE SOUS-PROGRAMME "VALPRO".
  12. *
  13. * MODE D'APPEL:
  14. * -------------
  15. *
  16. *     CALL VALPR1 (IPDIAG,IPSOUS,NBITER,PRECI1,PRECI2,IPVALP)
  17. *
  18. * PARAMETRES:   (E)=ENTREE   (S)=SORTIE
  19. * -----------
  20. *
  21. *     IPDIAG  ENTIER    (E)  POINTEUR SUR LE 'LISTREEL' REPRESENTANT LA
  22. *                            LA DIAGONALE DE LA MATRICE TRIDIAGONALE.
  23. *     IPSOUS  ENTIER    (E)  POINTEUR SUR LE 'LISTREEL' REPRESENTANT LA
  24. *                            SOUS-DIAGONALE DE LA MATRICE TRIDIAGONALE.
  25. *     NBITER  ENTIER    (E)  NOMBRE MAXIMAL D'ITERATIONS PERMISES.
  26. *                       (S)  VALEUR DE REMPLACEMENT FOURNIE SI "NBITER"
  27. *                            NON DONNE EN ENTREE.
  28. *     PRECI1  REEL DP   (E)  PRECISION ABSOLUE DE CONVERGENCE.
  29. *                       (S)  VALEUR DE REMPLACEMENT FOURNIE SI NI
  30. *                            "PRECI1", NI "PRECI2", N'ONT ETE DONNES EN
  31. *                            ENTREE.
  32. *     PRECI2  REEL DP   (E)  PRECISION RELATIVE DE CONVERGENCE.
  33. *                       (S)  VALEUR DE REMPLACEMENT FOURNIE SI NI
  34. *                            "PRECI1", NI "PRECI2", N'ONT ETE DONNES EN
  35. *                            ENTREE.
  36. *     IPVALP  ENTIER    (S)  POINTEUR SUR LE 'LISTREEL' CONTENANT LES
  37. *                            VALEURS PROPRES, DANS L'ORDRE DE CALCUL.
  38. *
  39. * SOUS-PROGRAMMES APPELES:
  40. * ------------------------
  41. *
  42. *     VALPR2
  43. *
  44. * AUTEUR, DATE DE CREATION:
  45. * -------------------------
  46. *
  47. *     PASCAL MANIGOT     27 FEVRIER 1985
  48. *
  49. * LANGAGE:
  50. * --------
  51. *
  52. *     ESOPE + FORTRAN77
  53. *
  54. ************************************************************************
  55. *
  56.       IMPLICIT INTEGER(I-N)
  57.  
  58. -INC PPARAM
  59. -INC CCOPTIO
  60. -INC SMLREEL
  61. *
  62.       REAL*8 PRECI1,PRECI2
  63. *
  64.       PARAMETER (NBRITE = 50)
  65.       REAL*8 PREC1,PREC2
  66.       PARAMETER (PREC1 = 1.D-30 , PREC2 = 0.D0)
  67. *
  68.       IF (NBITER .LE. 0) THEN
  69.          NBITER = NBRITE
  70.       END IF
  71.       IF (PRECI1 .LE. 0.D0 .AND. PRECI2 .LE. 0.D0) THEN
  72.          PRECI1 = PREC1
  73.          PRECI2 = PREC2
  74.       END IF
  75. *
  76.       IF (IIMPI .EQ. 732) THEN
  77.          WRITE (IOIMP,2000) NBITER,PRECI1,PRECI2
  78.  2000    FORMAT (//' LE CALCUL DE VALEURS PROPRES VA ETRE FAIT AVEC ',
  79.      &         ' UN NOMBRE MAXI D''ITERATIONS DE ',I3,', '/
  80.      &         ' UNE PRECISION ABSOLUE DE ',1PE10.3,' ET UNE ',
  81.      &         'PRECISION RELATIVE DE ',1PE10.3,' .'//)
  82.       END IF
  83. *
  84.       MLREEL = IPDIAG
  85.       MLREE1 = IPSOUS
  86.       SEGINI,MLREE2=MLREEL
  87.       IPVALP = MLREE2
  88.       SEGINI,MLREE3=MLREE1
  89. *
  90.       NTERME = MLREE2.PROG(/1)
  91. *
  92.       CALL VALPR2 (MLREE2.PROG,NTERME,MLREE3.PROG,NBITER,PRECI1,PRECI2)
  93.       IF (IERR .NE. 0) RETURN
  94. *
  95.       IF (IIMPI .EQ. 732) THEN
  96.          WRITE (IOIMP,2010) (MLREE3.PROG(I),I=1,(NTERME-1))
  97.  2010    FORMAT (//' TERMES DE LA SOUS-DIAGONALE APRES CONVERGENCE:'//
  98.      &          10(2X,1PE10.3))
  99.       END IF
  100. *
  101.       SEGSUP,MLREE3
  102.       SEGDES,MLREE2
  103. *
  104.       END
  105.  
  106.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales