Esope source files

   1 : C DBLECONT  SOURCE    CHAT      12/08/06    21:15:20     7479           
   2 :                  SUBROUTINE DBLECONT(A,B,C)
   3 : C====================================================================== 
   4 : C Routine permettant de calculer le produit doublement contracté de
   5 : C deux tenseurs C = A : B
   6 : C     OBJET                                                                     
   7 : C     -----                                                                     
   8 : C     CALCUL DE C = A : B                              
   9 : C                                                                               
  10 : C     ENTREES                                                                   
  11 : C     -------                                                                   
  12 : C     A(3,3) = MATRICE SYMETRIQUE          
  13 : C     B(3,3) = MATRICE SYMETRIQUE                                                                                    
  14 : C
  15 : C     SORTIES                                                                   
  16 : C     -------                                                                   
  17 : C     C      = PRODUIT DOUBLEMENT CONTRACTE
  18 : C                                                                               
  19 : C===============================================================                
  20 :       IMPLICIT INTEGER(I-N)
  21 :       IMPLICIT REAL*8 (A-H,O-Z)
  22 : 
  23 :       REAL*8 A(3,3),B(3,3)
  24 : 
  25 :       C = 0.0D0
  26 : 
  27 :       DO I=1,3
  28 :          DO J=1,3
  29 :             C = C + A(I,J)*B(I,J)
  30 :          ENDDO
  31 :       ENDDO
  32 : 
  33 :       END
  34 : 
  35 :  
  36 :  
  37 :  
  38 :  
  39 :