Esope source files

   1 : C VMULT4    SOURCE    CHAT      05/01/13    04:08:33     5004           
   2 :       SUBROUTINE VMULT4(A,AA,B)                                                 
   3 : C=======================================================================        
   4 : C     CALCULE LE PRODUIT AA = TRANSPOSE ( B ) * A * B                           
   5 : C   ENTREES                                                                     
   6 : C       A(3,3)    NON SYMETRIQUE                                                
   7 : C       B(3,3)    NON SYMETRIQUE                                                
   8 : C   SORTIES                                                                     
   9 : C       AA(9)     NON SYMETRIQUE  STOCKEE PAR COLONNE                           
  10 : C=======================================================================        
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  13 :       PARAMETER(XZER=0.D0)                                                      
  14 :       DIMENSION A(3,*),AA(*),B(3,*)                                             
  15 :       DIMENSION T(3)                                                            
  16 : C                                                                               
  17 :       IJ=1                                                                      
  18 :       DO 40 J=1,3                                                               
  19 :       DO 20 I1=1,3                                                              
  20 :       C=XZER                                                                    
  21 :       DO 10 J1=1,3                                                              
  22 : 10    C=C+A(I1,J1)*B(J1,J)                                                      
  23 : 20    T(I1)=C                                                                   
  24 :       DO 40  I=1,3                                                              
  25 :       C=XZER                                                                    
  26 :       DO 30  J1=1,3                                                             
  27 : 30    C=C+T(J1)*B(J1,I)                                                         
  28 :       AA(IJ)=C                                                                  
  29 :       IJ=IJ+1                                                                   
  30 :   40  CONTINUE                                                                  
  31 :       RETURN                                                                    
  32 :       END                                                                       
  33 :