Télécharger hinv33.eso

Retour à la liste

Numérotation des lignes :

  1. C HINV33 SOURCE AM 11/05/26 21:15:38 6982
  2. subroutine HINV33 (H,Hinv)
  3. C Routine permettant d'inverser une matrice 3,3
  4.  
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. DIMENSION H(3,3), Hinv(3,3), comatH(3,3),TcomatH(3,3)
  8. *
  9. detH=H(1,1)*H(2,2)*H(3,3)+H(2,1)*H(3,2)*H(1,3)
  10. & +H(3,1)*H(1,2)*H(2,3)-H(3,1)*H(2,2)*H(1,3)
  11. & -H(3,2)*H(2,3)*H(1,1)-H(3,3)*H(2,1)*H(1,2)
  12.  
  13. comatH(1,1)=H(2,2)*H(3,3)-H(3,2)*H(2,3)
  14. comatH(2,1)=-(H(1,2)*H(3,3)-H(3,2)*H(1,3))
  15. comatH(3,1)=H(1,2)*H(2,3)-H(2,2)*H(1,3)
  16. comatH(1,2)=-(H(2,1)*H(3,3)-H(3,1)*H(2,3))
  17. comatH(2,2)=H(1,1)*H(3,3)-H(3,1)*H(1,3)
  18. comatH(3,2)=-(H(1,1)*H(2,3)-H(2,1)*H(1,3))
  19. comatH(1,3)=H(2,1)*H(3,2)-H(3,1)*H(2,2)
  20. comatH(2,3)=-(H(1,1)*H(3,2)-H(3,1)*H(1,2))
  21. comatH(3,3)=H(1,1)*H(2,2)-H(2,1)*H(1,2)
  22.  
  23. call TRANSP (3,3,comatH,TcomatH)
  24. DO I=1,3
  25. DO J=1,3
  26. Hinv(I,J)=1/detH*TcomatH(I,J)
  27. ENDDO
  28. ENDDO
  29. return
  30. end
  31.  
  32.  

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