gettri.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE gettri
!                   *****************
!
     &(m,mdiff,teta,mesh3d,nplan,npoin2,nseg2d,ielm3,nelem2)
!
!***********************************************************************
! BIEF   V6P2                                   21/08/2010
!***********************************************************************
!
!brief    GETS THE TRIDIAGONAL PART OF A DIFFUSION MATRIX ON
!+                PRISMS AND REMOVES IT FROM THE INITIAL MATRIX.
!code
!+            IF MTRI IS THIS TRIDIAGONAL PART, M THE RESULT AND MDIF
!+            THE DIFFUSION MATRIX, THIS SUBROUTINE DOES:
!+
!+            M = TETA * MTRI
!+            MDIF CHANGED INTO (1-TETA) * MDIF
!
!warning  THE JACOBIAN MUST BE POSITIVE
!
!history  J-M HERVOUET (LNH)
!+        16/06/05
!+        V5P6
!+
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IELM3          |<--| TYPE OF ELEMENT
!| M              |<--| TRIDIAGONAL MATRIX
!| MDIFF          |-->| ORIGINAL DIFFUSION MATRIX
!| MESH3D         |-->| 3D MESH STRUCTURE
!| NELEM2         |-->| NUMBER OF TRIANGLES OF ORIGINAL 2D MESH
!| NPLAN          |-->| NUMBER OF PLANES
!| NPOIN2         |-->| NUMBER OF POINTS OF 2D MESH
!| NSEG2D         |-->| NUMBER OF SEGMENTS IN 2D
!| TETA           |-->| COEFFICIENT USED IN THE RESULT
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_gettri => gettri
!
      USE declarations_special
      IMPLICIT NONE
!
!-----------------------------------------------------------------------
!
      INTEGER, INTENT(IN) :: NPLAN,NPOIN2,NSEG2D,IELM3,NELEM2
!
      DOUBLE PRECISION, INTENT(IN)    :: TETA
      DOUBLE PRECISION, INTENT(INOUT) :: M(npoin2*nplan,*)
!
      TYPE(bief_obj),  INTENT(INOUT)  :: MDIFF
      TYPE(bief_mesh), INTENT(INOUT)  :: MESH3D
!
!-----------------------------------------------------------------------
!
      IF(mdiff%STO.EQ.1) THEN
!
        CALL gettriebe(m,mdiff%D%R,mdiff%X%R,teta,
     &                 mesh3d%IKLE%I,mesh3d%NPOIN,mesh3d%NELEM,
     &                 mesh3d%NELMAX,mesh3d,ielm3,nelem2,nplan,
     &                 mesh3d%KNOLG%I)
!
      ELSEIF(mdiff%STO.EQ.3) THEN
!
        CALL gettriseg(m,mdiff%D%R,mdiff%X%R,teta,
     &                 mesh3d%NPOIN,mesh3d,
     &                 mesh3d%NSEG,nseg2d,nplan,npoin2,ielm3)
!
      ELSE
!
        WRITE(lu,*) 'UNKNOWN STORAGE FOR MDIFF IN GETTRI'
        CALL plante(1)
        stop
!
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END