Télécharger pourot.eso

Retour à la liste

Numérotation des lignes :

  1. C POUROT SOURCE CHAT 05/01/13 02:20:19 5004
  2. SUBROUTINE POUROT(REL,LRE,P,RP)
  3. C-----------------------------------------------------------------------
  4. C ROUTINE DE PASSAGE D UNE MATRICE ELEMENTAIRE DES AXES LOCAUX AUX AXES
  5. C GLOBAUX
  6. C D APRES TEDEL
  7. C
  8. C ENTREE
  9. C LRE NOMBRE DE LIGNES DU TABLEAU REL
  10. C P MATRICE DE PASSAGE
  11. C RP MATRICE ELEMENTAIRE 12X12 EN AXES LOCAUX
  12. C
  13. C SORTIE
  14. C REL MATRICE ELEMENTAIRE EN AXES GLOBAUX
  15. C
  16. C VIENT DE TEDEL AOUT 84
  17. C-----------------------------------------------------------------------
  18. IMPLICIT INTEGER(I-N)
  19. IMPLICIT REAL*8(A-H,O-Z)
  20. DIMENSION REL(LRE,1),P(3,3),RP(12,12)
  21. DIMENSION NI(12),NJ(12)
  22. C
  23. DATA NI/1,1,1,4,4,4,7,7,7,10,10,10/
  24. DATA NJ/1,2,3,1,2,3,1,2,3,1,2,3/
  25. DO 10 K=1,12
  26. I3=NI(K)
  27. I4=I3+2
  28. LK=NJ(K)
  29. DO 11 L=K,12
  30. REL(K,L)= 0.
  31. I1=NI(L)
  32. I2=I1+2
  33. LL=NJ(L)
  34. LM=0
  35. DO 12 M=I3,I4
  36. LM=LM+1
  37. S=0.
  38. II=0
  39. DO 13 N=I1,I2
  40. II=II+1
  41. 13 S=S+ RP(M,N)*P(II,LL)
  42. 12 REL(K,L)= REL(K,L)+ S*P( LM,LK)
  43. 11 REL(L,K)= REL(K,L)
  44. 10 CONTINUE
  45. RETURN
  46. END
  47.  
  48.  

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