SUBROUTINE FKMUL
*-- Author : S.Burke / J.V. Morris
      SUBROUTINE FKMUL(Q1,D,Q)
**********************************************************************                                        
*                                                                    *                                        
* Project error matrix                                               *                                        
*                                                                    *                                        
**********************************************************************                                        
                                                                        
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)                               
                                                                        
*KEEP,FKNPL.                                                                                                  
      CHARACTER*5 CKDBG                                                 
      PARAMETER (CKDBG='FKDBG')                                         
      PARAMETER (NPL=72)                                                
      LOGICAL LTRUE,LFIRST,LTRPL,LTRPLD                                 
      DOUBLE PRECISION TRUE,RTRUE,CHITRU,SPRO,CPRO,SFIL,CFIL            
     &,                SSMT,CSMT,SSMTR,CSMTR,DPRO,CBPRO,QPRO,QGAIN      
     &,                RPRO,CRPRO,RFIL,CRFIL,RSMT,CRSMT,CHIFIL,CHISMT   
*                                                                                                             
* Per-track values can go in H1WORK; note that LTRUE and LFIRST must                                          
* be set at least per event.                                                                                  
*                                                                                                             
* This is about 36k words long; the remaining common blocks are                                               
* about 3.6k in total. Some of this could be in /H1WORK/, but the                                             
* blocks would have to be reorganised.                                                                        
*                                                                                                             
      COMMON /H1WORK/                                                   
* /FKPROJ/                                                                                                    
     &                SPRO(5,NPL),CPRO(5,5,NPL)                         
* /FKFILT/                                                                                                    
     &,               SFIL(5,NPL),CFIL(5,5,NPL)                         
* /FKSMTH/                                                                                                    
     &,               SSMT(5,NPL),CSMT(5,5,NPL)                         
     &,               SSMTR(5,NPL),CSMTR(5,5,NPL)                       
* /FKINT/                                                                                                     
     &,               DPRO(5,5,NPL),CBPRO(5,5,NPL),QPRO(5,5,NPL)        
     &,               QGAIN(5,5,NPL),IAPROX,LFIRST                      
* /FKRSID/                                                                                                    
     &,               RPRO(2,NPL),CRPRO(2,2,NPL),RFIL(2,NPL)            
     &,               CRFIL(2,2,NPL),RSMT(2,NPL),CRSMT(2,2,NPL)         
     &,               CHIFIL(NPL),CHISMT(NPL)                           
* /FKTRUE/                                                                                                    
     &,               TRUE(5,NPL),RTRUE(5,NPL),CHITRU(NPL),LTRUE        
* /FKDBG/                                                                                                     
     &,               LTRPL(NPL),LTRPLD(NPL)                            
*KEEP,FKINT.                                                                                                  
*KEND.                                                                                                        
                                                                        
      DIMENSION Q1(5,5),D(5,5),Q(5,5)                                   
                                                                        
