Sources Esope | Cast3M

pmvc
C PMVC      SOURCE    CHAT      05/01/13    02:15:26     5004
      SUBROUTINE PMVC(A,KA,X,Y,N,NNZ)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
C---  Calcul de y=A*x avec A stock{ en mode compress{, KA est le tableau
C     d'indirection. A et KA sont stock{s A(N,NNZ), NNZ nombre maxi de
C     valeurs non nulles par ligne de la matrice.
C
C---  Le code est optimis{ pour Cray (qui ne sait pas vectoriser les
C     boucles internes) et IBM (qui sait lui).
C
      DIMENSION A(N,NNZ),KA(N,NNZ),X(1),Y(1)
 
      YT=0.
      DO 10 I=1,N
         YT=0.D0
         DO 11 J=1,NNZ
            YT=YT+A(I,J)*X(KA(I,J))
 11      CONTINUE
         Y(I)=YT
 10   CONTINUE
      RETURN
      END