Sources Esope | Cast3M

pmvm
C PMVM      SOURCE    CHAT      05/01/13    02:15:29     5004
      SUBROUTINE PMVM(A,IA,JA,X,Y,N,LA)
C
C---  Calcul de y=A*x avec A stock{ en morse. IA ET JA sont les tableaux
C     d'indirection.
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION A(LA),X(1),Y(1)
      DIMENSION IA(N+1),JA(LA)
 
      YT=0.
      DO 10 I=1,N
         YT=0.
         DO 11 J=IA(I),IA(I+1)-1
            YT=YT+A(J)*X(JA(J))
 11      CONTINUE
         Y(I)=YT
 10   CONTINUE
      RETURN
      END