Esope source files

   1 : C MUCPLI    SOURCE    CHAT      09/10/09    21:21:17     6519           
   2 :       subroutine mucpli(lmchpo1,mrigid,lmchpo)
   3 : C                                                                               
   4 : C  ****  multiplication d'une matrice(mrigid) par une suite de champoin (lmchpo1)
   5 : C  ****  le resultat est une liste de champoin (lmchpoi).
   6 : C  ****  iret=ire2*ire1                                                         
   7 : C  ****  les champpoint resultat ont la meme dimension que la matrice.             
   8 : C          
   9 :       IMPLICIT INTEGER(I-N)
  10 :       IMPLICIT REAL*8 (A-H,O-Z)
  11 : -INC CCOPTIO
  12 : -INC SMLCHPO
  13 : -INC SMCHPOI
  14 : -INC SMRIGID
  15 : -INC SMCOORD
  16 : 
  17 :       MLCHP1 = lmchpo1 
  18 :       segact mlchp1
  19 :       n1 = mlchp1.ichpoi(/1)
  20 : 
  20 :       if (n1.le.0) call erreur(3)
  21 :       if (ierr.ne.0) return
  22 : 
  23 :       mchpo1 = mlchp1.ichpoi(1)
  24 :       segini mlchpo
  25 :       lmchpo = mlchpo
  25 :       call mucpri(mchpo1,mrigid,mchpoi)
  26 :       if (ierr.ne.0) return
  27 : 
  28 :       ichpoi(1) = mchpoi
  29 : 
  30 :       do ic = 2,n1
  31 :         mchpo1 = mlchp1.ichpoi(ic)
  31 :         call mucpr1(mchpo1,mrigid,mchpoi)
  32 :         if (ierr.ne.0) return
  33 :         ichpoi(ic) = mchpoi
  34 :       enddo
  35 : 
  36 :       segdes mlchp1, mlchpo
  37 :       return
  38 :       end                                                               
  39 :  
  40 :