Télécharger hmc1o2.eso

Retour à la liste

Numérotation des lignes :

  1. C HMC1O2 SOURCE CHAT 05/01/13 00:24:50 5004
  2.  
  3. SUBROUTINE HESSMAC1O2 (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),auxvecm(8),dx,typx(ndimx)
  8. do i=1,8
  9. auxx(i)=0.D0
  10. auxvecm(i)=0.D0
  11. enddo
  12. do i=1,ndimx
  13. auxx(i)=x(i)
  14. enddo
  15. do i=1,ndimx
  16. dx=deltax*max(abs(x(i)),typx(i))
  17. auxx(i)=x(i)+dx
  18. call vecflMAC(auxx,ndimx,Amat(1,i),nmodel)
  19. auxx(i)=x(i)-dx
  20. call vecflMAC(auxx,ndimx,auxvecm,nmodel)
  21. auxx(i)=x(i)
  22. do j=1,ndimx
  23. Amat(j,i)=(Amat(j,i)-auxvecm(j))/(2.D0*dx)
  24. enddo
  25. enddo
  26. return
  27. end
  28.  
  29.  
  30.  
  31.  

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