33     allquads_t* aq = calloc(1, sizeof(allquads_t));  in allquads_init()  local
34     aq->dimquads = 4;  in allquads_init()
35     return aq;  in allquads_init()
38 void allquads_free(allquads_t* aq) {  in allquads_free()  argument
39     free(aq);  in allquads_free()
43     allquads_t* aq = token;  in add_quad()  local
51     quad_write_const(aq->codes, aq->quads, quad, aq->starkd, qb->dimquads, aq->dimcodes);  in add_quad()
55     allquads_t* aq = token;  in check_AB()  local
56     debug("check_AB: iA=%i, iB=%i, aq->starA=%i\n", pq->iA, pq->iB, aq->starA);  in check_AB()
57     return (pq->iA == aq->starA);  in check_AB()
80 int allquads_create_quads(allquads_t* aq) {  in allquads_create_quads()  argument
87     qb->quadd2_high = aq->quad_d2_upper;  in allquads_create_quads()
88     qb->quadd2_low  = aq->quad_d2_lower;  in allquads_create_quads()
89     qb->check_scale_high = aq->use_d2_upper;  in allquads_create_quads()
90     qb->check_scale_low  = aq->use_d2_lower;  in allquads_create_quads()
92     qb->dimquads = aq->dimquads;  in allquads_create_quads()
94     qb->add_quad_token = aq;  in allquads_create_quads()
96     N = startree_N(aq->starkd);  in allquads_create_quads()
106         kdtree_copy_data_double(aq->starkd->tree, 0, N, xyz);  in allquads_create_quads()
126         nq = aq->quads->numquads;  in allquads_create_quads()
139             startree_get(aq->starkd, i, xyzA);  in allquads_create_quads()
141             startree_search_for(aq->starkd, xyzA, aq->quad_d2_upper,  in allquads_create_quads()
150             aq->starA = i;  in allquads_create_quads()
155             qb->check_AB_stars_token = aq;  in allquads_create_quads()
161 …i: wrote %i quads for this star, total %i so far.\n", i+1, N, aq->quads->numquads - nq, aq->quads-…  in allquads_create_quads()
175 int allquads_close(allquads_t* aq) {  in allquads_close()  argument
176     startree_close(aq->starkd);  in allquads_close()
179     if (quadfile_fix_header(aq->quads) ||  in allquads_close()
180         quadfile_close(aq->quads)) {  in allquads_close()
184     if (codefile_fix_header(aq->codes) ||  in allquads_close()
185         codefile_close(aq->codes)) {  in allquads_close()
192 int allquads_open_outputs(allquads_t* aq) {  in allquads_open_outputs()  argument
196     printf("Reading star kdtree %s ...\n", aq->skdtfn);  in allquads_open_outputs()
197     aq->starkd = startree_open(aq->skdtfn);  in allquads_open_outputs()
198     if (!aq->starkd) {  in allquads_open_outputs()
199         ERROR("Failed to open star kdtree %s\n", aq->skdtfn);  in allquads_open_outputs()
202     printf("Star tree contains %i objects.\n", startree_N(aq->starkd));  in allquads_open_outputs()
204     printf("Will write to quad file %s and code file %s\n", aq->quadfn, aq->codefn);  in allquads_open_outputs()
205     aq->quads = quadfile_open_for_writing(aq->quadfn);  in allquads_open_outputs()
206     if (!aq->quads) {  in allquads_open_outputs()
207         ERROR("Couldn't open file %s to write quads.\n", aq->quadfn);  in allquads_open_outputs()
210     aq->codes = codefile_open_for_writing(aq->codefn);  in allquads_open_outputs()
211     if (!aq->codes) {  in allquads_open_outputs()
212         ERROR("Couldn't open file %s to write codes.\n", aq->quadfn);  in allquads_open_outputs()
216     aq->quads->dimquads = aq->dimquads;  in allquads_open_outputs()
217     aq->codes->dimcodes = aq->dimcodes;  in allquads_open_outputs()
219     if (aq->id) {  in allquads_open_outputs()
220         aq->quads->indexid = aq->id;  in allquads_open_outputs()
221         aq->codes->indexid = aq->id;  in allquads_open_outputs()
225     hp = qfits_header_getint(startree_header(aq->starkd), "HEALPIX", -1);  in allquads_open_outputs()
229     aq->quads->healpix = hp;  in allquads_open_outputs()
230     aq->codes->healpix = hp;  in allquads_open_outputs()
232     hpnside = qfits_header_getint(startree_header(aq->starkd), "HPNSIDE", 1);  in allquads_open_outputs()
233     aq->quads->hpnside = hpnside;  in allquads_open_outputs()
234     aq->codes->hpnside = hpnside;  in allquads_open_outputs()
236     hdr = quadfile_get_header(aq->quads);  in allquads_open_outputs()
242     hdr = codefile_get_header(aq->codes);  in allquads_open_outputs()
248     if (quadfile_write_header(aq->quads)) {  in allquads_open_outputs()
249         ERROR("Couldn't write headers to quads file %s\n", aq->quadfn);  in allquads_open_outputs()
252     if (codefile_write_header(aq->codes)) {  in allquads_open_outputs()
253         ERROR("Couldn't write headers to code file %s\n", aq->codefn);  in allquads_open_outputs()
257     if (!aq->use_d2_lower)  in allquads_open_outputs()
258         aq->quad_d2_lower = 0.0;  in allquads_open_outputs()
259     if (!aq->use_d2_upper)  in allquads_open_outputs()
260         aq->quad_d2_upper = 10.0;  in allquads_open_outputs()
262     aq->codes->numstars = startree_N(aq->starkd);  in allquads_open_outputs()
263     aq->codes->index_scale_upper = distsq2rad(aq->quad_d2_upper);  in allquads_open_outputs()
264     aq->codes->index_scale_lower = distsq2rad(aq->quad_d2_lower);  in allquads_open_outputs()
266     aq->quads->numstars = aq->codes->numstars;  in allquads_open_outputs()
267     aq->quads->index_scale_upper = aq->codes->index_scale_upper;  in allquads_open_outputs()
268     aq->quads->index_scale_lower = aq->codes->index_scale_lower;  in allquads_open_outputs()