Esope source files

   1 : C DJ2       SOURCE    BR232186  12/08/01    21:15:24     7456           
   2 : C JACOBA    SOURCE    BR      06/01/10    21:15:29     5284           
   3 :       SUBROUTINE DJ2(A,B)
   4 : C===============================================================
   5 : C Routine permettant de calculer la dérivée du deuxième invariant
   6 : C du déviateur de A 
   7 : C     OBJET                                                                     
   8 : C     -----                                                                     
   9 : C     CALCUL DE B = DJ2(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      = DJ2(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),UNIT(3,3)
  25 : 
  26 :         DO I=1,3
  27 :          DO J=1,3
  28 :             IF (I.EQ.J) THEN
  29 :                UNIT(I,J) = 1.0D0
  30 :             ELSE
  31 :                UNIT(I,J) = 0.0D0
  32 :             ENDIF
  33 :          ENDDO
  34 :       ENDDO
  35 : 
  36 :       TRA_SE = 0.0D0
  37 :       TRA_SE = A(1,1) + A(2,2) + A(3,3)
  38 : 
  39 :       DO I=1,3
  40 :          DO J=1,3
  41 :             TRAV(I,J) = 0.0D0
  42 :             TRAV(I,J) = A(I,J)-(1.0D0/3.0D0)*TRA_SE*UNIT(I,J)
  43 :          ENDDO
  44 :       ENDDO
  45 : 
  46 :       C = 0.0D0
  46 :       CALL J2(A,C)
  47 : 
  48 :       DO I=1,3
  49 :          DO J=1,3
  50 :             B(I,J) = 0.0D0
  51 :             B(I,J) = (3.0D0/2.0D0)*TRAV(I,J)/C
  52 :          ENDDO
  53 :       ENDDO
  54 : 
  55 :       RETURN
  56 :       END
  57 : 
  58 :  
  59 :  
  60 :  
  61 :  
  62 :