Télécharger pmvm.eso

Retour à la liste

Numérotation des lignes :

pmvm
  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.  
  24.  

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