sys/vops/lz/assqd.x


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18

# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.

# ASSQ -- Vector sum of squares.

double procedure assqd (a, npix)
double	sum

double	a[ARB]
int	npix
int	i

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