C GRACO9    SOURCE    PV        09/11/05    21:15:19     6535           
      SUBROUTINE GRACO9(ith)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)                                                  
*                                                                               
*  extrait de graco7 pour pouvoir //
*                                                                               
-INC SMMATRI                                                                    
-INC SMVECTD
-INC CCOPTIO                                                                    
      common/cgraco9c/milign,mvect1,mvect2,inc,nbthr
      ithr=ith-1
      ino=ilign(/1)
      DO 2 I=1+ithr,INO,nbthr
      LLIGN=ILIGN(I)                                                            
      NA=IMMMM(/1)                   
      DO 121 JPA=1,IMMMM(/1)  
      P=0.D0 
      IPP=IPPO(JPA)                                       
      KPA=IPPO(JPA+1)-IPP                                 
      idli=0
      do ii=1,kpa
       idli=max(idli,linc(ippo(jpa)+ii))
      enddo
      Q=mvect2.VECTBb(IDLI)
      DO 122 MPA=1,KPA          
      LL=LINC(MPA+IPP) 
      if (ll.ne.idli)
     >mvect1.VECTBb(LL+inc*ithr)=mvect1.VECTBb(LL+inc*ithr)+ 
     >  Q*XXVA(MPA+IPP)
      P = P + mvect2.VECTBb(LL)*XXVA(MPA+IPP)
122   CONTINUE  
      mvect1.VECTBb(IDLI+inc*ithr)=mvect1.vectbb(idli+inc*ithr)+p
121   CONTINUE                                                                  
  2   CONTINUE                                                                                                                                                                                
      RETURN                                                                    
      END