**********************************************************************                                        
                                                                        
      DQ14 = Q1(4,1) + D(1,4)*Q1(4,4) + D(1,5)*Q1(5,4)                  
      DQ15 = Q1(5,1) + D(1,4)*Q1(5,4) + D(1,5)*Q1(5,5)                  
                                                                        
      DQ24 = Q1(4,2) + D(2,4)*Q1(4,4) + D(2,5)*Q1(5,4)                  
      DQ25 = Q1(5,2) + D(2,4)*Q1(5,4) + D(2,5)*Q1(5,5)                  
                                                                        
      IF (IAPROX.LE.2) THEN                                             
         Q(1,1) = Q1(1,1) + D(1,4)*(Q1(4,1)+DQ14)                       
     &                    + D(1,5)*(Q1(5,1)+DQ15)                       
         Q(2,1) = Q1(2,1) + D(1,4)*Q1(4,2) + D(1,5)*Q1(5,2)             
     &                    + D(2,4)*DQ14    + D(2,5)*DQ15                
         Q(2,2) = Q1(2,2) + D(2,4)*(Q1(4,2)+DQ24)                       
     &                    + D(2,5)*(Q1(5,2)+DQ25)                       
         Q(4,1) = DQ14                                                  
         Q(4,2) = DQ24                                                  
      ENDIF                                                             
                                                                        
      Q(3,3) = Q1(3,3)                                                  
                                                                        
      Q(4,3) = Q1(4,3)                                                  
      Q(4,4) = Q1(4,4)                                                  
                                                                        
      IF (IAPROX.EQ.1) THEN                                             
         Q(3,1) = Q1(3,1) + D(1,4)*Q1(4,3) + D(1,5)*Q1(5,3)             
         Q(3,2) = Q1(3,2) + D(2,4)*Q1(4,3) + D(2,5)*Q1(5,3)             
         Q(4,1) = DQ14                                                  
         Q(4,2) = DQ24                                                  
         Q(5,1) = DQ15                                                  
         Q(5,2) = DQ25                                                  
         Q(5,3) = Q1(5,3)                                               
         Q(5,4) = Q1(5,4)                                               
         Q(5,5) = Q1(5,5)                                               
         RETURN                                                         
      ENDIF                                                             
                                                                        
      DQ13 = Q1(3,1) + D(1,4)*Q1(4,3) + D(1,5)*Q1(5,3)                  
      DQ23 = Q1(3,2) + D(2,4)*Q1(4,3) + D(2,5)*Q1(5,3)                  
                                                                        
      IF (IAPROX.GE.3) THEN                                             
         DQ13 = DQ13 + D(1,3)*Q1(3,3)                                   
         DQ14 = DQ14 + D(1,3)*Q1(4,3)                                   
         DQ15 = DQ15 + D(1,3)*Q1(5,3)                                   
         DQ23 = DQ23 + D(2,3)*Q1(3,3)                                   
         DQ24 = DQ24 + D(2,3)*Q1(4,3)                                   
         DQ25 = DQ25 + D(2,3)*Q1(5,3)                                   
         Q(1,1) = Q1(1,1) + D(1,4)*(Q1(4,1)+DQ14)                       
     &                    + D(1,5)*(Q1(5,1)+DQ15)                       
     &                    + D(1,3)*(Q1(3,1)+DQ13)                       
         Q(2,1) = Q1(2,1) + D(1,4)*Q1(4,2) + D(1,5)*Q1(5,2)             
     &                    + D(2,4)*DQ14    + D(2,5)*DQ15                
     &                    + D(1,3)*Q1(3,2) + D(2,3)*DQ13                
         Q(2,2) = Q1(2,2) + D(2,4)*(Q1(4,2)+DQ24)                       
     &                    + D(2,5)*(Q1(5,2)+DQ25)                       
     &                    + D(2,3)*(Q1(3,2)+DQ23)                       
         Q(4,1) = DQ14                                                  
         Q(4,2) = DQ24                                                  
      ENDIF                                                             
                                                                        
      DQ53 = Q1(5,3) + D(5,3)*Q1(3,3) + D(5,4)*Q1(4,3)                  
      DQ54 = Q1(5,4) + D(5,3)*Q1(4,3) + D(5,4)*Q1(4,4)                  
                                                                        
      Q(3,1) = DQ13                                                     
      Q(3,2) = DQ23                                                     
                                                                        
      Q(5,1) =    DQ13*D(5,3) +    DQ14*D(5,4) +  DQ15                  
      Q(5,2) =    DQ23*D(5,3) +    DQ24*D(5,4) +  DQ25                  
      Q(5,3) = Q1(3,3)*D(5,3) + Q1(4,3)*D(5,4) +  Q1(5,3)               
      Q(5,4) = Q1(4,3)*D(5,3) + Q1(4,4)*D(5,4) +  Q1(5,4)               
      Q(5,5) = Q1(5,5) + (Q1(5,3)+DQ53)*D(5,3) + (Q1(5,4)+DQ54)*D(5,4)  
                                                                        
      RETURN                                                            
      END                                                               
*