mt01oo.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE mt01oo
!                   *****************
!
     &(a11,a12,a22,xmul,lgseg,nelem,nelmax)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    BUILDS THE MASS MATRIX FOR P1 SEGMENTS.
!code
!+    COMPUTES THE COEFFICIENTS OF THE FOLLOWING MATRIX:
!+
!+                              /
!+                    A    =   /  (P *P )*J(X,Y) DX
!+                     I J    /L    I  J
!+
!+    BY ELEMENTARY CELL
!+
!+    J(X,Y) : JACOBIAN OF THE ISOPARAMETRIC TRANSFORMATION
!
!warning  THE JACOBIAN MUST BE POSITIVE
!
!history  J-M HERVOUET (LNH)
!+        06/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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| A11            |<--| ELEMENTS OF MATRIX
!| A12            |<--| ELEMENTS OF MATRIX
!| A22            |<--| ELEMENTS OF MATRIX
!| LGSEG          |-->| LENGTH OF SEGMENTS.
!| NELEM          |-->| NUMBER OF ELEMENTS
!| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS
!| XMUL           |-->| MULTIPLICATION FACTOR
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_mt01oo => mt01oo
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN) :: NELEM,NELMAX
!
      DOUBLE PRECISION, INTENT(INOUT) :: A11(nelmax),A12(nelmax)
      DOUBLE PRECISION, INTENT(INOUT) :: A22(nelmax)
      DOUBLE PRECISION, INTENT(IN)    :: LGSEG(nelmax)
!
      DOUBLE PRECISION, INTENT(IN) :: XMUL
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER IELEM
      DOUBLE PRECISION XSUR6,DET1
!
!-----------------------------------------------------------------------
!
      xsur6  = xmul/6.d0
!
!-----------------------------------------------------------------------
!
      DO ielem = 1 , nelem
!
        det1 = lgseg(ielem) * xsur6
        a11(ielem) = det1 + det1
        a12(ielem) = det1
        a22(ielem) = det1 + det1
!
      ENDDO ! IELEM
!
!-----------------------------------------------------------------------
!
      RETURN
      END