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	utke	1.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