# invert.f

Go to the documentation of this file.
```00001 C:\opentelemac\v7p0\sources\tomawac\invert.f
00002 !
00047                         SUBROUTINE INVERT
00048 !                       *****************
00049 !
00050      &( RN    , N     , NP    )
00051 !
00052 !***********************************************************************
00053 ! TOMAWAC   V6P2                                   25/06/2012
00054 !***********************************************************************
00055 !
00056 !
00057 !
00058 !
00059 !
00060 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00061 !| N              |-->| MATRIX DIMENSION
00062 !| NP             |-->| MATRIX PHYSICAL DIMENSION
00063 !| RN             |<->| MATRIX TO INVERT
00064 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00065 !
00066       IMPLICIT NONE
00067 !
00068 !.....VARIABLES IN ARGUMENT
00069 !     """"""""""""""""""""
00070       INTEGER N, NP
00071       DOUBLE PRECISION RN(NP,NP)
00072 !
00073 !.....LOCAL VARIABLES
00074 !     """""""""""""""
00075       INTEGER INDX(NP),I,J
00076       DOUBLE PRECISION Y(NP,NP), A(NP,NP)
00077 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00078 !
00079       A=RN
00080       DO I=1,N
00081         DO J=1,N
00082           Y(I,J)=0.
00083         ENDDO
00084         Y(I,I)=1.
00085       ENDDO
00086       CALL LUDCMP(A,N,NP,INDX)
00087       DO J=1,N
00088         CALL LUBKSB(A,N,NP,INDX,Y(1,J))
00089       ENDDO
00090       RN=Y
00091 !
00092       RETURN
00093       END
```