appdots.f

Go to the documentation of this file.

!                   ******************
                    SUBROUTINE appdots
!                   ******************
!
     &( x , mesh )
!
!***********************************************************************
! BIEF   V6P2                                   11/05/2012
!***********************************************************************
!
!brief    Assemble and Prints P_DOTS of a vector which is not assembled.
!         This is for debugging parallelism and applies e.g. to
!         diagonals and right-hand sides of linear systems
!         P_DOTS is the dot product of a vector in parallel.
!         Only the master processor (IPID=0) prints the result
!
!history  J-M HERVOUET (LNHE)
!+        11/05/2012
!+        V6P2
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| X              |-->| THE VECTOR
!| MESH           |-->| THE MESH STRUCTURE
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      TYPE(bief_obj), INTENT(IN) :: X
      TYPE(bief_mesh), INTENT(INOUT) :: MESH
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DOUBLE PRECISION RESULT
!
!-----------------------------------------------------------------------
!
      CALL os('X=Y     ',x=mesh%T,y=x)
!
      IF(ncsize.GT.1) CALL parcom(mesh%T,2,mesh)
!
      result=p_dots(mesh%T,mesh%T,mesh)
!
      IF(ipid.EQ.0) THEN
        WRITE(lu,*) x%NAME,' DOTS=',result
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END