poly_rr Function

  1. special_fun_modules.f90
  2. nrutil
  3. poly_rr

public function poly_rr(x, coeffs)

Arguments

Type IntentOptional AttributesName
real(kind=SP), intent(in) :: x
real(kind=SP), intent(in), DIMENSION(:):: coeffs

Return Value real(kind=SP)


Called by

proc~~poly_rr~~CalledByGraph proc~poly_rr poly_rr proc~poly_rrv poly_rrv proc~poly_rrv->proc~poly_rr interface~poly poly interface~poly->proc~poly_rr interface~poly->proc~poly_rrv proc~poly_msk_rrv poly_msk_rrv interface~poly->proc~poly_msk_rrv proc~poly_msk_rrv->proc~poly_rrv

Contents

Source Code

poly_rr

Source Code

	FUNCTION poly_rr(x,coeffs)
	REAL(SP), INTENT(IN) :: x
	REAL(SP), DIMENSION(:), INTENT(IN) :: coeffs
	REAL(SP) :: poly_rr
	REAL(SP) :: pow
	REAL(SP), DIMENSION(:), ALLOCATABLE :: vec
	INTEGER(I4B) :: i,n,nn
	n=size(coeffs)
	if (n <= 0) then
		poly_rr=0.0_sp
	else if (n < NPAR_POLY) then
		poly_rr=coeffs(n)
		do i=n-1,1,-1
			poly_rr=x*poly_rr+coeffs(i)
		end do
	else
		allocate(vec(n+1))
		pow=x
		vec(1:n)=coeffs
		do
			vec(n+1)=0.0_sp
			nn=ishft(n+1,-1)
			vec(1:nn)=vec(1:n:2)+pow*vec(2:n+1:2)
			if (nn == 1) exit
			pow=pow*pow
			n=nn
		end do
		poly_rr=vec(1)
		deallocate(vec)
	end if
	END FUNCTION poly_rr