cg1112.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE cg1112
!                   *****************
!
     &(x,dim1,dim2,ikle,nelem,nelmax)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    CHANGES THE DISCRETISATION OF A VECTOR
!+                FROM 11 TO 12 HERE.
!
!history  J-M HERVOUET (LNH)
!+        09/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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| DIM1           |-->| FIRST DIMENSION OF X
!| DIM2           |-->| SECOND DIMENSION OF X
!| IKLE           |-->| CONNECTIVITY TABLE
!| NELEM          |-->| NUMBER OF ELEMENTS IN THE MESH
!| NELMAX         |-->| FIRST DIMENSION OF IKLE AND W.
!| X              |<--| VECTOR TO BE MODIFIED.
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER         , INTENT(IN)    :: NELEM,NELMAX,DIM1,DIM2
      DOUBLE PRECISION, INTENT(INOUT) :: X(dim1,dim2)
      INTEGER         , INTENT(IN)    :: IKLE(nelmax,4)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER IELEM,IDIM1
!
      DOUBLE PRECISION TIERS
!
!-----------------------------------------------------------------------
!
      tiers = 1.d0/3.d0
!
!-----------------------------------------------------------------------
!
!VOCL LOOP,NOVREC
!DIR$ IVDEP
      DO idim1  = 1 , dim2
      DO ielem = 1 , nelem
!
        x(ikle(ielem,4),idim1) = tiers * ( x(ikle(ielem,1),idim1)
     &                                  + x(ikle(ielem,2),idim1)
     &                                  + x(ikle(ielem,3),idim1) )
!
      ENDDO ! IELEM
      ENDDO ! IDIM1
!
!-----------------------------------------------------------------------
!
      RETURN
      END