cg1113.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE cg1113
!                   *****************
!
     &(x,dim1,dim2,ikle,nelem,nelmax)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    CHANGES THE DISCRETISATION OF A VECTOR
!+                FROM 11 TO 13 HERE (LINEAR TO QUADRATIC).
!
!history  J-M HERVOUET (LNH)
!+        06/02/08
!+        V5P9
!+
!
!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,6)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER IELEM,IDIM
!
!-----------------------------------------------------------------------
!
!VOCL LOOP,NOVREC
!DIR$ IVDEP
      DO idim  = 1 , dim2
      DO ielem = 1 , nelem
!
        x(ikle(ielem,4),idim) = 0.5d0 * ( x(ikle(ielem,1),idim)
     &                                  + x(ikle(ielem,2),idim) )
        x(ikle(ielem,5),idim) = 0.5d0 * ( x(ikle(ielem,2),idim)
     &                                  + x(ikle(ielem,3),idim) )
        x(ikle(ielem,6),idim) = 0.5d0 * ( x(ikle(ielem,3),idim)
     &                                  + x(ikle(ielem,1),idim) )
!
      ENDDO
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END