Esope source files

   1 : C I2        SOURCE    BR232186  12/08/03    22:02:48     7473           
   2 :       SUBROUTINE I2(A,B)
   3 : C====================================================================== 
   4 : C Routine permettant de calculer le deuxieme invariant 
   5 : C de A 
   6 : C     OBJET                                                                     
   7 : C     -----                                                                     
   8 : C     CALCUL DE B = I2(A)                            
   9 : C                                                                               
  10 : C     ENTREES                                                                   
  11 : C     -------                                                                   
  12 : C     A(3,3) = MATRICE SYMETRIQUE          
  13 : C
  14 : C     SORTIES                                                                   
  15 : C     -------                                                                   
  16 : C     B      = I2(A) 
  17 : C                                                                               
  18 : C===============================================================                
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8 (A-H,O-Z)
  21 :       
  22 :       REAL*8 A(3,3),TRAV(3,3)
  23 : 
  24 :       B = 0.0D0
  25 : 
  26 :       DO I=1,3
  27 :          DO J=1,3
  28 :             TRAV(I,J) = A(I,J)
  29 :          ENDDO
  30 :       ENDDO
  31 : 
  32 :       B = 0.0D0
  33 :       
  34 :       DO I=1,3
  35 :          DO J=1,3
  36 :             B = B + TRAV(I,J)*TRAV(I,J)
  37 :          ENDDO
  38 :       ENDDO
  39 : 
  40 :       B = (3.0D0/2.0D0)*B
  41 :       B = DSQRT(B)
  42 : 
  43 :       RETURN
  44 :       END
  45 : 
  46 :  
  47 :  
  48 :  
  49 :  
  50 :  
  51 :