Télécharger la plaquette Cast3M


fr en

Télécharger icgprm.eso

Retour à la liste

Numérotation des lignes : 

icgprm

  1. C ICGPRM    SOURCE    CHAT      05/01/13    00:28:12     5004
  2.       SUBROUTINE ICGPRM(X,Y,P,IA,JA,NEL)
  3. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  4. C                                                                   C
  5. C  Effectue une r‚solution par descente remont‚e de la matrice      C
  6. C  de pr‚conditionnement.                                           C
  7. C  La matrice est stock‚e en morse dans P.                          C
  8. C                                                                   C
  9. C +---------------------------------------------------------------+ C
  10. C ! DESCRIPTION DES ARGUMENTS                                     ! C
  11. C !---------------------------------------------------------------! C
  12. C !                                                               ! C
  13. C +----------+------+---------------------------------------------+ C
  14. C ! NOM      ! MODE ! DESCRIPTION                                 ! C
  15. C +----------+------+---------------------------------------------+ C
  16. C ! X        ! <--- ! Second membre                               ! C
  17. C +----------+------+---------------------------------------------+ C
  18. C ! Y        ! ---> ! Solution                                    ! C
  19. C +----------+------+---------------------------------------------+ C
  20. C ! P        ! <--- ! La matrice de pr‚conditionnement            ! C
  21. C +----------+------+---------------------------------------------+ C
  22. C ! IA       ! <--- ! Le tableau de pointeurs                     ! C
  23. C +----------+------+---------------------------------------------+ C
  24. C ! JA       ! <--- ! Le tableau d'indexage correspondant         ! C
  25. C +----------+------+---------------------------------------------+ C
  26. C ! NEL      ! <--- ! Nombre de lignes de la matrice              ! C
  27. C +----------+------+---------------------------------------------+ C
  28. C !                                                               ! C
  29. C ! MODE = ---> VARIABLE ECRITE (SORTIE)                          ! C
  30. C ! MODE = <--- VARIABLE LUE    (ENTREE)                          ! C
  31. C ! MODE = <--> VARIABLE LUE ET ECRITE                            ! C
  32. C !                                                               ! C
  33. C +---------------------------------------------------------------+ C
  34. C !                                                               ! C
  35. C ! CONTACT : Fran‡ois Thomas  69 08 24 64 ou 76 96               ! C
  36. C !                                                               ! C
  37. C +---------------------------------------------------------------+ C
  38. C !                                                               ! C
  39. C ! RELATIONS AVEC D'AUTRES MODULES                               ! C
  40. C !                                                               ! C
  41. C +---------------------------------------------------------------+ C
  42. C !                                                               ! C
  43. C ! MODIFICATIONS                                                 ! C
  44. C !                                                               ! C
  45. C ! 11/03/92 : Cr‚ation                                           ! C
  46. C !                                                               ! C
  47. C +---------------------------------------------------------------+ C
  48. C                                                                   C
  49. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  50.       IMPLICIT INTEGER(I-N)
  51.       IMPLICIT REAL*8 (A-H,O-Z)
  52.       DIMENSION X(NEL)
  53.       DIMENSION Y(NEL)
  54.       DIMENSION P(*)
  55.       DIMENSION IA(NEL+1)
  56.       DIMENSION JA(*)
  57. C
  58.       DO 8 I=1,NEL
  59.          X(I)=Y(I)
  60.  8    CONTINUE
  61.  
  62.       DO 10 I=1,NEL
  63.          S=0.D0
  64.          DO 11 IKA=IA(I)+1,IA(I+1)-1
  65.             K=JA(IKA)
  66.             IF(K.LT.I) S=S+P(IKA)*X(K)
  67.  11      CONTINUE
  68.          X(I)=X(I)-S
  69.  10   CONTINUE
  70.  
  71.       DO 20 I=1,NEL
  72.          X(I)=X(I)*P(IA(I))
  73.  20   CONTINUE
  74.  
  75.       DO 30 II=0,NEL-1
  76.          I=NEL-II
  77.          S=0.D0
  78.          IF(I.GT.0) THEN
  79.          DO 31 IKA=IA(I)+1,IA(I+1)-1
  80.             K=JA(IKA)
  81.             IF(K.GT.I) S=S+P(IKA)*X(K)
  82.  31      CONTINUE
  83.          ENDIF
  84.          X(I)=X(I)-S
  85.  30   CONTINUE
  86.  
  87.       RETURN
  88.       END
  89.  
  90.  
  91.  

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