cumsum_i Function

  1. special_fun_modules.f90
  2. nrutil
  3. cumsum_i

public recursive function cumsum_i(arr, seed) result(ans)

Arguments

Type IntentOptional AttributesName
integer(kind=I4B), intent(in), DIMENSION(:):: arr
integer(kind=I4B), intent(in), optional :: seed

Return Value integer(kind=I4B), DIMENSION(size(arr))


Called by

proc~~cumsum_i~~CalledByGraph proc~cumsum_i cumsum_i proc~cumsum_i->proc~cumsum_i interface~cumsum cumsum interface~cumsum->proc~cumsum_i

Contents

Source Code

cumsum_i

Source Code

	RECURSIVE FUNCTION cumsum_i(arr,seed) RESULT(ans)
	INTEGER(I4B), DIMENSION(:), INTENT(IN) :: arr
	INTEGER(I4B), OPTIONAL, INTENT(IN) :: seed
	INTEGER(I4B), DIMENSION(size(arr)) :: ans
	INTEGER(I4B) :: n,j,sd
	n=size(arr)
	if (n == 0_i4b) RETURN
	sd=0_i4b
	if (present(seed)) sd=seed
	ans(1)=arr(1)+sd
	if (n < NPAR_CUMSUM) then
		do j=2,n
			ans(j)=ans(j-1)+arr(j)
		end do
	else
		ans(2:n:2)=cumsum_i(arr(2:n:2)+arr(1:n-1:2),sd)
		ans(3:n:2)=ans(2:n-1:2)+arr(3:n:2)
	end if
	END FUNCTION cumsum_i