Télécharger la plaquette Cast3M


fr en

Télécharger trsydr.eso

Retour à la liste

Numérotation des lignes : 

trsydr

  1. C TRSYDR    SOURCE    CHAT      05/01/13    03:53:03     5004
  2.       SUBROUTINE TRSYDR(XE,N, PT1,PT2)
  3. C-----------------------------------------------------------------
  4. C     Calcul de la transfomee d'un ensemble de point par rapport
  5. C     a un point (TRSYPT), une droite (TRSYDR) ou un plan (TRSYPL)
  6. C     de symetrie ou par rapport a une translation (TRTRVE)
  7. C
  8. C     Entree   XE(3,N)     Tableau de point
  9. C
  10. C              PT1(3)      point (DISYPT), point de la droite (DISYDR)
  11. C                          normale NORMEE au plan (DISYPL), vecteur de
  12. C                          la translation
  13. C              PT2(3)      vecteur directeur NORME de la droite (DISYDR)
  14. C              D           distance a l'origine (DISYPL)
  15. C
  16. C     Sortie   XXL2M       Carre de la distance minimale
  17. C
  18. C     Appele par CONNE1 NLOCA1
  19. C
  20. C     PP 23/10/92
  21. C-----------------------------------------------------------------
  22.       IMPLICIT INTEGER(I-N)
  23.       IMPLICIT REAL*8(A-H,O-Z)
  24.       DIMENSION XE(3,N),PT1(3),PT2(3)
  25. C
  26.       DO 3 I=1,N
  27.         A=0.D0
  28.         DO 1 J=1,3
  29.           A=A+PT2(J)*(XE(J,I)-PT1(J))
  30.  1      CONTINUE
  31.         DO 2 J=1,3
  32.           XE(J,I)=XE(J,I)+2*(PT1(J)-XE(J,I)+A*PT2(J))
  33.  2      CONTINUE
  34.  3    CONTINUE
  35.       RETURN
  36. C-----------------------------------------------------------------
  37.       E N T R Y  TRSYPT(XE,N, PT1)
  38. C
  39.       DO 5 I=1,N
  40.         DO 4 J=1,3
  41.           XE(J,I)=XE(J,I)+2*(PT1(J)-XE(J,I))
  42.  4      CONTINUE
  43.  5    CONTINUE
  44.       RETURN
  45. C-----------------------------------------------------------------
  46.       E N T R Y  TRSYPL(XE,N, PT1,D)
  47. C
  48.       DO 8 I=1,N
  49.         B=D
  50.         DO 6 J=1,3
  51.           B=B+XE(J,I)*PT1(J)
  52.  6      CONTINUE
  53.         B=B*2
  54.         DO 7 J=1,3
  55.           XE(J,I)=XE(J,I)-B*PT1(J)
  56.  7      CONTINUE
  57.  8    CONTINUE
  58.       RETURN
  59. C-----------------------------------------------------------------
  60.       E N T R Y  TRTRVE(XE,N, PT1)
  61. C
  62.       DO 10 I=1,N
  63.         DO 9 J=1,3
  64.           XE(J,I)=XE(J,I)+PT1(J)
  65.  9      CONTINUE
  66.  10   CONTINUE
  67.       RETURN
  68. C-----------------------------------------------------------------
  69.       END
  70.  
  71.  

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