Télécharger hmc1o1.eso

Retour à la liste

Numérotation des lignes :

hmc1o1
  1. C HMC1O1 SOURCE CHAT 05/01/13 00:24:46 5004
  2.  
  3. SUBROUTINE HESSMAC1O1 (X,NDIMX,AMAT,NDIMA,NMODEL,DELTAX,TYPX)
  4. IMPLICIT INTEGER(I-N)
  5. integer ndimx,ndima,nmodel,i,j
  6. real*8 x(ndimx),Amat(ndima,ndima),deltax
  7. real*8 auxx(8),vecm(8),dx,typx(ndimx)
  8. do i=1,8
  9. auxx(i)=0.D0
  10. vecm(i)=0.D0
  11. enddo
  12. call vecflMAC(x,ndimx,vecm,nmodel)
  13. do i=1,ndimx
  14. auxx(i)=x(i)
  15. enddo
  16. do i=1,ndimx
  17. dx=sign(deltax*max(abs(x(i)),typx(i)),x(i))
  18. auxx(i)=x(i)+dx
  19. call vecflMAC(auxx,ndimx,Amat(1,i),nmodel)
  20. auxx(i)=x(i)
  21. do j=1,ndimx
  22. Amat(j,i)=(Amat(j,i)-vecm(j))/dx
  23. enddo
  24. enddo
  25. return
  26. end
  27.  
  28.  
  29.  
  30.  

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