Télécharger mucpr2.eso

Retour à la liste

Numérotation des lignes :

  1. C MUCPR2 SOURCE PV 17/06/16 14:33:49 9460
  2. subroutine mucpr2(inpl,indl,re,xbuffp,xbuffd)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. * appelé par mucpr1
  6. REAL*8 re(indl,inpl),xbuffp(inpl),xbuffd(indl)
  7. do 10 ind=1,indl
  8. * remplace par appel ddot car bug gfortran 7.1?
  9. * aux=0.d0
  10. * do 20 inp=1,inpl
  11. * aux=+aux+xbuffp(inp)*re(ind,inp)
  12. * 20 continue
  13. * xbuffd(ind)=aux
  14. xbuffd(ind)=ddot(inpl,xbuffp(1),1,re(ind,1),indl)
  15. 10 continue
  16. end
  17.  
  18.  
  19.  
  20.  
  21.  
  22.  

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