Esope source files

   1 : C DI2       SOURCE    BR232186  12/08/01    21:15:23     7456           
   2 : C JACOBA    SOURCE    BR      06/01/10    21:15:29     5284           
   3 :       SUBROUTINE DI2(A,B)
   4 : C===============================================================
   5 : C Routine permettant de calculer la dérivée du deuxième invariant
   6 : C de A 
   7 : C     OBJET                                                                     
   8 : C     -----                                                                     
   9 : C     CALCUL DE B = DI2(A)                            
  10 : C                                                                               
  11 : C     ENTREES                                                                   
  12 : C     -------                                                                   
  13 : C     A(3,3) = MATRICE SYMETRIQUE          
  14 : C
  15 : C     SORTIES                                                                   
  16 : C     -------                                                                   
  17 : C     B      = DI2(A) 
  18 : C                                                                               
  19 : C===============================================================                
  20 :       IMPLICIT INTEGER(I-N)
  21 :       IMPLICIT REAL*8 (A-H,O-Z)
  22 :       
  23 :       INTEGER NCAS
  24 :       REAL*8 A(3,3),TRAV(3,3),B(3,3)
  25 : 
  26 :       DO I=1,3
  27 :          DO J=1,3
  28 :             TRAV(I,J) = 0.0D0
  29 :             TRAV(I,J) = A(I,J)
  30 :          ENDDO
  31 :       ENDDO
  32 : 
  33 :       C = 0.0D0
  33 :       CALL I2(A,C)
  34 : 
  35 :       DO I=1,3
  36 :          DO J=1,3
  37 :             B(I,J) = 0.0D0
  38 :             B(I,J) = (3.0D0/2.0D0)*TRAV(I,J)/C
  39 :          ENDDO
  40 :       ENDDO
  41 : 
  42 :       RETURN
  43 :       END
  44 : 
  45 :  
  46 :  
  47 :  
  48 :  
  49 :