5 * Revision 1.1.1.1 1996/04/01 15:02:26 mclareni
10 FUNCTION LKKSPL(X,T,M)
12 #include "gen/imp64.inc"
15 ************************************************************************
16 * NORBAS, VERSION: 15.03.1993
17 ************************************************************************
19 * FOR GIVEN X AND A NON-DECREASING SEQUENCE T(1) , T(2) ,..., T(M)
20 * THE FUNCTION LKKSPL COMPUTES AND RETURNS THE UNIQUELY DETERMINED
21 * INDEX KK SUCH THAT THE RELATION
22 * T(KK-1) <= X < T(KK)
24 * IT RETURNS KK = 1 , IF X < T(1) , AND
25 * KK = M+1 , IF X >= T(M) .
27 * (THE SEARCH ALGORITHM IS OF REGULA FALSI - TYPE).
29 ************************************************************************
33 ELSEIF(X .LT. T(1)) THEN
40 10 M3=M1+(M2-M1)*(X-T1)/(T2-T1)
44 IF((T3 .GT. X .OR. X .GE. T4) .AND.
45 + (T3 .NE. X .OR. X .NE. T4)) THEN