Télécharger matvec.eso

Retour à la liste

Numérotation des lignes :

matvec
  1. C MATVEC SOURCE CHAT 05/01/13 01:35:46 5004
  2. SUBROUTINE MATVEC(U,V,P,NB)
  3. C=======================================================================
  4. C ROUTINE DE PASSAGE DU VECTEUR DEPLACEMENT GLOBAL AU LOCAL
  5. C
  6. C
  7. C ENTREE
  8. C U(3*NB) = DEPLACEMENTS GLOBAUX AUX 2 NOEUDS
  9. C P(3,3) = MATRICE DE PASSAGE
  10. C NB = NOMBRE DE VECTEURS A 3 LIGNES CONTENUS DANS U ET V
  11. C SORTIE
  12. C V(3*NB) = DEPLACEMENTS LOCAUX AUX 2 NOEUDS
  13. C=======================================================================
  14. IMPLICIT INTEGER(I-N)
  15. IMPLICIT REAL*8(A-H,O-Z)
  16. DIMENSION U(*),V(*),P(3,*)
  17. DATA XZER/0.D0/
  18. C
  19. DO 600 J=1,3
  20. KINC=0
  21. DO 700 K=1,NB
  22. V(J+KINC)=XZER
  23. C
  24. DO 500 I=1,3
  25. V(J+KINC)=V(J+KINC)+P(J,I)*U(I+KINC)
  26. 500 CONTINUE
  27. KINC=KINC+3
  28. 700 CONTINUE
  29. 600 CONTINUE
  30. RETURN
  31. END
  32.  
  33.  

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