Esope source files

   1 : C MUCPR2    SOURCE    CHAT      05/01/13    01:54:55     5004           
   2 :       subroutine mucpr2(inpl,indl,re,xbuffp,xbuffd,itypel)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)
   5 : *  appelé par mucpr1
   6 : *  si itypel=22 on impose l'egalite des multiplicateurs en supprimant
   7 : *  Juillet 2003 ==> simple multiplicateur de Lagrange
   8 : *  leur matrice de couplage
   9 :       dimension re(indl,inpl),xbuffp(inpl),xbuffd(indl)
  10 :       do 10 ind=1,indl
  11 :        aux=0.d0
  12 :        inpd=1
  13 : *      if (itypel.eq.22.and.ind.le.1) inpd=2
  14 :        do 20 inp=inpd,inpl
  15 :         aux=aux+re(ind,inp)*xbuffp(inp)
  16 :   20   continue
  17 :        xbuffd(ind)=aux
  18 :   10  continue
  19 :       end
  20 :  
  21 :  
  22 :