asmvec.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE asmvec
!                   *****************
!
     &(x, ikle,npoin,nelem,nelmax,ndp,w,init,lv)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    MULTIPLICATIVE ASSEMBLY FOR A VECTOR.
!
!warning  THIS VECTOR IS INITIALISED TO 1 IF INIT = .TRUE.
!
!history  J-M HERVOUET (LNH)    ; F  LEPEINTRE (LNH)
!+        17/08/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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IKLE           |-->| CONNECTIVITY TABLE
!| INIT           |-->| LOGICAL : IF TRUE X IS INITIALISED TO 0.
!| LV             |-->| VECTOR LENGTH OF THE COMPUTER
!| NDP            |-->| SECOND DIMENSION OF IKLE
!| NELEM          |-->| NUMBER OF ELEMENTS IN THE MESH
!| NELMAX         |-->| FIRST DIMENSION OF IKLE AND W.
!| NPOIN          |-->| NUMBER OF POINTS IN X
!| W              |-->| WORK ARRAY WITH A NON ASSEMBLED FORM OF THE
!|                |   | RESULT
!|                |   | W HAS DIMENSION NELMAX * NDP(IELM)
!|                |   | NDP IS THE NUMBER OF POINTS IN THE ELEMENT
!| X              |<->| ASSEMBLED VECTOR
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_asmvec => asmvec
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER         , INTENT(IN)    :: NELMAX,NPOIN,NELEM,NDP,LV
      DOUBLE PRECISION, INTENT(INOUT) :: X(npoin)
      INTEGER         , INTENT(IN)    :: IKLE(nelmax,ndp)
      DOUBLE PRECISION, INTENT(IN)    :: W(nelmax,ndp)
      LOGICAL         , INTENT(IN)    :: INIT
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER IDP
!
!-----------------------------------------------------------------------
!   INITIALISES VECTOR X TO 1 IF(INIT)
!-----------------------------------------------------------------------
!
      IF(init) CALL ov('X=C     ', x=x, c=1.d0, dim1=npoin)
!
!-----------------------------------------------------------------------
!   ASSEMBLES, CONTRIBUTION OF LOCAL POINTS 1,... TO NDP
!-----------------------------------------------------------------------
!
      DO idp = 1 , ndp
!
        CALL asmve1(x, ikle(1,idp),w(1,idp),npoin,nelem,nelmax,lv)
!
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END