ident.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	ident (x, px, n, id)

Function/Subroutine Documentation

ident()

subroutine ident	(	real*8, dimension(n), intent(in)	x,
		real*8, intent(in)	px,
		integer, intent(in)	n,
		integer, intent(out)	id
	)

      implicit none
!
      integer n,id,n2,m
!
      real*8 x(n),px
!
      intent(in) x,px,n
!
      intent(out) id
!
      id=0
      if(n.eq.0) return
      n2=n+1
      do                                                           
         m=(n2+id)/2
         if(px.ge.x(m)) then
            id=m   
         else
            n2=m    
         endif
         if((n2-id).eq.1) return
      enddo