Télécharger icgprm.eso

Retour à la liste

Numérotation des lignes :

  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