Esope source files

   1 : C HINV22    SOURCE    AM        11/05/26    21:15:37     6982           
   2 :       subroutine HINV22 (H,Hinv)
   3 : C     Routine permettant d'inverser une matrice 2,2
   4 : 
   5 :         IMPLICIT INTEGER(I-N)
   6 :         IMPLICIT REAL*8(A-H,O-Z)
   7 :         DIMENSION H(2,2), Hinv(2,2)
   8 : *
   9 :         detH=H(1,1)*H(2,2)-H(1,2)*H(2,1)
  10 :         Hinv(1,1)=H(2,2)
  11 :         Hinv(2,2)=H(1,1)
  12 :         Hinv(1,2)=-H(1,2)
  13 :         Hinv(2,1)=-H(2,1)
  14 :        DO I=1,2
  15 :         DO J=1,2
  16 :          Hinv(I,J)=1/detH*Hinv(I,J)
  17 :         ENDDO
  18 :        ENDDO 
  19 :         return
  20 :         end
  21 :