https://github.com/cran/fields
Tip revision: 6769ffc81115fbf0bf7d9c566cf7ac81be0049dc authored by Doug Nychka on 25 July 2005, 00:00:00 UTC
version 3.04
version 3.04
Tip revision: 6769ffc
ifind.f
INTEGER FUNCTION IFIND(X,XK,N)
C FIND I SUCH THAT XK(I) LE X LT XK(I+1)
C IFIND=0 IF X LT XK(1)
C IFIND=N IF X GT XK(N)
C J F MONAHAN JAN 1982 DEPT OF STAT, N C S U, RALEIGH, NC 27650
REAL*8 X,XK(N)
IFIND=0
IF(X.LT.XK(1)) RETURN
IFIND=N
IF(X.GE.XK(N)) RETURN
IL=1
IU=N
1 IF(IU-IL.LE.1) GO TO 4
I=(IU+IL)/2
IF(X-XK(I)) 2,5,3
2 IU=I
GO TO 1
3 IL=I
GO TO 1
4 IFIND=IL
RETURN
5 IFIND=I
RETURN
END