1 : C PMVM SOURCE CHAT 05/01/13 02:15:29 5004
2 : SUBROUTINE PMVM(A,IA,JA,X,Y,N,LA)
3 : C
4 : C--- Calcul de y=A*x avec A stock{ en morse. IA ET JA sont les tableaux
5 : C d'indirection.
6 : C
7 : IMPLICIT INTEGER(I-N)
8 : IMPLICIT REAL*8 (A-H,O-Z)
9 : DIMENSION A(LA),X(1),Y(1)
10 : DIMENSION IA(N+1),JA(LA)
11 :
12 : YT=0.
13 : DO 10 I=1,N
14 : YT=0.
15 : DO 11 J=IA(I),IA(I+1)-1
16 : YT=YT+A(J)*X(JA(J))
17 : 11 CONTINUE
18 : Y(I)=YT
19 : 10 CONTINUE
20 : RETURN
21 : END
22 :
23 :
© Cast3M 2003 - All rights reserved.
Disclaimer