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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
|
include "../../lib/ptkeysdef.h"
include "pexamine.h"
# PT_WTFILE -- Write out the catalogs of selected and rejected objects.
procedure pt_wtfile (apd, key, apout, aprej, deleted, nstars)
int apd # input catalog file descriptor
pointer key # key structure for textfiles
int apout # output catalog file descriptor
int aprej # rejections catalog file descriptor
int deleted[ARB] # deletions array
int nstars # number of stars in the catalog
int i, nselect, ndelete
int pt_selrej()
begin
# ST tables format.
if (key == NULL) {
# Write out the good data.
if (apout != NULL) {
call tbtcre (apout)
call tbhcal (apd, apout)
nselect = 0
do i = 1, nstars {
if (deleted[i] == PX_DELETE)
next
nselect = nselect + 1
call tbrcpy (apd, apout, i, nselect)
}
}
# Write out the deletions.
if (aprej != NULL) {
call tbtcre (aprej)
call tbhcal (apd, aprej)
ndelete = 0
do i = 1, nstars {
if (deleted[i] != PX_DELETE)
next
ndelete = ndelete + 1
call tbrcpy (apd, aprej, i, ndelete)
}
}
# Write out a text file.
} else {
if (pt_selrej (apd, apout, aprej, deleted, nstars) < nstars)
;
}
end
define LEN_LONGLINE 10
# PT_SELREJ -- Select and/or reject records based on evaluating a logical
# expression.
int procedure pt_selrej (tp_in, tp_out, tp_rej, deleted, nstars)
int tp_in # the input catalog file descriptor
int tp_out # the output catalog file descriptor
int tp_rej # the rejections catalog file descriptor
int deleted[ARB] # the deletions array
int nstars # maximum number of stars
int record, nchars, buflen, lenrecord
pointer line, lline
int getline()
begin
# Check that output has been requested.
if (tp_out == NULL && tp_rej == NULL)
return (0)
# Rewind the input file.
call seek (tp_in, BOF)
# Initialize the file read.
record = 0
lenrecord = 0
buflen = LEN_LONGLINE * SZ_LINE
call malloc (line, SZ_LINE, TY_CHAR)
call malloc (lline, buflen, TY_CHAR)
# Loop over the text file records.
repeat {
# Read in a line of the text file.
nchars = getline (tp_in, Memc[line])
if (nchars == EOF)
break
# Determine the type of record.
if (Memc[line] == KY_CHAR_POUND || Memc[line] == KY_CHAR_NEWLINE) {
if (tp_out != NULL)
call putline (tp_out, Memc[line])
if (tp_rej != NULL)
call putline (tp_rej, Memc[line])
} else {
# Reallocate the temporary record space.
if (lenrecord > buflen) {
buflen = buflen + SZ_LINE
call realloc (lline, buflen, TY_CHAR)
}
# Store the record.
call amovc (Memc[line], Memc[lline+lenrecord], nchars)
lenrecord = lenrecord + nchars
Memc[lline+lenrecord] = EOS
# Do the record bookkeeping.
if (Memc[line+nchars-2] != KY_CHAR_CONT) {
# Increment the record counter.
record = record + 1
# Write out the expression.
if ((tp_out != NULL) && (deleted[record] != PX_DELETE))
call putline (tp_out, Memc[lline])
if ((tp_rej != NULL) && (deleted[record] == PX_DELETE))
call putline (tp_rej, Memc[lline])
if (record >= nstars)
break
# Reinitialize the record read.
lenrecord = 0
}
}
}
# Cleanup.
call mfree (line, TY_CHAR)
call mfree (lline, TY_CHAR)
return (record)
end
|
