1 : C PMVC SOURCE CHAT 05/01/13 02:15:26 5004
2 : SUBROUTINE PMVC(A,KA,X,Y,N,NNZ)
3 : IMPLICIT INTEGER(I-N)
4 : IMPLICIT REAL*8 (A-H,O-Z)
5 : C
6 : C--- Calcul de y=A*x avec A stock{ en mode compress{, KA est le tableau
7 : C d'indirection. A et KA sont stock{s A(N,NNZ), NNZ nombre maxi de
8 : C valeurs non nulles par ligne de la matrice.
9 : C
10 : C--- Le code est optimis{ pour Cray (qui ne sait pas vectoriser les
11 : C boucles internes) et IBM (qui sait lui).
12 : C
13 : DIMENSION A(N,NNZ),KA(N,NNZ),X(1),Y(1)
14 :
15 : YT=0.
16 : DO 10 I=1,N
17 : YT=0.D0
18 : DO 11 J=1,NNZ
19 : YT=YT+A(I,J)*X(KA(I,J))
20 : 11 CONTINUE
21 : Y(I)=YT
22 : 10 CONTINUE
23 : RETURN
24 : END
25 :
26 :
© Cast3M 2003 - All rights reserved.
Disclaimer