parcom2_comp.f

Go to the documentation of this file.

!                   ***********************
                    SUBROUTINE parcom2_comp
!                   ***********************
!
     &( x1 , x2 , x3 ,errx, npoin , nplan , icom , ian , mesh )
!
!***********************************************************************
! BIEF   V6P1                                   24/02/2016
!***********************************************************************
!
!brief    COMPLEMENTS A VECTOR AND ERROR VECTOR AT THE INTERFACES BETWEEN
!+                SUB-DOMAINS.
!+
!+            X CAN BE A BLOCK OF VECTORS. IN THIS CASE, ALL THE
!+                VECTORS IN THE BLOCK ARE TREATED.
!
!note     IMPORTANT : FROM RELEASE 5.9 ON, IDENTICAL VALUES ARE
!+                     ENSURED AT INTERFACE POINTS SO THAT DIFFERENT
!+                     PROCESSORS WILL ALWAYS MAKE THE SAME DECISION
!+                     IN TESTS ON REAL NUMBERS.
!
!warning  IF THE VECTORS HAVE A SECOND DIMENSION, IT IS
!+            IGNORED FOR THE TIME BEING
!
!history  R.NHEILI (Univerte de Perpignan, DALI)
!+        24/02/2016
!+        V7P3
!+
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IAN            |-->| NUMBER OF VECTORS TO BE CONDIDERED (1, 2 OR 3)
!| ICOM           |-->| COMMUNICATION MODE
!|                |   | = 1 : VALUE WITH MAXIMUM ABSOLUTE VALUE TAKEN
!|                |   | = 2 : CONTRIBUTIONS ADDED
!|                |   | = 3 : MAXIMUM CONTRIBUTION RETAINED
!|                |   | = 4 : MINIMUM CONTRIBUTION RETAINED
!| MESH           |-->| MESH STRUCTURE
!| NPLAN          |-->| NUMBER OF PLANES IN THE 3D MESH OF PRISMS
!| NPOIN          |-->| NUMBER OF POINTS
!| X1             |<->| VECTOR TO BE COMPLETED
!| X2             |<->| VECTOR TO BE COMPLETED
!| X3             |<->| VECTOR TO BE COMPLETED
!| ERRX           |<->| ERRORS VECTOR TO BE COMPLETED
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_parcom2_comp => parcom2_comp
!
      USE declarations_special
      IMPLICIT NONE
!
!-----------------------------------------------------------------------
!
      INTEGER, INTENT(IN) :: ICOM,NPOIN,NPLAN,IAN
!
!     STRUCTURES: VECTORS OR BLOCKS
!
      TYPE(bief_mesh) , INTENT(INOUT) :: MESH
      DOUBLE PRECISION, INTENT(INOUT) :: X1(npoin,nplan)
      DOUBLE PRECISION, INTENT(INOUT) :: X2(npoin,nplan)
      DOUBLE PRECISION, INTENT(INOUT) :: X3(npoin,nplan)
      DOUBLE PRECISION, INTENT(INOUT) :: ERRX(*)
!
!-----------------------------------------------------------------------
!
      CALL paraco_comp(x1,x2,x3,errx,npoin,icom,ian,nplan,
     &            mesh%NB_NEIGHB,mesh%NB_NEIGHB_PT%I,mesh%LIST_SEND%I,
     &            mesh%NH_COM%I,mesh%NH_COM%DIM1,mesh%BUF_SEND%R,
     &            mesh%BUF_RECV%R,mesh%BUF_SEND%DIM1,
     &            mesh%BUF_SEND_ERR%R,mesh%BUF_RECV_ERR%R)
!
!-----------------------------------------------------------------------
!
      RETURN
      END
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!