aboutsummaryrefslogtreecommitdiff
path: root/sys/vops/lz/asumd.x
blob: 24e4e7a94f81903c6fede88f8889fae324596397 (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.

# ASUM -- Vector sum.  Take care to prevent integer overflow by returning
# a floating point sum.

double procedure asumd (a, npix)

double	a[ARB]
int	npix
int	i

double	sum

begin
	sum = 0.0D0
	do i = 1, npix
	    sum = sum + a[i]
	
	return (sum)
end