as3_1212_s.f

Go to the documentation of this file.

!                   *********************
                    SUBROUTINE as3_1212_s
!                   *********************
!
     &(xm,nseg11,nseg12,xmt,nelmax,nelem,eltseg1,eltseg2,eltseg3,
     &                                   eltseg4,eltseg5,eltseg6)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    ASSEMBLES MATRICES EXTRA-DIAGONAL TERMS
!+                IN THE CASE OF EDGE-BASED STORAGE.
!+
!+            CASE OF QUASIBUBBLE-QUASIBUBBLE ELEMENT
!+                AND SYMMETRICAL MATRIX.
!
!history  J-M HERVOUET (LNH)
!+        30/06/99
!+        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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| ELTSEG1        |-->| FIRST SEGMENT OF A TRIANGLE
!| ELTSEG2        |-->| SECOND SEGMENT OF A TRIANGLE
!| ELTSEG3        |-->| THIRD SEGMENT OF A TRIANGLE
!| ELTSEG4        |-->| FOURTH SEGMENT OF A QUADRATIC TRIANGLE
!| ELTSEG5        |-->| FIFTH SEGMENT OF A QUADRATIC TRIANGLE
!| ELTSEG6        |-->| SIXTH SEGMENT OF A QUADRATIC TRIANGLE
!| NELEM          |-->| NUMBER OF ELEMENTS IN THE MESH
!| NELMAX         |-->| FIRST DIMENSION OF IKLE AND W.
!| NSEG11         |-->| NUMBER OF SEGMENTS (HERE JOINING LINEAR POINTS)
!| NSEG12         |-->| NUMBER OF SEGMENTS (HERE QUADRATIC)
!| XM             |<--| ASSEMBLED OFF-DIAGONAL TERMS XA12,23,31
!| XMT            |-->| ELEMENT BY ELEMENT STORAGE OF MATRIX
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER         , INTENT(IN)    :: NELMAX,NELEM,NSEG11,NSEG12
      INTEGER         , INTENT(IN)    :: ELTSEG1(nelmax)
      INTEGER         , INTENT(IN)    :: ELTSEG2(nelmax)
      INTEGER         , INTENT(IN)    :: ELTSEG3(nelmax)
      INTEGER         , INTENT(IN)    :: ELTSEG4(nelmax)
      INTEGER         , INTENT(IN)    :: ELTSEG5(nelmax)
      INTEGER         , INTENT(IN)    :: ELTSEG6(nelmax)
      DOUBLE PRECISION, INTENT(INOUT) :: XMT(nelmax,*)
      DOUBLE PRECISION, INTENT(INOUT) :: XM(nseg12)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER ISEG,IELEM
!
!-----------------------------------------------------------------------
!
!  INITIALISES
!
      DO iseg = 1 , nseg11
        xm(iseg) = 0.d0
      ENDDO
!
!  ASSEMBLES PART P1
!
      DO ielem = 1,nelem
!       TERM 12
        xm(eltseg1(ielem)) = xm(eltseg1(ielem)) + xmt(ielem,1)
!       TERM 23
        xm(eltseg2(ielem)) = xm(eltseg2(ielem)) + xmt(ielem,4)
!       TERM 31
        xm(eltseg3(ielem)) = xm(eltseg3(ielem)) + xmt(ielem,2)
      ENDDO
!
!  ASSEMBLES QUASIBUBBLE PART
!
      DO ielem = 1,nelem
!       TERM 14
        xm(eltseg4(ielem)) = xmt(ielem,3)
!       TERM 24
        xm(eltseg5(ielem)) = xmt(ielem,5)
!       TERM 34
        xm(eltseg6(ielem)) = xmt(ielem,6)
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END