OpenAD/code_ad_moc/f2sDefs.F

       Subroutine oad_s_LAGRAN(i,x,a,sp,lagranRet)
#include "CPP_OPTIONS.h"
        INTEGER i
        _RS x
        _RL a(4)
        _RL lagranRet
        INTEGER sp

C-      local variables:
        INTEGER k
        _RL numer,denom

        numer = 1. _d 0
        denom = 1. _d 0

#ifdef TARGET_NEC_SX
!CDIR UNROLL=8
#endif /* TARGET_NEC_SX */
        do k=1,sp
         if ( k .ne. i) then
          denom = denom*(a(i) - a(k))
          numer = numer*(x    - a(k))
         endif
        enddo

        lagranRet = numer/denom

       END

1	Subroutine oad_s_LAGRAN(i,x,a,sp,lagranRet)
2	#include "CPP_OPTIONS.h"
3	INTEGER i
4	_RS x
5	_RL a(4)
6	_RL lagranRet
7	INTEGER sp
8
9	C- local variables:
10	INTEGER k
11	_RL numer,denom
12
13	numer = 1. _d 0
14	denom = 1. _d 0
15
16	#ifdef TARGET_NEC_SX
17	!CDIR UNROLL=8
18	#endif /* TARGET_NEC_SX */
19	do k=1,sp
20	if ( k .ne. i) then
21	denom = denom*(a(i) - a(k))
22	numer = numer*(x - a(k))
23	endif
24	enddo
25
26	lagranRet = numer/denom
27
28	END
29