Esope source files

   1 : C INVERT    SOURCE    CHAT      05/01/13    00:42:55     5004           
   2 :       SUBROUTINE INVERT(A,B,DET,ICRI)
   3 : C     CALCULDE L'INVERSE DE LA MATRICE A
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8 (A-H,O-Z)  
   6 :       DIMENSION A(3,3),B(3,3)
   7 :       DATA EPSI/1.D-9/
   8 :       ICRI=0
   8 :       CALL DETERM(A,DET)
   9 :       IF(DABS(DET).LE.EPSI) THEN
  10 :         ICRI=1
  11 :         RETURN
  12 :       END IF
  13 :       B(1,1)= A(2,2)*A(3,3)-A(3,2)*A(2,3)
  14 :       B(1,2)=-A(1,2)*A(3,3)+A(3,2)*A(1,3)
  15 :       B(1,3)= A(1,2)*A(2,3)-A(2,2)*A(1,3)
  16 :       B(2,1)=-A(2,1)*A(3,3)+A(3,1)*A(2,3)
  17 :       B(2,2)= A(1,1)*A(3,3)-A(3,1)*A(1,3)
  18 :       B(2,3)=-A(1,1)*A(2,3)+A(2,1)*A(1,3)
  19 :       B(3,1)= A(2,1)*A(3,2)-A(2,2)*A(3,1)
  20 :       B(3,2)=-A(1,1)*A(3,2)+A(3,1)*A(1,2)
  21 :       B(3,3)= A(1,1)*A(2,2)-A(2,1)*A(1,2)
  22 :       DO 10 I=1,3
  23 :       DO 10 J=1,3
  24 :   10  B(I,J)=B(I,J)/DET
  25 :       END
  26 :  
  27 :