errmax.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\utils\bief\errmax.f
00002 !
00054                      SUBROUTINE ERRMAX
00055 !                    *****************
00056 !
00057      &(X1,X2,ERR,IERR)
00058 !
00059 !***********************************************************************
00060 ! BIEF   V6P1                                   21/08/2010
00061 !***********************************************************************
00062 !
00063 !
00064 !
00065 !
00066 !
00067 !
00068 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00069 !| ERR            |<--| MAXIMUM ABSOLUTE DIFFERENCE
00070 !| IERR           |<--| POINT WHERE THE DIFFERENCE OCCURS
00071 !| X1             |-->| ARRAY TO COMPARE WITH X2
00072 !| X2             |-->| ARRAY TO COMPARE WITH X1
00073 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00074 !
00075       USE BIEF_DEF
00076 !
00077       IMPLICIT NONE
00078       INTEGER LNG,LU
00079       COMMON/INFO/LNG,LU
00080 !
00081 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00082 !
00083       INTEGER          , INTENT(OUT) :: IERR
00084       DOUBLE PRECISION , INTENT(OUT) :: ERR
00085       TYPE (BIEF_OBJ)  , INTENT(IN)  :: X1,X2
00086 !
00087 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00088 !
00089       INTEGER I
00090 !
00091       INTRINSIC ABS
00092 !
00093 !---------------------------------------------------------------------
00094 !
00095       IERR=1
00096       ERR=-1.D0
00097       DO I=1,X1%DIM1
00098 !
00099         IF(ABS(X1%R(I)-X2%R(I)).GT.ERR) THEN
00100            ERR=ABS(X1%R(I)-X2%R(I))
00101            IERR=I
00102         ENDIF
00103 !
00104       ENDDO
00105 !
00106 !---------------------------------------------------------------------
00107 !
00108       RETURN
00109       END

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