Esope source files

   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 :         
  22 :         call TRANSP (3,3,comatH,TcomatH)
  23 :        DO I=1,3
  24 :         DO J=1,3
  25 :          Hinv(I,J)=1/detH*TcomatH(I,J)
  26 :         ENDDO
  27 :        ENDDO
  28 :         return
  29 :         end
  30 :