Esope source files

   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 :