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
60
61
62
63
64
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <plset.h>
include <plio.h>
# PL_SSIZE -- Set the size of a mask, i.e., given an existing open mask
# descriptor, create an empty mask with the given dimensionality and size.
procedure pl_ssize (pl, naxes, axlen, depth)
pointer pl #I mask descriptor
int naxes #I number of axes (dimensionality of mask)
long axlen[ARB] #I length of each axis
int depth #I mask depth, bits
int npix, i
pointer sp, px, lp
int pl_p2ls()
errchk malloc, calloc, mfree
begin
npix = axlen[1]
call smark (sp)
call salloc (px, npix, TY_SHORT)
# Initialize the old descriptor.
if (PL_LPP(pl) != NULL)
call mfree (PL_LPP(pl), TY_INT)
if (PL_LLBP(pl) != NULL)
call mfree (PL_LLBP(pl), TY_SHORT)
call amovki (1, PL_PLANE(pl,1), PL_MAXDIM)
# Set up the empty descriptor.
PL_NAXES(pl) = naxes
if (depth > 0)
PL_MAXVAL(pl) = MV(depth)
else
PL_MAXVAL(pl) = MV(PL_MAXDEPTH)
call amovl (axlen, PL_AXLEN(pl,1), naxes)
do i = naxes + 1, PL_MAXDIM
PL_AXLEN(pl,i) = 1
# Allocate the line list buffer.
PL_MAXLINE(pl) = (axlen[1] * 3) + LL_CURHDRLEN
PL_LLLEN(pl) = PL_LLBUFLEN
call malloc (PL_LLBP(pl), PL_LLBUFLEN, TY_SHORT)
lp = PL_LLBP(pl)
# Encode the empty line line-list.
call aclrs (Mems[px], npix)
PL_LLOP(pl) = pl_p2ls (Mems[px], 1, Mems[lp], npix)
# Set up the initial line list index (all lines empty).
PL_NLP(pl) = 1
do i = 2, naxes
PL_NLP(pl) = PL_NLP(pl) * axlen[i]
call calloc (PL_LPP(pl), PL_NLP(pl), TY_INT)
# Set up the LL header for the empty line.
LP_NREFS(lp) = PL_NLP(pl)
LP_SETBLEN(lp, PL_LLOP(pl))
call sfree (sp)
end
|
