Esope source files

   1 : C INVER2    SOURCE    CHAT      05/01/13    00:42:40     5004           
   2 :       SUBROUTINE INVER2(A,B,DET,ICRI)
   3 : C     CALCUL DE L'INVERSE DE LA MATRICE A
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8 (A-H,O-Z)  
   6 :       DIMENSION A(2,2),B(2,2)
   7 :       DATA EPSI/1.D-9/
   8 :       ICRI=0
   9 :       DET=A(1,1)*A(2,2)-A(2,1)*A(1,2)
  10 :       IF(ABS(DET).LE.EPSI) THEN
  11 :         ICRI=1
  12 :         RETURN
  13 :       END IF
  14 :       B(1,1)= A(2,2)/DET
  15 :       B(1,2)=-A(1,2)/DET
  16 :       B(2,1)=-A(2,1)/DET
  17 :       B(2,2)= A(1,1)/DET
  18 :       END
  19 :  
  20 :