This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / neutron / tbsplt.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1995/10/24 10:21:59  cernlib
6 * Geant
7 *
8 *
9 #include "geant321/pilot.h"
10 *CMZ :  3.21/02 29/03/94  15.41.48  by  S.Giani
11 *-- Author :
12       SUBROUTINE TBSPLT(A,E,NP,Y)
13 C       THIS ROUTINE DETERMINES A CROSS SECTION AT A GIVEN
14 C       ENERGY FROM A CROSS SECTION VERSUS ENERGY TABLE USING
15 C       A TABLE SPLITTING METHOD
16       DIMENSION A(1)
17       SAVE
18       IPP=1
19       IF(E.LE.A(1))GO TO 40
20       IF(E.GE.A(2*NP-1))GO TO 50
21       INDXH=NP
22       INDXL=0
23    10 IF(INDXL+1.EQ.INDXH)GO TO 30
24       J=(INDXH+INDXL)/2
25       N=2*J-1
26       IF(E.LE.A(N))GO TO 20
27       INDXL=J
28       GO TO 10
29    20 INDXH=J
30       GO TO 10
31    30 N=2*INDXH-1
32       Y=A(N-1)+(E-A(N-2))*(A(N+1)-A(N-1))/(A(N)-A(N-2))
33       RETURN
34    40 Y=A(IPP+1)
35       RETURN
36    50 Y=A(2*NP)
37       RETURN
38       END