Sources Esope | Cast3M

hinv22
C HINV22    SOURCE    AM        11/05/26    21:15:37     6982
      subroutine HINV22 (H,Hinv)
C     Routine permettant d'inverser une matrice 2,2
 
        IMPLICIT INTEGER(I-N)
        IMPLICIT REAL*8(A-H,O-Z)
        DIMENSION H(2,2), Hinv(2,2)
*
        detH=H(1,1)*H(2,2)-H(1,2)*H(2,1)
        Hinv(1,1)=H(2,2)
        Hinv(2,2)=H(1,1)
        Hinv(1,2)=-H(1,2)
        Hinv(2,1)=-H(2,1)
       DO I=1,2
        DO J=1,2
         Hinv(I,J)=1/detH*Hinv(I,J)
        ENDDO
       ENDDO
        return
        end