Esope source files

   1 : C CROSS2    SOURCE    CHAT      05/01/12    22:31:17     5004           
   2 :       SUBROUTINE CROSS2 (V1,V2,V3,IER)                                          
   3 : C=======================================================================        
   4 : C                                                                               
   5 : C     PRODUITS VECTORIELS ET NORMALISATIONS                                     
   6 : C    ENTREES                                                                    
   7 : C       V1(3),V2(3)= VECTEURS                                                   
   8 : C    SORTIES                                                                    
   9 : C       V3(3)= V1 VECTORIEL V2                                                  
  10 : C       IER    = INDICATEUR DE SUCCES 1 SI SUCCES 0 SINON                       
  11 : C=======================================================================        
  12 :       IMPLICIT INTEGER(I-N)
  13 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  14 :       DIMENSION V1(*),V2(*),V3(*)                                               
  15 : C                                                                               
  16 : C                                                                               
  17 :       IER = 1                                                                   
  18 :       V3(1) = V1(2)*V2(3) - V2(2)*V1(3)                                         
  19 :       V3(2) = V1(3)*V2(1) - V2(3)*V1(1)                                         
  20 :       V3(3) = V1(1)*V2(2) - V2(1)*V1(2)                                         
  21 :       RMOD=SQRT(V3(1)*V3(1)+V3(2)*V3(2)+V3(3)*V3(3))                            
  22 : C                                                                               
  23 :       IF (RMOD.NE.0) GO TO 10                                                   
  24 :       IER = 0                                                                   
  25 :       RETURN                                                                    
  26 : C                                                                               
  27 :    10 CONTINUE                                                                  
  28 :       DO 20 I = 1,3                                                             
  29 :    20 V3(I) = V3(I)/RMOD                                                        
  30 :       RETURN                                                                    
  31 :       END                                                                       
  32 :