summaryrefslogtreecommitdiffstats
path: root/funtools/search.c
diff options
context:
space:
mode:
authorWilliam Joye <wjoye@cfa.harvard.edu>2016-10-26 21:13:00 (GMT)
committerWilliam Joye <wjoye@cfa.harvard.edu>2016-10-26 21:13:00 (GMT)
commitda2e3d212171bbe64c1af39114fd067308656990 (patch)
tree9601f7ed15fa1394762124630c12a792bc073ec2 /funtools/search.c
parent76b109ad6d97d19ab835596dc70149ef379f3733 (diff)
downloadblt-da2e3d212171bbe64c1af39114fd067308656990.zip
blt-da2e3d212171bbe64c1af39114fd067308656990.tar.gz
blt-da2e3d212171bbe64c1af39114fd067308656990.tar.bz2
rm funtools for update
Diffstat (limited to 'funtools/search.c')
-rw-r--r--funtools/search.c1107
1 files changed, 0 insertions, 1107 deletions
diff --git a/funtools/search.c b/funtools/search.c
deleted file mode 100644
index 68c01b4..0000000
--- a/funtools/search.c
+++ /dev/null
@@ -1,1107 +0,0 @@
-/* search.c
- */
-/* Starbase Data Tables - An Ascii Database for UNIX
-
- copyright 1991, 1993, 1995, 1999, 2003 John B. Roll jr.
- */
-
-extern int N;
-
-#include <math.h>
-#include <ctype.h>
-#include <except.h>
-
-#include <table.h>
-#include <tablesearch.h>
-#include <filelist.h>
-
-double slaDsep();
-void radecbox();
-
-exception ex_search = "search";
-
-typedef void *(*vvector)();
-
-TableSearch table_search();
-TableSearch table_searchsetup();
-char *table_searchcheckindexpattern();
-Range *table_srange2();
-
-typedef struct Axis {
- int type;
- int once;
- int c1;
- int c2;
- char *value;
- double lo;
- double hi;
- double wp;
- Range *R;
- int method;
- TableSearch search;
- struct Axis *pair;
-} Axis;
-
-char *table_searchcheckindexpattern();
-
-
-void table_searchpairs(table, axes, naxes, filecache)
- char *table;
- Axis **axes;
- int naxes;
- void *filecache;
-{
- char *here;
- char directo[1024];
- int i, j;
-
- TableSearch search = table_searchsetup(table);
- int Numeric;
- int Method;
-
- strcpy(directo, table);
- if ( (here = strrchr(directo, '/')) ) {
- *here = '\0';
- here++;
- } else {
- strcpy(directo, ".");
- here = table;
- }
-
- for ( i = 0; i < naxes; i++ ) {
- for ( j = 0; j < naxes; j++ ) {
- char *name;
- char *col1;
- char *col2;
- filelist list;
- int len;
- char pattern[1024];
-
- if ( i == j ) continue;
- if ( axes[i]->search != NULL ) continue;
- if ( axes[j]->search != NULL ) continue;
-
- col1 = table_colnam(search->THead, axes[i]->c2);
- col2 = table_colnam(search->THead, axes[j]->c2);
-
- sprintf(pattern, "%s.%s.%s[-.:][sbhir]", here, col1, col2);
-
- if ( (name = table_searchcheckindexpattern(search, directo, pattern)) ) {
- table_searchsetupindex(search, name, col1);
- table_searchsetupcache(search, NULL);
-
- axes[i]->search = search;
- axes[j]->search = search;
- axes[i]->search->Index = name;
- axes[i]->search->Tcolumn = table_colnum(search->THead, col1);
- axes[i]->search->Indxcol = table_colnum(search->IHead, col1);
- axes[i]->search->Tcolum2 = table_colnum(search->THead, col2);
- axes[i]->search->Indxco2 = table_colnum(search->IHead, col2);
- axes[i]->type = 'S';
- axes[j]->type = '4';
- axes[j]->pair = axes[i];
-
- search = table_searchsetup(table); /* Leak one search descriptor */
- }
- } }
-}
-
-
-typedef struct CsysDef {
- void (*setup)();
- int (*filte)();
- int nrange;
- double *range;
- int naxes;
- Axis *axis;
- char *data;
-} CsysDef;
-
-typedef struct CsysType {
- char option;
- char *name;
- int (*filte)();
- void (*setup)();
- char *(*parse)();
-} CsysType;
-
-typedef struct RCache {
- struct _CNode cache;
-} *RCache;
-
-typedef struct JoinData {
- CsysDef *csys;
- int nsys;
- Range *r;
- File ofile;
- int print;
- int setcnt;
- RCache cache;
- int posit;
- TableHead T1;
- TableRow R1;
-} JoinData;
-
-rowcmp(a, b)
- CNode a;
- CNode b;
-{
- return a->hash == b->hash;
-}
-
-char *cartparse(csys, option)
- CsysDef *csys;
- char *option;
-{
- int i;
-
- if ( *option )
- csys->naxes = atoi(option);
- else csys->naxes = 1;
- csys->nrange = 1;
-
- if ( !csys->naxes )
- except(ex_search, "error in Cartesian coordinate system option");
-
- Malloc(csys->axis, sizeof(Axis) * csys->naxes);
-
- for ( i = 0; i < csys->naxes; i++ ) {
- csys->axis[i].type = '2';
- csys->axis[i].once = 0;
- csys->axis[i].method = TSEARCH_SEXAGAS;
- }
- return "C";
-}
-
-void cartsetup(csys, row)
- CsysDef *csys;
- TableRow row;
-{
- int i;
-
- for ( i = 0; i < csys->naxes; i++ ) {
- csys->axis[i].lo = table_colvald(row, csys->axis[i].c1) - csys->range[0];
- csys->axis[i].hi = table_colvald(row, csys->axis[i].c1) + csys->range[0];
- }
-}
-
-cartfilte(csys, row1, row2)
- CsysDef *csys;
- TableRow row1;
- TableRow row2;
-{
- int i;
- double square = 0.0;
-
- for ( i = 0; i < csys->naxes; i++ ) {
- double d = table_colvald(row1, csys->axis[i].c1)
- - table_colvald(row2, csys->axis[i].c2);
-
- square += d * d;
- }
-
- if ( sqrt(square) <= csys->range[0] ) return 1;
- else return 0;
-}
-
-char *sphrparse(csys, option)
- CsysDef *csys;
- char *option;
-{
- int i;
-
- if ( *option == '1' || *option == '2' ) {
- csys->naxes = *option++ - '0';
- } else
- csys->naxes = 2;
-
- if ( csys->naxes == 1 )
- csys->filte = NULL;
-
- csys->nrange = 1;
-
- Malloc(csys->axis, sizeof(Axis) * csys->naxes);
-
- for ( i = 0; i < csys->naxes; i++ ) {
- csys->axis[i].type = '2';
- csys->axis[i].once = 0;
- csys->axis[i].method = TSEARCH_SEXAGAS;
- }
-
- /* Parse units options. Set up wrap values.
- */
- NewString(csys->data, "hdd");
- for ( i = 0; i < 3 && option[i]; i++ )
- csys->data[i] = option[i];
-
- for ( i = 0; i < csys->naxes; i++ ) {
- switch ( csys->data[i] ) {
- case 'h': csys->axis[i].wp = 24.0; break;
- case 'd': csys->axis[i].wp = 360.0; break;
- case 'r': csys->axis[i].wp = X_2PI; break;
- default:
- except(ex_search, "bad units option to spherical system\n");
- }
- }
-
- return "C";
-}
-
-void sphrsetup(csys, row)
- CsysDef *csys;
- TableRow row;
-{
- double ra, dec;
- double width;
- double r0, r1, d0, d1;
-
- if ( csys->naxes == 1 ) {
- csys->axis[0].lo = table_colvald(row, csys->axis[0].c1) -
- csys->range[0];
- if ( csys->axis[0].lo < 0 ) csys->axis[0].lo += csys->axis[0].wp;
-
- csys->axis[0].hi = table_colvald(row, csys->axis[0].c1) +
- csys->range[0];
- if ( csys->axis[0].hi > csys->axis[0].wp )
- csys->axis[0].hi -= csys->axis[0].wp;
-
- return;
- }
-
-
- /* Input units conversion.
- */
- switch ( csys->data[0] ) {
- case 'h': ra = table_colvald(row, csys->axis[0].c1); break;
- case 'd': ra = d2h(table_colvald(row, csys->axis[0].c1)); break;
- case 'r': ra = r2h(table_colvald(row, csys->axis[0].c1)); break;
- }
- switch ( csys->data[1] ) {
- case 'h': dec = h2d(table_colvald(row, csys->axis[1].c1)); break;
- case 'd': dec = table_colvald(row, csys->axis[1].c1); break;
- case 'r': dec = r2d(table_colvald(row, csys->axis[1].c1)); break;
- }
- switch ( csys->data[2] ) {
- case 'h': width = h2d(csys->range[0]); break;
- case 'd': width = csys->range[0]; break;
- case 'r': width = r2d(csys->range[0]); break;
- }
-
- radecbox(ra, dec, width * 2, &r0, &r1, &d0, &d1);
-
- /* Output units conversion
- */
- switch ( csys->data[0] ) {
- case 'h': csys->axis[0].lo = r0;
- csys->axis[0].hi = r1;
- break;
- case 'd': csys->axis[0].lo = h2d(r0);
- csys->axis[0].hi = h2d(r1);
- break;
- case 'r': csys->axis[0].lo = h2r(r0);
- csys->axis[0].hi = h2r(r1);
- break;
- }
- switch ( csys->data[1] ) {
- case 'h': csys->axis[1].lo = d2h(d0);
- csys->axis[1].hi = d2h(d1);
- break;
- case 'd': csys->axis[1].lo = d0;
- csys->axis[1].hi = d1;
- break;
- case 'r': csys->axis[1].lo = d2r(d0);
- csys->axis[1].hi = d2r(d1);
- break;
- }
-}
-
-int sphrfilte(csys, row1, row2)
- CsysDef *csys;
- TableRow row1;
- TableRow row2;
-{
- double r1, r2, d1, d2;
- double radius;
-
- switch ( csys->data[0] ) {
- case 'h': r1 = h2r(table_colvald(row1, csys->axis[0].c1));
- r2 = h2r(table_colvald(row2, csys->axis[0].c2));
- break;
- case 'd': r1 = d2r(table_colvald(row1, csys->axis[0].c1));
- r2 = d2r(table_colvald(row2, csys->axis[0].c2));
- break;
- case 'r': r1 = table_colvald(row1, csys->axis[0].c1);
- r2 = table_colvald(row2, csys->axis[0].c2);
- break;
- }
- switch ( csys->data[1] ) {
- case 'h': d1 = h2r(table_colvald(row1, csys->axis[1].c1));
- d2 = h2r(table_colvald(row2, csys->axis[1].c2));
- break;
- case 'd': d1 = d2r(table_colvald(row1, csys->axis[1].c1));
- d2 = d2r(table_colvald(row2, csys->axis[1].c2));
- break;
- case 'r': d1 = table_colvald(row1, csys->axis[1].c1);
- d2 = table_colvald(row2, csys->axis[1].c2);
- break;
- }
- switch ( csys->data[2] ) {
- case 'h': radius = h2r(csys->range[0]); break;
- case 'd': radius = d2r(csys->range[0]); break;
- case 'r': radius = csys->range[0]; break;
- }
-
- if ( slaDsep(r1, d1, r2, d2) <= radius ) return 1;
- else return 0;
-}
-
-char *singparse(csys, option)
- CsysDef *csys;
- char *option;
-{
-
- csys->naxes = 1;
- csys->nrange = 0;
-
- Malloc(csys->axis, sizeof(Axis) * csys->naxes);
-
- csys->axis[0].type = '1';
- csys->axis[0].once = 0;
- csys->axis[0].method = TSEARCH_ASCII;
- return "C";
-}
-
-void singsetup(csys, row)
- CsysDef *csys;
- TableRow row;
-{
- csys->axis[0].value = table_colval(row, csys->axis[0].c1);
-}
-
-
-char *rangparse(csys, option)
- CsysDef *csys;
- char *option;
-{
-
- csys->naxes = 1;
- csys->nrange = 0;
-
- Malloc(csys->axis, sizeof(Axis) * 3);
-
- if ( *option == ':' ) {
- double wrap = SAOstrtod(&option[1], NULL);
-
- csys->axis[0].wp = wrap;
- }
-
- csys->axis[0].type = '2';
- csys->axis[0].once = 0;
- csys->axis[0].method = TSEARCH_SEXAGAS;
- return "211";
-}
-
-void rangsetup(csys, row)
- CsysDef *csys;
- TableRow row;
-{
- csys->naxes = 1;
-
- csys->axis[0].lo = table_colvald(row, csys->axis[1].c1);
- csys->axis[0].hi = table_colvald(row, csys->axis[2].c1);
-}
-
-char *valuparse(csys, option)
- CsysDef *csys;
- char *option;
-{
- csys->naxes = 1;
- csys->nrange = 0;
-
- Malloc(csys->axis, sizeof(Axis) * (csys->naxes + 0));
- csys->axis[0].type = '1';
- csys->axis[0].once = 1;
- csys->axis[0].method = TSEARCH_ASCII;
-
- return "2V";
-}
-
-void valusetup(csys, row)
- CsysDef *csys;
- TableRow row;
-{
- TableRow X = table_rowloc(NULL, 0, 0);
-
- csys->axis[0].R = range_union(table_svalue(csys->axis[0].search
- , csys->axis[0].value
- , (vvector) range_add, NULL, X));
-
- table_rowfree(X);
- csys->setup = NULL;
-}
-
-char *limtparse(csys, option)
- CsysDef *csys;
- char *option;
-{
- csys->naxes = 1;
- csys->nrange = 2;
-
- Malloc(csys->axis, sizeof(Axis) * csys->naxes);
- csys->axis[0].type = '2';
- csys->axis[0].once = 1;
- csys->axis[0].method = TSEARCH_SEXAGAS;
-
- return "2";
-}
-
-
-void limtsetup(csys, row)
- CsysDef *csys;
- TableRow row;
-{
- csys->axis[0].lo = csys->range[0];
- csys->axis[0].hi = csys->range[1];
-
- csys->setup = NULL;
- csys->naxes = 0;
-}
-
-
-CsysType CsysTypes[] = {
- { 'C' , "cartesian", cartfilte, cartsetup, cartparse }
- , { 'S', "spherical", sphrfilte, sphrsetup, sphrparse }
- , { 'V', "col value", NULL, NULL , valuparse }
- , { 'Z', "col range", NULL, rangsetup, rangparse }
- , { 'R', "col limit", NULL, limtsetup, limtparse }
- , { '\0', NULL, NULL, NULL, NULL }
- , { '\0', "singleaxs", NULL, singsetup, singparse }
-#define N_CsysTypes 6
-};
-
-
-void *joinem(data, start, end, row, table, file)
- void *data;
- off_t start;
- off_t end;
- TableRow row;
- TableHead table;
- File file;
-{
- JoinData *d = (JoinData *) data;
- int i;
-
- struct _CNode node;
-
- node.hash = start;
-
- /* If the row won't pass the filter - bail.
- */
- for ( i = 0; i < d->nsys; i++ ) {
- if ( ! d->csys[i].filte) { continue; }
- if ( !(*d->csys[i].filte)(&d->csys[i], d->R1, row) ) {
- return data;
- }
- }
-
- if ( !d->print ) {
- d->r = range_add(d->r, start, end);
- return data;
- }
-
- if ( d->posit ) {
- fprintf(d->ofile, OffFormat " %d\n", start, end-start);
- return data;
- }
-
- if ( d->setcnt ) {
- if ( cache_hashe(d->cache, &node, NULL) != NULL )
- return data;
-
- cache_cache(d->cache, &node);
- fprintf(d->ofile, "%d\t", d->setcnt);
- }
-
- if ( d->T1 )
- table_rowput(d->ofile, d->T1, d->R1, NULL, NULL, 0);
-
- table_rowput(d->ofile, table, row, NULL, NULL, 0);
- return data;
-}
-
-char *trim();
-
-main(argc, argv)
- int argc;
- char *argv[];
-{
- int n, i, j, k;
- int nsys = 0;
- int asys = 0;
-
- Axis **axes;
- int naxes = 0;
-
- int *lcol = NULL;
- int ncol = 0;
-
- CsysDef *csys = NULL;
-
- File ifile = Stdin;
- File ofile = Stdout;
- void *obuff;
-
- File F1, F2 = NULL;
- TableHead T1 = NULL, T2;
- char *file1 = "-";
- char *file2 = NULL;
-
- TableRow row, row1 = NULL, X;
-
- JoinData jdata;
-
- char *zeval = "0.0";
- char *args = NULL;
-
- int header = 1;
- int invert = 0;
- int joiner = 0;
- int count = 0;
- int setcnt = 0;
- int posit = 0;
-
- off_t Start, End;
-
- if ( argc < 2 ) {
- FPrint(Stderr
- , "usage: search table [-p] [-s] [-j] [-i keytable] [-m[bhisr]]\n");
- FPrint(Stderr, " [Column]\n");
- FPrint(Stderr, " [-R Column <Min> <Max>]\n");
- FPrint(Stderr, " [-S[uuu] Column Column Range]\n");
- FPrint(Stderr, " [-C[n] Column [... Column<n>] Range]\n");
- FPrint(Stderr, " [-Z[wrap] Column Column]\n");
- exit(1);
- }
-
- /* Look ahead for file option.
- */
- for ( i = 1; i < argc - 1 ; i++ ) {
- if ( argv[i][0] == '-' && argv[i][1] == 'i' ) {
- if ( (ifile = (!strcmp("-", argv[i+1]))
- ? Stdin
- : Open(argv[i+1], "r")) == NULL ) {
- except(ex_search
- , "search: can't open input file: %s : %s"
- , argv[i+1], exc_oserror());
- } else {
- for ( ; i < argc - 2; i++ )
- argv[i] = argv[i+2]; /* Ooch these over */
- argc -= 2;
- break;
- }
- }
- }
-
- F1 = ifile;
-
- for ( n = 1; n < argc; n++ ) {
- char *arg;
-
- if ( argv[n][0] == '-' ) {
- switch ( argv[n][1] ) {
- case 'c' : count = 1; continue;
- case 'h' : header = 0; continue;
- case 'I' : invert = 1; continue;
- case 'p' : posit = 1; continue;
- case 'j' : joiner = 1; continue;
- case 's' : setcnt = 1; continue;
- case 'o' :
- if ( (ofile = (!strcmp("-", argv[n+1]))
- ? Stdout : Open(argv[n+1], "w")) == NULL ) {
- except(ex_search
- , "can't open output file: %s : %s"
- , argv[n+1], exc_oserror());
- }
- n++;
- continue;
- default:
-
- for ( i = 0; CsysTypes[i].name != NULL; i++ ) {
- if ( CsysTypes[i].option == argv[n][1] )
- break;
- }
-
- if ( CsysTypes[i].name == NULL )
- except(ex_search
- , "search: unknown option: %c", argv[n][1]);
-
- arg = &argv[n][2];
- n++;
- break;
- }
- } else {
- if ( F2 ) {
- i = N_CsysTypes;
- arg = NULL;
- } else {
- if ( (F2 = Open(argv[n], "r")) == NULL )
- except(ex_search, "can't open table: %s : %s", argv[n]
- , exc_oserror());
- if ( (T2 = table_header(F2, 0)) == NULL )
- except(ex_search, "can't read table header from: %s\n"
- , argv[n]);
- file2 = argv[n];
-
- if ( (Start = Seek(F2, 0, SEEK_CUR)) == -1 ) {
- perror("search:");
- exit(1);
- }
- End = Seek(F2, 0 , SEEK_END);
- Seek(F2, Start, SEEK_SET);
-
- continue;
- }
- }
-
- if ( file2 == NULL ) except(ex_search, "no table to search");
-
- if ( nsys >= asys ) {
- asys *= 2;
- asys += 4;
- ReAlloc(csys, asys * sizeof(CsysDef));
- }
-
- csys[nsys].setup = CsysTypes[i].setup;
- csys[nsys].filte = CsysTypes[i].filte;
-
- args = CsysTypes[i].parse(&csys[nsys], arg);
-
- if ( csys[nsys].nrange )
- Malloc(csys[nsys].range, sizeof(double) * csys[nsys].nrange);
-
- for ( i = 0; i < csys[nsys].naxes || (*args && *args != 'C')
- ; n++, i++ ) {
- char *c1;
- char *c2;
- char *nchr;
-
- if ( n >= argc )
- except(ex_search
- , "not enough columns for this coordinate system");
-
- c1 = argv[n];
-
- if ( *args != 'V' && (nchr = strchr(c1, '.')) )
- *nchr++ = '\0';
- else nchr = NULL;
-
- if ( *args == 'C' && (c2 = strchr(c1, ':')) )
- *c2++ = '\0';
- else c2 = c1;
-
- if ( *args == 'C' || *args == '1' ) {
- if ( T1 == NULL && (T1 = table_header(F1, TABLE_PARSE)) == NULL )
- except(ex_search
- , "can't read table header from: %s\n", file1);
-
- if ( !(csys[nsys].axis[i].c1 = table_colnum(T1, c1)) )
- except(ex_search, "no column %s in keytable", c1);
- }
-
- if ( *args == 'C' || *args == '2' ) {
- if ( !(csys[nsys].axis[i].c2 = table_colnum(T2, c2)) )
- except(ex_search, "no column %s in table %s\n", c2, file2);
-
- }
-
- if ( *args == 'V' ) {
- csys[nsys].axis[i-1].value = argv[n];
- csys[nsys].axis[i-1].hi = SAOstrtod(argv[n], NULL);
- csys[nsys].axis[i-1].lo = SAOstrtod(argv[n], NULL);
- } else {
- csys[nsys].axis[i].value = NULL;
- csys[nsys].axis[i].R = NULL;
- }
-
- if ( *args != 'C' ) args++;
- }
-
- for ( i = 0; i < csys[nsys].nrange; i++, n++ ) {
- char *here = NULL;
-
- if ( n >= argc )
- except(ex_search
- , "not enough ranges given for this coordinate system.\n");
-
- csys[nsys].range[i] = SAOstrtod(argv[n], &here);
- if ( *here )
- except(ex_search, "this doesn't look like a range: %s"
- , argv[n]);
- }
-
- naxes += csys[nsys].naxes;
- nsys++;
- n--;
- }
-
- Malloc(axes, naxes * sizeof(Axis *));
-
- for ( i = k = 0; i < nsys; i++ ) {
- for ( j = 0; j < csys[i].naxes; j++, k++ ) {
- axes[k] = &csys[i].axis[j];
- axes[k]->search = NULL;
- }
- }
-
- table_searchpairs(file2, axes, naxes, NULL);
-
- for ( k = 0; k < naxes; k++ ) {
- if ( axes[k]->search == NULL ) {
- if ( (axes[k]->search
- = table_search(file2, table_colnam(T2, axes[k]->c2)
- , axes[k]->method, NULL)) == NULL ) {
-
- except(ex_search
- , "can't access index for: %s column: %s\n", file2
- , table_colnam(T2, axes[k]->c2));
- }
- }
- }
-
- row = table_rowloc(NULL, 0, 0);
- X = table_rowloc(NULL, 0, 0);
-
-
- if ( posit ) {
- fprintf(ofile, "Offset Length\n");
- fprintf(ofile, "------ ------\n");
- } else
- if ( T1 && joiner ) {
- if ( header )
- table_hdrput(ofile, T1, TBLHDR_VALUES | TBLHDR_TEXT
- , NULL, NULL, 0);
- table_ors(T1, '\t');
-
- if ( header ) {
- /* rename cols and generate the headline
- */
- Malloc(lcol, (T1->header->ncol + T2->header->ncol)
- * sizeof(int));
-
- if ( setcnt ) fprintf(ofile, "SetNumber\t");
-
- for ( j = 1; j <= T1->header->ncol; j++ )
- /* look up column name in the other table.
- */
- lcol[ncol++] = hdrfix(ofile, table_colval(T1->header, j)
- , table_colnum(T2
- , trim(table_colval(T1->header, j))) ? 1 : 0
- , '\t');
-
- for ( j = 1; j <= T2->header->ncol; j++ )
- lcol[ncol++] = hdrfix(ofile, table_colval(T2->header, j)
- , table_colnum(T1
- , trim(table_colval(T2->header, j))) ? 2 : 0
- , j == T2->header->ncol ? '\n' : '\t');
-
- if ( setcnt ) fprintf(ofile, "---------\t");
-
- /* Output dashline
- */
- for ( i = 0; i < ncol; i++ ) {
- table_dashes(ofile , lcol[i] , (i == ncol - 1) ? '\n'
- : '\t');
- }
- }
- } else {
- table_hdrput(ofile, T2, TBLHDR_VALUES | TBLHDR_TEXT, NULL, NULL, 0);
- if ( setcnt ) fprintf(ofile, "SetNumber\t");
- table_hdrput(ofile, T2, TBLHDR_HEADLINE, NULL, NULL, 0);
- if ( setcnt ) fprintf(ofile, "---------\t");
- table_hdrput(ofile, T2, TBLHDR_DASHLINE, NULL, NULL, 0);
- }
-
-
- {
- Range *U = NULL; /* Union of all keytable rows */
-
- jdata.ofile = ofile;
-
- jdata.csys = csys;
- jdata.nsys = nsys;
-
- jdata.T1 = joiner ? T1 : NULL;
- jdata.r = NULL;
- jdata.print = joiner;
- jdata.setcnt= setcnt;
- jdata.posit = posit;
- if ( setcnt && jdata.print == 0 ) jdata.print = 1;
-
- if ( setcnt ) {
- jdata.cache = (RCache) cache_init(sizeof(struct RCache), 7967, 0
- , NULL, NULL, rowcmp);
- } else {
- jdata.cache = NULL;
- }
-
- /* For each row in the input.
- */
- while ( T1 == NULL
- || (row1 = table_rowget(F1, T1, row, NULL, NULL, 0)) ) {
- Range *R = NULL; /* Intersection of all axes */
- Range *r; /* Search of one axis */
-
- if ( count != 0 ) { progress(count); count++; }
-
- jdata.R1 = row1;
-
-
- /* Call the setup function on each csys
- */
- for ( i = 0; i < nsys; i++ ) {
- if ( csys[i].setup ) csys[i].setup(&csys[i], row1);
- }
-
- /* Build the range set
- */
- for ( i = 0; i < naxes; i++ ) {
- switch ( axes[i]->type ) {
- case 'R': {
- r = range_copy(axes[i]->R);
- break;
- }
- case '1': {
- r = table_svalue(axes[i]->search
- , axes[i]->value
- , (vvector) range_add, NULL, X);
- break;
- }
- case '2': {
- double hi = axes[i]->hi;
- double lo = axes[i]->lo;
- double wp = axes[i]->wp ? axes[i]->wp
- : 0.0 ;
-
- char hival[64];
- char loval[64];
- char wpval[64];
-
- sprintf(hival, "%.21g", hi);
- sprintf(loval, "%.21g", lo);
- sprintf(wpval, "%.21g", wp);
-
- if ( lo > hi ) {
- r = table_srange(axes[i]->search
- , loval, wpval, (vvector) range_add, NULL, X);
- r = table_srange(axes[i]->search
- , zeval, hival, (vvector) range_add, r, X);
- } else {
- r = table_srange(axes[i]->search
- , loval, hival, (vvector) range_add, NULL, X);
- }
- break;
- }
- case '4': {
- double hi1 = axes[i]->pair->hi;
- double lo1 = axes[i]->pair->lo;
- double wp1 = axes[i]->pair->wp ? axes[i]->pair->wp : 0.0 ;
-
- double hi2 = axes[i]->hi;
- double lo2 = axes[i]->lo;
- double wp2 = axes[i]->wp ? axes[i]->wp : 0.0 ;
-
- char hival1[64];
- char loval1[64];
- char wpval1[64];
-
- char hival2[64];
- char loval2[64];
- char wpval2[64];
-
- sprintf(hival1, "%.21g", hi1);
- sprintf(loval1, "%.21g", lo1);
- sprintf(wpval1, "%.21g", wp1);
-
- sprintf(hival2, "%.21g", hi2);
- sprintf(loval2, "%.21g", lo2);
- sprintf(wpval2, "%.21g", wp2);
-
- if ( lo1 > hi1 ) {
- if ( lo2 > hi2 ) {
- r = table_srange2(axes[i]->pair->search
- , loval1, wpval1
- , loval2, wpval2
- , (vvector) range_add, NULL, X);
- r = table_srange2(axes[i]->pair->search
- , loval1, wpval1
- , wpval2, hival2
- , (vvector) range_add, NULL, X);
- r = table_srange2(axes[i]->pair->search
- , wpval1, hival1
- , loval2, wpval2
- , (vvector) range_add, r, X);
- r = table_srange2(axes[i]->pair->search
- , wpval1, hival1
- , wpval2, hival2
- , (vvector) range_add, r, X);
- } else {
- r = table_srange2(axes[i]->pair->search
- , wpval1, hival1
- , loval2, hival2
- , (vvector) range_add, r, X);
- r = table_srange2(axes[i]->pair->search
- , wpval1, hival1
- , loval2, hival2
- , (vvector) range_add, r, X);
- }
- } else {
- if ( lo2 > hi2 ) {
- r = table_srange2(axes[i]->pair->search
- , loval1, hival1
- , loval2, wpval2
- , (vvector) range_add, r, X);
- r = table_srange2(axes[i]->pair->search
- , loval1, hival1
- , wpval2, hival2
- , (vvector) range_add, r, X);
- } else {
- r = table_srange2(axes[i]->pair->search
- , loval1, hival1
- , loval2, hival2
- , (vvector) range_add, NULL, X);
- }
- }
- break;
- }
- case 'S': { continue; }
- default: {
- fprintf(stderr, "search error unknown axis type %c\n",
- axes[i]->type);
- exit(1);
- }
- }
- r = range_union(r);
-
- if ( axes[i]->once ) {
- range_copy(r);
- axes[i]->type = 'R';
- }
-
- if ( R != NULL ) R = range_inter(R, r);
- else R = r;
-
- if ( R == NULL ) break;
- }
-
- /* Filter/Output the Range.
- */
- if ( R != NULL ) {
- table_rangerows(F2, T2, X, R, joinem, &jdata);
- Free(R);
-
- if ( !joiner ) {
- U = range_union(range_cat(U, jdata.r));
- jdata.r = NULL;
- }
- if ( jdata.setcnt ) jdata.setcnt++;
- }
-
- if ( T1 == NULL ) break;
- }
-
- if ( invert ) {
- U = range_invrt(U, Start, End);
- }
-
- if ( U ) {
- jdata.print = 1;
- jdata.nsys = 0;
- table_rangerows(F2, T2, X, U, joinem, &jdata);
- }
- }
- exit(0);
-}
-
-
-void radecbox(ra, dec, width, r1, r2, d1, d2)
- double ra;
- double dec;
- double width;
- double *r1, *r2, *d1, *d2;
-{
- double cosdec;
-
- *d1 = dec - width / 2.0;
- if ( *d1 <= -90.0 ) {
- *d1 = -90.0;
- *d2 = dec + width / 2.0;
- *r1 = 0.0;
- *r2 = 24.0;
- } else {
- *d2 = dec + width / 2.0;
- if ( *d2 >= 90.0 ) {
- *d1 = dec - width / 2.0;
- *d2 = 90.0;
- *r1 = 0.0;
- *r2 = 24.0;
- } else {
- if ( dec > 0.0 ) cosdec = Abs(cos(d2r(*d1)));
- else cosdec = Abs(cos(d2r(*d2)));
-
- *r1 = ra - d2h(width) / 2 / cosdec;
- *r2 = ra + d2h(width) / 2 / cosdec;
-
- if ( *r1 < 0.0 ) *r1 += 24;
- if ( *r2 > 24.0 ) *r2 -= 24;
- }
- }
-}
-
-
-void slaDcs2c ( double a, double b, double v[3] )
-{
- double cosb;
-
- cosb = cos ( b );
- v[0] = cos ( a ) * cosb;
- v[1] = sin ( a ) * cosb;
- v[2] = sin ( b );
-}
-
-double slaDsep ( double a1, double b1, double a2, double b2 )
-{
- int i;
- double d, v1[3], v2[3], s2, c2;
-
- slaDcs2c ( a1, b1, v1 );
- slaDcs2c ( a2, b2, v2 );
-
- s2 = 0.0;
- for ( i = 0; i < 3; i++ ) {
- d = v1[i] - v2[i];
- s2 += d * d;
- }
- s2 /= 4.0;
-
- c2 = 1.0 - s2;
- return 2.0 * atan2 ( sqrt ( s2 ), sqrt ( Max ( 0.0, c2 )));
-}
-
-
-progress(count)
- int count;
-{
- if ( (count % 1000) == 0 ) {
- if ( (count % 100000) == 0 ) {
- FPrint(Stderr, "=");
- } else {
- if ( (count % 10000) == 0 ) {
- FPrint(Stderr, "+");
- } else {
- FPrint(Stderr, "*");
- }
- }
-
- Flush(Stderr);
- }
-}