Esope source files

   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 :