Sources Esope | Cast3M

cross2
C CROSS2    SOURCE    CHAT      05/01/12    22:31:17     5004
      SUBROUTINE CROSS2 (V1,V2,V3,IER)
C=======================================================================
C
C     PRODUITS VECTORIELS ET NORMALISATIONS
C    ENTREES
C       V1(3),V2(3)= VECTEURS
C    SORTIES
C       V3(3)= V1 VECTORIEL V2
C       IER    = INDICATEUR DE SUCCES 1 SI SUCCES 0 SINON
C=======================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION V1(*),V2(*),V3(*)
C
C
      IER = 1
      V3(1) = V1(2)*V2(3) - V2(2)*V1(3)
      V3(2) = V1(3)*V2(1) - V2(3)*V1(1)
      V3(3) = V1(1)*V2(2) - V2(1)*V1(2)
      RMOD=SQRT(V3(1)*V3(1)+V3(2)*V3(2)+V3(3)*V3(3))
C
      IF (RMOD.NE.0) GO TO 10
      IER = 0
      RETURN
C
   10 CONTINUE
      DO 20 I = 1,3
   20 V3(I) = V3(I)/RMOD
      RETURN
      END