diff options
Diffstat (limited to 'sys/vops/ak/absul.x')
| -rw-r--r-- | sys/vops/ak/absul.x | 35 |
1 files changed, 35 insertions, 0 deletions
diff --git a/sys/vops/ak/absul.x b/sys/vops/ak/absul.x new file mode 100644 index 00000000..ff803cc6 --- /dev/null +++ b/sys/vops/ak/absul.x @@ -0,0 +1,35 @@ +# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc. + +# ABSU -- Vector block sum. Each pixel in the output vector is the +# sum of the input vector over a block of pixels. The input vector must +# be at least (nblocks * npix_per_block) pixels in length. + +procedure absul (a, b, nblocks, npix_per_block) + +long a[ARB] # input vector +long b[nblocks] # output vector +int nblocks # number of blocks (pixels in output vector) +int npix_per_block # number of input pixels per block + +real sum + +int i, j +int block_offset, next_block, block_length + +begin + block_offset = 1 + block_length = npix_per_block + + if (block_length <= 1) + call amovl (a[block_offset], b, nblocks) + else { + do j = 1, nblocks { + next_block = block_offset + block_length + sum = 0 + do i = block_offset, next_block - 1 + sum = sum + a[i] + b[j] = sum + block_offset = next_block + } + } +end |