Télécharger pmvc.eso

Retour à la liste

Numérotation des lignes :

  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.  
  27.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales