bief_dim2_ext.f

Go to the documentation of this file.

!                   ******************************
                    INTEGER FUNCTION bief_dim2_ext
!                   ******************************
!
     &(ielm1,ielm2,sto,typext,mesh)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    GIVES THE SECOND DIMENSION OF A MATRICE'S EXTRA-DIAGONAL
!+                TERMS.
!
!history  J-M HERVOUET (LNH)
!+        12/12/94
!+        V5P1
!+
!
!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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IELM1          |-->| TYPE OF LINE ELEMENT
!| IELM2          |-->| TYPE OF COLUMN ELEMENT
!| MESH           |-->| MESH STRUCTURE
!| STO            |-->| TYPE OF STORAGE
!| TYPEXT         |-->| TYPE OF OFF-DIAGONAL TERMS
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_bief_dim2_ext => bief_dim2_ext
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER         , INTENT(IN) :: IELM1,IELM2,STO
      CHARACTER(LEN=1), INTENT(IN) :: TYPEXT
      TYPE(bief_mesh) , INTENT(IN) :: MESH
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER NDIM
!
!-----------------------------------------------------------------------
!
      IF(typext.EQ.'0') THEN
!
        bief_dim2_ext = 1
!
      ELSEIF(sto.EQ.1) THEN
!
!       NDIM IS HERE THE SECOND DIMENSION OF VECTOR
        ndim =     bief_nbpel(ielm1,mesh)*bief_nbpel(ielm2,mesh)
     &       - min(bief_nbpel(ielm1,mesh),bief_nbpel(ielm2,mesh))
!       SYMMETRICAL MATRIX
        IF(ielm1.EQ.ielm2.AND.typext.EQ.'S') THEN
          ndim = ndim / 2
        ENDIF
        bief_dim2_ext = ndim
!
      ELSEIF(sto.EQ.3) THEN
!
!       ASSEMBLED EBE STORAGE OR EDGE BASED
        bief_dim2_ext = 1
!
      ELSE
!
        WRITE(lu,101) sto
101     FORMAT(1x,'BIEF_DIM2_EXT : UNKNOWN TYPE OF STORAGE: ',1i6)
        CALL plante(1)
        stop
!
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END