Télécharger mucpr2.eso

Retour à la liste

Numérotation des lignes :

  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.  
  23.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales