Esope source files

   1 : C ELFREE    SOURCE    CHAT      05/01/12    23:32:33     5004           
   2 :       SUBROUTINE ELFREE(V,A,D)                                                  
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 : C                                                                               
   6 : C  ====================================================================         
   7 : C    SOUS-PROGRAMME FORTRAN APPELE PAR ELFRES POUR LE CALCUL DU PRODUIT         
   8 : C    DE CONVOLUTION PARTIEL SUR 1 ELEMENT ET 1 PAS DE TEMPS :                   
   9 : C                                                                               
  10 : C          VN=AN*DN+BN*CN                                                       
  11 : C                                                                               
  12 : C          V  : VECTEUR SECOND MEMBRE VN                                        
  13 : C          A  : TERMES DES MATRICES AN ET BN                                    
  14 : C          D  : VECTEUR DN ET CN                                                
  15 : C                                                                               
  16 : C    CREATION    :  3/11/86                                                     
  17 : C    PROGRAMMEUR : GUILBAUD                                                     
  18 : C  ====================================================================         
  19 : C                                                                               
  20 :       DIMENSION V(*),A(*),D(*)                                                  
  21 : C                                                                               
  22 :       V( 1)=V( 1)+A( 2)*D( 1)            +A(16)*D(13)                           
  23 :       V( 2)=V( 2)+A( 8)*D( 2)+A( 7)*D( 6)+A(22)*D(14)-A(21)*D(18)               
  24 :       V( 3)=V( 3)+A(13)*D( 3)-A(12)*D( 5)+A(27)*D(15)+A(26)*D(17)               
  25 :       V( 4)=V( 4)            +A( 4)*D( 4)            +A(18)*D(16)               
  26 :       V( 5)=V( 5)-A(14)*D( 3)+A(13)*D( 5)+A(28)*D(15)+A(27)*D(17)               
  27 :       V( 6)=V( 6)+A( 9)*D( 2)+A( 8)*D( 6)-A(23)*D(14)+A(22)*D(18)               
  28 :       V( 7)=V( 7)+A(16)*D( 1)            +A( 2)*D(13)                           
  29 :       V( 8)=V( 8)+A(22)*D( 2)+A(21)*D( 6)+A( 8)*D(14)-A( 7)*D(18)               
  30 :       V( 9)=V( 9)+A(27)*D( 3)-A(26)*D( 5)+A(13)*D(15)+A(12)*D(17)               
  31 :       V(10)=V(10)            +A(18)*D( 4)            +A( 4)*D(16)               
  32 :       V(11)=V(11)-A(28)*D( 3)+A(27)*D( 5)+A(14)*D(15)+A(13)*D(17)               
  33 :       V(12)=V(12)+A(23)*D( 2)+A(22)*D( 6)-A( 9)*D(14)+A( 8)*D(18)               
  34 : C                                                                               
  35 :       V( 1)=V( 1)+A( 1)*D( 7)            +A(15)*D(19)                           
  36 :       V( 2)=V( 2)+A( 5)*D( 8)            +A(19)*D(20)+A(20)*D(24)               
  37 :       V( 3)=V( 3)+A(10)*D( 9)            +A(24)*D(21)-A(25)*D(23)               
  38 :       V( 4)=V( 4)            +A( 3)*D(10)            +A(17)*D(22)               
  39 :       V( 5)=V( 5)            -A(12)*D(11)+A(25)*D(21)-A(26)*D(23)               
  40 :       V( 6)=V( 6)            -A( 7)*D(12)-A(20)*D(20)-A(21)*D(24)               
  41 :       V( 7)=V( 7)+A(15)*D( 7)            +A( 1)*D(19)                           
  42 :       V( 8)=V( 8)+A(19)*D( 8)-A(20)*D(12)+A( 5)*D(20)                           
  43 :       V( 9)=V( 9)+A(24)*D( 9)+A(25)*D(11)+A(10)*D(21)                           
  44 :       V(10)=V(10)            +A(17)*D(10)            +A( 3)*D(22)               
  45 :       V(11)=V(11)-A(25)*D( 9)-A(26)*D(11)            -A(12)*D(23)               
  46 :       V(12)=V(12)+A(20)*D( 8)-A(21)*D(12)            -A( 7)*D(24)               
  47 : C                                                                               
  48 : C     WRITE(6,*) ' AN '                                                         
  49 : C     WRITE(6,1001)(A(J),J=1,28)                                                
  50 : C     WRITE(6,*) ' DN '                                                         
  51 : C     WRITE(6,1001)(D(J),J=1,24)                                                
  52 : C     WRITE(6,*) ' VN '                                                         
  53 : C     WRITE(6,1001)(V(J),J=1,12)                                                
  54 : C1001 FORMAT(1X,12(1PE9.2,1X))                                                  
  55 :       RETURN                                                                    
  56 :       END                                                                       
  57 :