blob: 92cb036de597dd9014382ead7d867e8911a1f7cb (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <syserr.h>
include <plset.h>
include <plio.h>
# PL_ACCESS -- Return a pointer (type short) to the encoded line list data
# for the indicated mask image line. A valid pointer is always returned;
# if the mask line is empty, the pointer will point to "empty line" linelist.
pointer procedure pl_access (pl, v)
pointer pl #I mask descriptor
long v[PL_MAXDIM] #I coordinates of desired line
int pl_reference()
begin
return (Ref (pl, pl_reference(pl,v)))
end
# PL_REFERENCE -- Return a reference (llbuf offset) to the indicated mask
# image line. A valid offset is always returned; if the mask line is empty,
# the offset will be that of the "empty line" linelist.
int procedure pl_reference (pl, v)
pointer pl #I mask descriptor
long v[PL_MAXDIM] #I coordinates of desired line
int index, i
int totlen, axlen
define oob_ 91
begin
# Compute the index of the line in the line pointer array.
if (PL_NAXES(pl) == 2) {
# Optimized for case naxes=2.
index = v[2]
if (index < 1 || index > PL_AXLEN(pl,2))
goto oob_
} else {
# General case.
index = 1
totlen = 1
do i = 2, PL_NAXES(pl) {
axlen = PL_AXLEN(pl,i)
if (v[i] < 1 || v[i] > axlen)
goto oob_
index = index + totlen * (v[i] - 1)
totlen = totlen * axlen
}
}
return (PL_LP(pl,index))
oob_
call syserr (SYS_PLREFOOB)
end
|
