anti_dirac.f

Go to the documentation of this file.

!                   ************************************
                    DOUBLE PRECISION FUNCTION anti_dirac
!                   ************************************
!
     &( x,eps )
!
!***********************************************************************
! TELEMAC2D   V7P1                                    29/09/2015
!***********************************************************************
!
!brief    GIVES 0 IF (ABS(X)<EPS) AND 1 ELSWHERE
!+
!
!history  R.ATA (LNHE)
!+        29/09/2015
!+        V7P1
!+
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| X              |-->| INPUT VARIABLE
!| EPS            |-->| TOLERENCE
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_anti_dirac => anti_dirac
      USE declarations_special
!
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DOUBLE PRECISION, INTENT(IN) :: X,EPS
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      IF(abs(x).LT.eps) THEN
        anti_dirac=0.d0
      ELSE
        anti_dirac=1.0d0
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END