errmax.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE errmax
!                   *****************
!
     &(x1,x2,err,ierr)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    COMPUTES MAX DIFFERENCES BETWEEN 2 COMPUTED ARRAYS.
!
!history  E. BARROS
!+        27/04/93
!+
!+
!
!history  A. LEOPARDI (UNINA)
!+        02/10/00
!+        V5P1
!+   UPGRADE
!
!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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| ERR            |<--| MAXIMUM ABSOLUTE DIFFERENCE
!| IERR           |<--| POINT WHERE THE DIFFERENCE OCCURS
!| X1             |-->| ARRAY TO COMPARE WITH X2
!| X2             |-->| ARRAY TO COMPARE WITH X1
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief_def
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER          , INTENT(OUT) :: IERR
      DOUBLE PRECISION , INTENT(OUT) :: ERR
      type(bief_obj)  , INTENT(IN)  :: x1,x2
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER I
!
      INTRINSIC abs
!
!---------------------------------------------------------------------
!
      ierr=1
      err=-1.d0
      DO i=1,x1%DIM1
!
        IF(abs(x1%R(i)-x2%R(i)).GT.err) THEN
          err=abs(x1%R(i)-x2%R(i))
          ierr=i
        ENDIF
!
      ENDDO
!
!---------------------------------------------------------------------
!
      RETURN
      END