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