Télécharger vpast.eso

Retour à la liste

Numérotation des lignes :

  1. C VPAST SOURCE CHAT 05/01/13 04:09:46 5004
  2. SUBROUTINE VPAST(VCORE,B)
  3. C=======================================================================
  4. C
  5. C CONSTRUCTION DE LA MATRICE DE PASSAGE
  6. C FAISANT PASSER DES COORDONNEES GLOBALES AUX COORDONNEES
  7. C LOCALES POUR UN TRIANGLE OU UN QUADRANGLE( DANS CE CAS ON PREND
  8. C LE PLAN S APPUYANT SUR LES 3 PREMIERS NOEUDS
  9. C ENTREES
  10. C VCORE(3,3)=COORD DU TRIANGLE OU DU QUADRANGLE (1 COL = 1 PT)
  11. C SORTIES
  12. C B(3,3) =MATRICE DE PASSAGE
  13. C=======================================================================
  14. IMPLICIT INTEGER(I-N)
  15. IMPLICIT REAL*8(A-H,O-Z)
  16. DIMENSION VCORE(3,*),B(3,*)
  17. C
  18. C
  19. X21=VCORE(1,2)-VCORE(1,1)
  20. Y21=VCORE(2,2)-VCORE(2,1)
  21. Z21=VCORE(3,2)-VCORE(3,1)
  22. X31=VCORE(1,3)-VCORE(1,1)
  23. Y31=VCORE(2,3)-VCORE(2,1)
  24. Z31=VCORE(3,3)-VCORE(3,1)
  25. VL12=SQRT(X21*X21+Y21*Y21+Z21*Z21)
  26. A1=Y21*Z31-Y31*Z21
  27. A2=X31*Z21-X21*Z31
  28. A3=X21*Y31-X31*Y21
  29. TWOA=SQRT(A1*A1+A2*A2+A3*A3)
  30. C
  31. VL12=1.D0/VL12
  32. TWOA=1.D0/TWOA
  33. B(1,1)=X21*VL12
  34. B(1,2)=Y21*VL12
  35. B(1,3)=Z21*VL12
  36. B(3,1)=A1*TWOA
  37. B(3,2)=A2*TWOA
  38. B(3,3)=A3*TWOA
  39. B(2,1)=B(3,2)*B(1,3)-B(3,3)*B(1,2)
  40. B(2,2)=B(3,3)*B(1,1)-B(1,3)*B(3,1)
  41. B(2,3)=B(3,1)*B(1,2)-B(3,2)*B(1,1)
  42. RETURN
  43. END
  44.  
  45.  

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