appdots.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\appdots.f
00002 !
00048                      SUBROUTINE APPDOTS
00049 !                    ******************
00050 !
00051      &( X , MESH )
00052 !
00053 !***********************************************************************
00054 ! BIEF   V6P2                                   11/05/2012
00055 !***********************************************************************
00056 !
00057 !         This is for debugging parallelism and applies e.g. to
00058 !         diagonals and right-hand sides of linear systems
00059 !         P_DOTS is the dot product of a vector in parallel.
00060 !         Only the master processor (IPID=0) prints the result
00061 !
00062 !
00063 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00064 !| X              |-->| THE VECTOR
00065 !| MESH           |-->| THE MESH STRUCTURE
00066 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00067 !
00068       USE BIEF
00069 !
00070       IMPLICIT NONE
00071       INTEGER LNG,LU
00072       COMMON/INFO/LNG,LU
00073 !
00074 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00075 !
00076       TYPE(BIEF_OBJ), INTENT(IN) :: X
00077       TYPE(BIEF_MESH), INTENT(INOUT) :: MESH
00078 !
00079 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00080 !
00081       DOUBLE PRECISION RESULT
00082 !
00083 !-----------------------------------------------------------------------
00084 !
00085       CALL OS('X=Y     ',X=MESH%T,Y=X)
00086 !
00087       IF(NCSIZE.GT.1) CALL PARCOM(MESH%T,2,MESH)
00088 !
00089       RESULT=P_DOTS(MESH%T,MESH%T,MESH)
00090 !
00091       IF(IPID.EQ.0) THEN
00092         WRITE(LU,*) X%NAME,' DOTS=',RESULT
00093       ENDIF
00094 !
00095 !-----------------------------------------------------------------------
00096 !
00097       RETURN
00098       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0