Télécharger la plaquette Cast3M


fr en

Télécharger alpha.eso

Retour à la liste

Numérotation des lignes : 

alpha

  1. C ALPHA     SOURCE    JC220346  16/11/29    21:15:01     9221           
  2. C                                                                     |
  3.        FUNCTION ALPHA(II,KK,JJ,LL)
  4. C                                                                     |
  5. C      CETTE FONCTION REELLE RENVOIE LA VALEUR DE L'ANGLE ENTRE       |
  6. C      LES 2 FACETTES (KK II JJ) ET (LL II JJ)                        |
  7. C      POUR DES RAISONS DE COUT ON RENVOIE LA COTANGENTE CE QUI       |
  8. C      SUFFIT POUR LES TESTS QUE L'ON EN FERA                         |
  9. C                                                                     |
  10. C---------------------------------------------------------------------|
  11. C
  12.       IMPLICIT INTEGER(I-N)
  13.       IMPLICIT REAL*8(A-H,O-Z)
  14. -INC TDEMAIT
  15.  
  16. -INC PPARAM
  17. -INC CCOPTIO
  18.        DIMENSION V1(3),V2(3),V3(3),VN1(3),VN2(3),VN3(3)
  19. C
  20. C
  21.        DO 10 I=1,3
  22.              V1(I)=XYZ(I,KK)-XYZ(I,II)
  23.              V2(I)=XYZ(I,JJ)-XYZ(I,II)
  24.              V3(I)=XYZ(I,LL)-XYZ(I,II)
  25. 10     CONTINUE
  26. C
  27. C
  28. C      CALCUL DE LA NORME DU VECTEUR V2
  29.        V2NORM=SQRT(V2(1)**2+V2(2)**2+V2(3)**2)
  30.        IF (V2NORM.NE.0.d0) GOTO 15
  31.        IF (IVERB.EQ.1) WRITE (6,1000)
  32. 1000   FORMAT(' UNE NORME EST EGALE A ZERO')
  33.        IF (IVERB.EQ.1) WRITE (6,*) ' II,JJ,KK,LL ',II,JJ,KK,LL
  34.        ALPHA=-1D6
  35.        RETURN
  36. C
  37. C      LE VECTEUR V2 EST NORME
  38. 15     DO 20 I=1,3
  39.              V2(I)=V2(I)/V2NORM
  40. 20     CONTINUE
  41. C
  42. C      VN3=V2^V1
  43.        VN3(1)=V2(2)*V1(3)-V2(3)*V1(2)
  44.        VN3(2)=V2(3)*V1(1)-V2(1)*V1(3)
  45.        VN3(3)=V2(1)*V1(2)-V2(2)*V1(1)
  46. C
  47. C      VN1=V2^V3
  48.        VN1(1)=V2(2)*V3(3)-V2(3)*V3(2)
  49.        VN1(2)=V2(3)*V3(1)-V2(1)*V3(3)
  50.        VN1(3)=V2(1)*V3(2)-V2(2)*V3(1)
  51. C
  52. C      SCALAIRE(VN3,VN1)-----> COS(V1,V3)
  53.        SCAL=VN3(1)*VN1(1)+VN3(2)*VN1(2)+VN3(3)*VN1(3)
  54. C
  55. C      PRODUIT MIXTE (VN3,VN1,V2)---->  SIN(V1,V3)
  56.        PV=     (VN3(2)*VN1(3)-VN3(3)*VN1(2))*V2(1)
  57.      #        +(VN3(3)*VN1(1)-VN3(1)*VN1(3))*V2(2)
  58.      #        +(VN3(1)*VN1(2)-VN3(2)*VN1(1))*V2(3)
  59. C
  60. C
  61.        IF (PV.EQ.0.d0) GOTO 100
  63.        IF(ABS(ALPHA).GT.1d4) GOTO 150
  64.        IF (PV.LT.0.) ALPHA=-2d6+ALPHA
  65.        RETURN
  66. 100    CONTINUE
  67. C  A PRIORI VN2 EST NUL
  68.        ALPHA=-1d6
  69. 150    CONTINUE
  70.        IF (SCAL.GT.0..AND.PV.GT.0.) ALPHA=1E6
  71.        IF (SCAL.GT.0..AND.PV.LE.0.) ALPHA=-3E6
  72.        IF (SCAL.LE.0) ALPHA=-1d6
  73. *      WRITE (6,1111) KK,II,LL
  74. *1111   FORMAT(' ALPHA POINTS ALIGNES ',3I8)
  75. *      WRITE (6,1112) II,KK,JJ,LL
  76. *1112   FORMAT(' ALPHA APPELEE AVEC ',4I6)
  77. *      WRITE (6,1113) (V1(I),I=1,3)
  78. *1113   FORMAT(' V1 ',3G12.5)
  79. *      WRITE (6,1114) (V2(I),I=1,3)
  80. *1114   FORMAT(' V2 ',3G12.5)
  81. *      WRITE (6,1115) (V3(I),I=1,3)
  82. *1115   FORMAT(' V3 ',3G12.5)
  83.        RETURN
  84. C
  85. C      FIN DE LA FONCTION ALPHA
  86.        END
  87.  
  88.  
  89.  
  90.  

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