123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968969970971972973974975976977978979980981982983984985986987988989990991992993994995996997998999100010011002100310041005100610071008100910101011101210131014101510161017101810191020102110221023102410251026102710281029103010311032103310341035103610371038103910401041104210431044104510461047104810491050105110521053105410551056105710581059106010611062106310641065106610671068106910701071107210731074107510761077107810791080108110821083108410851086108710881089109010911092109310941095109610971098109911001101110211031104110511061107110811091110111111121113111411151116111711181119112011211122112311241125112611271128112911301131113211331134113511361137113811391140114111421143114411451146114711481149115011511152115311541155115611571158115911601161116211631164116511661167116811691170117111721173117411751176117711781179118011811182118311841185118611871188118911901191119211931194119511961197119811991200120112021203120412051206120712081209121012111212121312141215121612171218121912201221122212231224122512261227122812291230123112321233123412351236123712381239124012411242124312441245124612471248124912501251125212531254125512561257125812591260126112621263126412651266126712681269127012711272127312741275127612771278127912801281128212831284128512861287128812891290129112921293129412951296129712981299130013011302130313041305130613071308130913101311131213131314131513161317131813191320132113221323132413251326132713281329133013311332133313341335133613371338133913401341134213431344134513461347134813491350135113521353135413551356135713581359136013611362136313641365136613671368136913701371137213731374137513761377137813791380138113821383138413851386138713881389139013911392139313941395139613971398139914001401140214031404140514061407140814091410141114121413141414151416141714181419142014211422142314241425142614271428142914301431143214331434143514361437143814391440144114421443144414451446144714481449145014511452145314541455145614571458145914601461146214631464146514661467146814691470147114721473147414751476147714781479148014811482148314841485148614871488148914901491149214931494149514961497149814991500150115021503150415051506150715081509151015111512151315141515151615171518151915201521152215231524152515261527152815291530153115321533153415351536153715381539154015411542154315441545154615471548154915501551155215531554155515561557155815591560156115621563156415651566156715681569157015711572157315741575157615771578157915801581158215831584158515861587158815891590159115921593159415951596159715981599160016011602160316041605160616071608160916101611161216131614161516161617161816191620162116221623162416251626162716281629163016311632163316341635163616371638163916401641164216431644164516461647164816491650165116521653165416551656165716581659166016611662166316641665166616671668166916701671167216731674167516761677 |
- /*
- * This file is part of the UCB release of Plan 9. It is subject to the license
- * terms in the LICENSE file found in the top-level directory of this
- * distribution and at http://akaros.cs.berkeley.edu/files/Plan9License. No
- * part of the UCB release of Plan 9, including this file, may be copied,
- * modified, propagated, or distributed except according to the terms contained
- * in the LICENSE file.
- */
- #include <u.h>
- #include <libc.h>
- #include <bio.h>
- #include <draw.h>
- #include <event.h>
- #include "sky.h"
- #include "strings.C"
- enum
- {
- NNGC=7840, /* number of NGC numbers [1..NNGC] */
- NIC = 5386, /* number of IC numbers */
- NNGCrec=NNGC+NIC, /* number of records in the NGC catalog (including IC's, starting at NNGC */
- NMrec=122, /* number of M records */
- NM=110, /* number of M numbers */
- NAbell=2712, /* number of records in the Abell catalog */
- NName=1000, /* number of prose names; estimated maximum (read from editable text file) */
- NBayer=1517, /* number of bayer entries */
- NSAO=258998, /* number of SAO stars */
- MAXcon=1932, /* maximum number of patches in a constellation */
- Ncon=88, /* number of constellations */
- Npatch=92053, /* highest patch number */
- };
- char ngctype[NNGCrec];
- Mindexrec mindex[NMrec];
- Namerec name[NName];
- Bayerec bayer[NBayer];
- int32_t con[MAXcon];
- uint16_t conindex[Ncon+1];
- int32_t patchaddr[Npatch+1];
- Record *rec;
- Record *orec;
- Record *cur;
- char *dir=DIR;
- int saodb;
- int ngcdb;
- int abelldb;
- int ngctypedb;
- int mindexdb;
- int namedb;
- int bayerdb;
- int condb;
- int conindexdb;
- int patchdb;
- char parsed[3];
- int32_t nrec;
- int32_t nreca;
- int32_t norec;
- int32_t noreca;
- Biobuf bin;
- Biobuf bout;
- int
- main(int argc, char *argv[])
- {
- char *line;
- Binit(&bin, 0, OREAD);
- Binit(&bout, 1, OWRITE);
- if(argc != 1)
- dir = argv[1];
- astro("", 1);
- while((line = Brdline(&bin, '\n'))){
- line[Blinelen(&bin)-1] = 0;
- lookup(line, 1);
- Bflush(&bout);
- }
- if(display != nil){
- closedisplay(display);
- /* automatic refresh of rio window is triggered by mouse */
- close(open("/dev/mouse", OREAD));
- }
- return 0;
- }
- void
- reset(void)
- {
- nrec = 0;
- cur = rec;
- }
- void
- grow(void)
- {
- nrec++;
- if(nreca < nrec){
- nreca = nrec+50;
- rec = realloc(rec, nreca*sizeof(Record));
- if(rec == 0){
- fprint(2, "scat: realloc fails\n");
- exits("realloc");
- }
- }
- cur = rec+nrec-1;
- }
- void
- copy(void)
- {
- if(noreca < nreca){
- noreca = nreca;
- orec = realloc(orec, nreca*sizeof(Record));
- if(orec == 0){
- fprint(2, "scat: realloc fails\n");
- exits("realloc");
- }
- }
- memmove(orec, rec, nrec*sizeof(Record));
- norec = nrec;
- }
- int
- eopen(char *s)
- {
- char buf[128];
- int f;
- sprint(buf, "%s/%s.scat", dir, s);
- f = open(buf, 0);
- if(f<0){
- fprint(2, "scat: can't open %s\n", buf);
- exits("open");
- }
- return f;
- }
- void
- Eread(int f, char *name, void *addr, int32_t n)
- {
- if(read(f, addr, n) != n){ /* BUG! */
- fprint(2, "scat: read error on %s\n", name);
- exits("read");
- }
- }
- char*
- skipbl(char *s)
- {
- while(*s!=0 && (*s==' ' || *s=='\t'))
- s++;
- return s;
- }
- char*
- skipstr(char *s, char *t)
- {
- while(*s && *s==*t)
- s++, t++;
- return skipbl(s);
- }
- /* produce little-endian long at address l */
- int32_t
- Long(int32_t *l)
- {
- uint8_t *p;
- p = (uint8_t*)l;
- return (int32_t)p[0]|((int32_t)p[1]<<8)|((int32_t)p[2]<<16)|((int32_t)p[3]<<24);
- }
- /* produce little-endian long at address l */
- int
- Short(int16_t *s)
- {
- uint8_t *p;
- p = (uint8_t*)s;
- return p[0]|(p[1]<<8);
- }
- void
- nameopen(void)
- {
- Biobuf b;
- int i;
- char *l, *p;
- if(namedb == 0){
- namedb = eopen("name");
- Binit(&b, namedb, OREAD);
- for(i=0; i<NName; i++){
- l = Brdline(&b, '\n');
- if(l == 0)
- break;
- p = strchr(l, '\t');
- if(p == 0){
- Badformat:
- Bprint(&bout, "warning: name.scat bad format; line %d\n", i+1);
- break;
- }
- *p++ = 0;
- strcpy(name[i].name, l);
- if(strncmp(p, "ngc", 3) == 0)
- name[i].ngc = atoi(p+3);
- else if(strncmp(p, "ic", 2) == 0)
- name[i].ngc = atoi(p+2)+NNGC;
- else if(strncmp(p, "sao", 3) == 0)
- name[i].sao = atoi(p+3);
- else if(strncmp(p, "abell", 5) == 0)
- name[i].abell = atoi(p+5);
- else
- goto Badformat;
- }
- if(i == NName)
- Bprint(&bout, "warning: too many names in name.scat (max %d); extra ignored\n", NName);
- close(namedb);
- bayerdb = eopen("bayer");
- Eread(bayerdb, "bayer", bayer, sizeof bayer);
- close(bayerdb);
- for(i=0; i<NBayer; i++)
- bayer[i].sao = Long(&bayer[i].sao);
- }
- }
- void
- saoopen(void)
- {
- if(saodb == 0){
- nameopen();
- saodb = eopen("sao");
- }
- }
- void
- ngcopen(void)
- {
- if(ngcdb == 0){
- nameopen();
- ngcdb = eopen("ngc2000");
- ngctypedb = eopen("ngc2000type");
- Eread(ngctypedb, "ngctype", ngctype, sizeof ngctype);
- close(ngctypedb);
- }
- }
- void
- abellopen(void)
- {
- /* nothing extra to do with abell: it's directly indexed by number */
- if(abelldb == 0)
- abelldb = eopen("abell");
- }
- void
- patchopen(void)
- {
- Biobuf *b;
- int32_t l, m;
- char buf[100];
- if(patchdb == 0){
- patchdb = eopen("patch");
- sprint(buf, "%s/patchindex.scat", dir);
- b = Bopen(buf, OREAD);
- if(b == 0){
- fprint(2, "can't open %s\n", buf);
- exits("open");
- }
- for(m=0,l=0; l<=Npatch; l++)
- patchaddr[l] = m += Bgetc(b)*4;
- Bterm(b);
- }
- }
- void
- mopen(void)
- {
- int i;
- if(mindexdb == 0){
- mindexdb = eopen("mindex");
- Eread(mindexdb, "mindex", mindex, sizeof mindex);
- close(mindexdb);
- for(i=0; i<NMrec; i++)
- mindex[i].ngc = Short(&mindex[i].ngc);
- }
- }
- void
- constelopen(void)
- {
- int i;
- if(condb == 0){
- condb = eopen("con");
- conindexdb = eopen("conindex");
- Eread(conindexdb, "conindex", conindex, sizeof conindex);
- close(conindexdb);
- for(i=0; i<Ncon+1; i++)
- conindex[i] = Short((int16_t*)&conindex[i]);
- }
- }
- void
- lowercase(char *s)
- {
- for(; *s; s++)
- if('A'<=*s && *s<='Z')
- *s += 'a'-'A';
- }
- int
- loadngc(int32_t index)
- {
- static int failed;
- int32_t j;
- ngcopen();
- j = (index-1)*sizeof(NGCrec);
- grow();
- cur->type = NGC;
- cur->index = index;
- seek(ngcdb, j, 0);
- /* special case: NGC data may not be available */
- if(read(ngcdb, &cur->ngc, sizeof(NGCrec)) != sizeof(NGCrec)){
- if(!failed){
- fprint(2, "scat: NGC database not available\n");
- failed++;
- }
- cur->type = NONGC;
- cur->ngc.ngc = 0;
- cur->ngc.ra = 0;
- cur->ngc.dec = 0;
- cur->ngc.diam = 0;
- cur->ngc.mag = 0;
- return 0;
- }
- cur->ngc.ngc = Short(&cur->ngc.ngc);
- cur->ngc.ra = Long(&cur->ngc.ra);
- cur->ngc.dec = Long(&cur->ngc.dec);
- cur->ngc.diam = Long(&cur->ngc.diam);
- cur->ngc.mag = Short(&cur->ngc.mag);
- return 1;
- }
- int
- loadabell(int32_t index)
- {
- int32_t j;
- abellopen();
- j = index-1;
- grow();
- cur->type = Abell;
- cur->index = index;
- seek(abelldb, j*sizeof(Abellrec), 0);
- Eread(abelldb, "abell", &cur->abell, sizeof(Abellrec));
- cur->abell.abell = Short(&cur->abell.abell);
- if(cur->abell.abell != index){
- fprint(2, "bad format in abell catalog\n");
- exits("abell");
- }
- cur->abell.ra = Long(&cur->abell.ra);
- cur->abell.dec = Long(&cur->abell.dec);
- cur->abell.glat = Long(&cur->abell.glat);
- cur->abell.glong = Long(&cur->abell.glong);
- cur->abell.rad = Long(&cur->abell.rad);
- cur->abell.mag10 = Short(&cur->abell.mag10);
- cur->abell.pop = Short(&cur->abell.pop);
- cur->abell.dist = Short(&cur->abell.dist);
- return 1;
- }
- int
- loadsao(int index)
- {
- if(index<=0 || index>NSAO)
- return 0;
- saoopen();
- grow();
- cur->type = SAO;
- cur->index = index;
- seek(saodb, (index-1)*sizeof(SAOrec), 0);
- Eread(saodb, "sao", &cur->sao, sizeof(SAOrec));
- cur->sao.ra = Long(&cur->sao.ra);
- cur->sao.dec = Long(&cur->sao.dec);
- cur->sao.dra = Long(&cur->sao.dra);
- cur->sao.ddec = Long(&cur->sao.ddec);
- cur->sao.mag = Short(&cur->sao.mag);
- cur->sao.mpg = Short(&cur->sao.mpg);
- cur->sao.hd = Long(&cur->sao.hd);
- return 1;
- }
- int
- loadplanet(int index, Record *r)
- {
- if(index<0 || index>NPlanet || planet[index].name[0]=='\0')
- return 0;
- grow();
- cur->type = Planet;
- cur->index = index;
- /* check whether to take new or existing record */
- if(r == nil)
- memmove(&cur->planet, &planet[index], sizeof(Planetrec));
- else
- memmove(&cur->planet, &r->planet, sizeof(Planetrec));
- return 1;
- }
- int
- loadpatch(int32_t index)
- {
- int i;
- patchopen();
- if(index<=0 || index>Npatch)
- return 0;
- grow();
- cur->type = Patch;
- cur->index = index;
- seek(patchdb, patchaddr[index-1], 0);
- cur->patch.nkey = (patchaddr[index]-patchaddr[index-1])/4;
- Eread(patchdb, "patch", cur->patch.key, cur->patch.nkey*4);
- for(i=0; i<cur->patch.nkey; i++)
- cur->patch.key[i] = Long(&cur->patch.key[i]);
- return 1;
- }
- int
- loadtype(int t)
- {
- int i;
- ngcopen();
- for(i=0; i<NNGCrec; i++)
- if(t == (ngctype[i])){
- grow();
- cur->type = NGCN;
- cur->index = i+1;
- }
- return 1;
- }
- void
- flatten(void)
- {
- int i, j, notflat;
- Record *or;
- int32_t key;
- loop:
- copy();
- reset();
- notflat = 0;
- for(i=0,or=orec; i<norec; i++,or++){
- switch(or->type){
- default:
- fprint(2, "bad type %d in flatten\n", or->type);
- break;
- case NONGC:
- break;
- case Planet:
- case Abell:
- case NGC:
- case SAO:
- grow();
- memmove(cur, or, sizeof(Record));
- break;
- case NGCN:
- if(loadngc(or->index))
- notflat = 1;
- break;
- case NamedSAO:
- loadsao(or->index);
- notflat = 1;
- break;
- case NamedNGC:
- if(loadngc(or->index))
- notflat = 1;
- break;
- case NamedAbell:
- loadabell(or->index);
- notflat = 1;
- break;
- case PatchC:
- loadpatch(or->index);
- notflat = 1;
- break;
- case Patch:
- for(j=1; j<or->patch.nkey; j++){
- key = or->patch.key[j];
- if((key&0x3F) == SAO)
- loadsao((key>>8)&0xFFFFFF);
- else if((key&0x3F) == Abell)
- loadabell((key>>8)&0xFFFFFF);
- else
- loadngc((key>>16)&0xFFFF);
- }
- break;
- }
- }
- if(notflat)
- goto loop;
- }
- int
- ism(int index)
- {
- int i;
- for(i=0; i<NMrec; i++)
- if(mindex[i].ngc == index)
- return 1;
- return 0;
- }
- char*
- alpha(char *s, char *t)
- {
- int n;
- n = strlen(t);
- if(strncmp(s, t, n)==0 && (s[n]<'a' || 'z'<s[n]))
- return skipbl(s+n);
- return 0;
- }
- char*
- text(char *s, char *t)
- {
- int n;
- n = strlen(t);
- if(strncmp(s, t, n)==0 && (s[n]==0 || s[n]==' ' || s[n]=='\t'))
- return skipbl(s+n);
- return 0;
- }
- int
- cull(char *s, int keep, int dobbox)
- {
- int i, j, nobj, keepthis;
- Record *or;
- char *t;
- int dogrtr, doless, dom, dosao, dongc, doabell;
- int mgrtr, mless;
- char obj[100];
- memset(obj, 0, sizeof(obj));
- nobj = 0;
- dogrtr = 0;
- doless = 0;
- dom = 0;
- dongc = 0;
- dosao = 0;
- doabell = 0;
- mgrtr = mless= 0;
- if(dobbox)
- goto Cull;
- for(;;){
- if(s[0] == '>'){
- dogrtr = 1;
- mgrtr = 10 * strtod(s+1, &t);
- if(mgrtr==0 && t==s+1){
- fprint(2, "bad magnitude\n");
- return 0;
- }
- s = skipbl(t);
- continue;
- }
- if(s[0] == '<'){
- doless = 1;
- mless = 10 * strtod(s+1, &t);
- if(mless==0 && t==s+1){
- fprint(2, "bad magnitude\n");
- return 0;
- }
- s = skipbl(t);
- continue;
- }
- if((t = text(s, "m")) ){
- dom = 1;
- s = t;
- continue;
- }
- if((t = text(s, "sao"))){
- dosao = 1;
- s = t;
- continue;
- }
- if((t = text(s, "ngc"))){
- dongc = 1;
- s = t;
- continue;
- }
- if((t = text(s, "abell"))){
- doabell = 1;
- s = t;
- continue;
- }
- for(i=0; names[i].name; i++)
- if((t = alpha(s, names[i].name))){
- if(nobj > 100){
- fprint(2, "too many object types\n");
- return 0;
- }
- obj[nobj++] = names[i].type;
- s = t;
- goto Continue;
- }
- break;
- Continue:;
- }
- if(*s){
- fprint(2, "syntax error in object list\n");
- return 0;
- }
- Cull:
- flatten();
- copy();
- reset();
- if(dom)
- mopen();
- if(dosao)
- saoopen();
- if(dongc || nobj)
- ngcopen();
- if(doabell)
- abellopen();
- for(i=0,or=orec; i<norec; i++,or++){
- keepthis = !keep;
- if(dobbox && inbbox(or->ngc.ra, or->ngc.dec))
- keepthis = keep;
- if(doless && or->ngc.mag <= mless)
- keepthis = keep;
- if(dogrtr && or->ngc.mag >= mgrtr)
- keepthis = keep;
- if(dom && (or->type==NGC && ism(or->ngc.ngc)))
- keepthis = keep;
- if(dongc && or->type==NGC)
- keepthis = keep;
- if(doabell && or->type==Abell)
- keepthis = keep;
- if(dosao && or->type==SAO)
- keepthis = keep;
- for(j=0; j<nobj; j++)
- if(or->type==NGC && or->ngc.type==obj[j])
- keepthis = keep;
- if(keepthis){
- grow();
- memmove(cur, or, sizeof(Record));
- }
- }
- return 1;
- }
- int
- compar(const void *va, const void *vb)
- {
- Record const *a=va, *b=vb;
- if(a->type == b->type)
- return a->index - b->index;
- return a->type - b->type;
- }
- void
- sort(void)
- {
- int i;
- Record *r, *s;
- if(nrec == 0)
- return;
- qsort(rec, nrec, sizeof(Record), compar);
- r = rec+1;
- s = rec;
- for(i=1; i<nrec; i++,r++){
- /* may have multiple instances of a planet in the scene */
- if(r->type==s->type && r->index==s->index && r->type!=Planet)
- continue;
- memmove(++s, r, sizeof(Record));
- }
- nrec = (s+1)-rec;
- }
- char greekbuf[128];
- char*
- togreek(char *s)
- {
- char *t;
- int i, n;
- Rune r;
- t = greekbuf;
- while(*s){
- for(i=1; i<=24; i++){
- n = strlen(greek[i]);
- if(strncmp(s, greek[i], n)==0 && (s[n]==' ' || s[n]=='\t')){
- s += n;
- t += runetochar(t, &greeklet[i]);
- goto Cont;
- }
- }
- n = chartorune(&r, s);
- for(i=0; i<n; i++)
- *t++ = *s++;
- Cont:;
- }
- *t = 0;
- return greekbuf;
- }
- char*
- fromgreek(char *s)
- {
- char *t;
- int i, n;
- Rune r;
- t = greekbuf;
- while(*s){
- n = chartorune(&r, s);
- for(i=1; i<=24; i++){
- if(r == greeklet[i]){
- strcpy(t, greek[i]);
- t += strlen(greek[i]);
- s += n;
- goto Cont;
- }
- }
- for(i=0; i<n; i++)
- *t++ = *s++;
- Cont:;
- }
- *t = 0;
- return greekbuf;
- }
- #ifdef OLD
- /*
- * Old version
- */
- int
- coords(int deg)
- {
- int i;
- int x, y;
- Record *or;
- int32_t dec, ra, ndec, nra;
- int rdeg;
- flatten();
- copy();
- reset();
- deg *= 2;
- for(i=0,or=orec; i<norec; i++,or++){
- if(or->type == Planet) /* must keep it here */
- loadplanet(or->index, or);
- dec = or->ngc.dec/MILLIARCSEC;
- ra = or->ngc.ra/MILLIARCSEC;
- rdeg = deg/cos((dec*PI)/180);
- for(y=-deg; y<=+deg; y++){
- ndec = dec*2+y;
- if(ndec/2>=90 || ndec/2<=-90)
- continue;
- /* fp errors hurt here, so we round 1' to the pole */
- if(ndec >= 0)
- ndec = ndec*500*60*60 + 60000;
- else
- ndec = ndec*500*60*60 - 60000;
- for(x=-rdeg; x<=+rdeg; x++){
- nra = ra*2+x;
- if(nra/2 < 0)
- nra += 360*2;
- if(nra/2 >= 360)
- nra -= 360*2;
- /* fp errors hurt here, so we round up 1' */
- nra = nra/2*MILLIARCSEC + 60000;
- loadpatch(patcha(angle(nra), angle(ndec)));
- }
- }
- }
- sort();
- return 1;
- }
- #endif
- /*
- * New version attempts to match the boundaries of the plot better.
- */
- int
- coords(int deg)
- {
- int i;
- int x, y, xx;
- Record *or;
- int32_t min, circle;
- double factor;
- flatten();
- circle = 360*MILLIARCSEC;
- deg *= MILLIARCSEC;
- /* find center */
- folded = 0;
- bbox(0, 0, 0);
- /* now expand */
- factor = cos(angle((decmax+decmin)/2));
- if(factor < .2)
- factor = .2;
- factor = floor(1/factor);
- folded = 0;
- bbox(factor*deg, deg, 1);
- Bprint(&bout, "%s to ", hms(angle(ramin)));
- Bprint(&bout, "%s\n", hms(angle(ramax)));
- Bprint(&bout, "%s to ", dms(angle(decmin)));
- Bprint(&bout, "%s\n", dms(angle(decmax)));
- copy();
- reset();
- for(i=0,or=orec; i<norec; i++,or++)
- if(or->type == Planet) /* must keep it here */
- loadplanet(or->index, or);
- min = ramin;
- if(ramin > ramax)
- min -= circle;
- for(x=min; x<=ramax; x+=250*60*60){
- xx = x;
- if(xx < 0)
- xx += circle;
- for(y=decmin; y<=decmax; y+=250*60*60)
- if(-circle/4 < y && y < circle/4)
- loadpatch(patcha(angle(xx), angle(y)));
- }
- sort();
- cull(nil, 1, 1);
- return 1;
- }
- void
- pplate(char *flags)
- {
- int i;
- int32_t c;
- int na, rah, ram, d1, d2;
- double r0;
- int ra, dec;
- int32_t ramin, ramax, decmin, decmax; /* all in degrees */
- Record *r;
- int folded;
- Angle racenter, deccenter, rasize, decsize, a[4];
- Picture *pic;
- rasize = -1.0;
- decsize = -1.0;
- na = 0;
- for(;;){
- while(*flags==' ')
- flags++;
- if(('0'<=*flags && *flags<='9') || *flags=='+' || *flags=='-'){
- if(na >= 3)
- goto err;
- a[na++] = getra(flags);
- while(*flags && *flags!=' ')
- flags++;
- continue;
- }
- if(*flags){
- err:
- Bprint(&bout, "syntax error in plate\n");
- return;
- }
- break;
- }
- switch(na){
- case 0:
- break;
- case 1:
- rasize = a[0];
- decsize = rasize;
- break;
- case 2:
- rasize = a[0];
- decsize = a[1];
- break;
- case 3:
- case 4:
- racenter = a[0];
- deccenter = a[1];
- rasize = a[2];
- if(na == 4)
- decsize = a[3];
- else
- decsize = rasize;
- if(rasize<0.0 || decsize<0.0){
- Bprint(&bout, "negative sizes\n");
- return;
- }
- goto done;
- }
- folded = 0;
- /* convert to milliarcsec */
- c = 1000*60*60;
- Again:
- if(nrec == 0){
- Bprint(&bout, "empty\n");
- return;
- }
- ramin = 0x7FFFFFFF;
- ramax = -0x7FFFFFFF;
- decmin = 0x7FFFFFFF;
- decmax = -0x7FFFFFFF;
- for(r=rec,i=0; i<nrec; i++,r++){
- if(r->type == Patch){
- radec(r->index, &rah, &ram, &dec);
- ra = 15*rah+ram/4;
- r0 = c/cos(RAD(dec));
- ra *= c;
- dec *= c;
- if(dec == 0)
- d1 = c, d2 = c;
- else if(dec < 0)
- d1 = c, d2 = 0;
- else
- d1 = 0, d2 = c;
- }else if(r->type==SAO || r->type==NGC || r->type==Abell){
- ra = r->ngc.ra;
- dec = r->ngc.dec;
- d1 = 0, d2 = 0, r0 = 0;
- }else if(r->type==NGCN){
- loadngc(r->index);
- continue;
- }else if(r->type==NamedSAO){
- loadsao(r->index);
- continue;
- }else if(r->type==NamedNGC){
- loadngc(r->index);
- continue;
- }else if(r->type==NamedAbell){
- loadabell(r->index);
- continue;
- }else
- continue;
- if(dec+d2 > decmax)
- decmax = dec+d2;
- if(dec-d1 < decmin)
- decmin = dec-d1;
- if(folded){
- ra -= 180*c;
- if(ra < 0)
- ra += 360*c;
- }
- if(ra+r0 > ramax)
- ramax = ra+r0;
- if(ra < ramin)
- ramin = ra;
- }
- if(!folded && ramax-ramin>270*c){
- folded = 1;
- goto Again;
- }
- racenter = angle(ramin+(ramax-ramin)/2);
- deccenter = angle(decmin+(decmax-decmin)/2);
- if(rasize<0 || decsize<0){
- rasize = angle(ramax-ramin)*cos(deccenter);
- decsize = angle(decmax-decmin);
- }
- done:
- if(DEG(rasize)>1.1 || DEG(decsize)>1.1){
- Bprint(&bout, "plate too big: %s", ms(rasize));
- Bprint(&bout, " x %s\n", ms(decsize));
- Bprint(&bout, "trimming to 30'x30'\n");
- rasize = RAD(0.5);
- decsize = RAD(0.5);
- }
- Bprint(&bout, "%s %s ", hms(racenter), dms(deccenter));
- Bprint(&bout, "%s", ms(rasize));
- Bprint(&bout, " x %s\n", ms(decsize));
- Bflush(&bout);
- flatten();
- pic = image(racenter, deccenter, rasize, decsize);
- if(pic == 0)
- return;
- Bprint(&bout, "plate %s locn %d %d %d %d\n", pic->name, pic->minx, pic->miny, pic->maxx, pic->maxy);
- Bflush(&bout);
- displaypic(pic);
- }
- void
- lookup(char *s, int doreset)
- {
- int i, j, k;
- int rah, ram, deg;
- char *starts, *inputline=s, *t, *u;
- Record *r;
- int32_t n;
- double x;
- Angle ra;
- lowercase(s);
- s = skipbl(s);
- if(*s == 0)
- goto Print;
- if((t = alpha(s, "flat"))){
- if(*t){
- fprint(2, "flat takes no arguments\n");
- return;
- }
- if(nrec == 0){
- fprint(2, "no records\n");
- return;
- }
- flatten();
- goto Print;
- }
- if((t = alpha(s, "print"))){
- if(*t){
- fprint(2, "print takes no arguments\n");
- return;
- }
- for(i=0,r=rec; i<nrec; i++,r++)
- prrec(r);
- return;
- }
- if((t = alpha(s, "add"))){
- lookup(t, 0);
- return;
- }
- if((t = alpha(s, "sao"))){
- n = strtoul(t, &u, 10);
- if(n<=0 || n>NSAO)
- goto NotFound;
- t = skipbl(u);
- if(*t){
- fprint(2, "syntax error in sao\n");
- return;
- }
- if(doreset)
- reset();
- if(!loadsao(n))
- goto NotFound;
- goto Print;
- }
- if((t = alpha(s, "ngc"))){
- n = strtoul(t, &u, 10);
- if(n<=0 || n>NNGC)
- goto NotFound;
- t = skipbl(u);
- if(*t){
- fprint(2, "syntax error in ngc\n");
- return;
- }
- if(doreset)
- reset();
- if(!loadngc(n))
- goto NotFound;
- goto Print;
- }
- if((t = alpha(s, "ic"))){
- n = strtoul(t, &u, 10);
- if(n<=0 || n>NIC)
- goto NotFound;
- t = skipbl(u);
- if(*t){
- fprint(2, "syntax error in ic\n");
- return;
- }
- if(doreset)
- reset();
- if(!loadngc(n+NNGC))
- goto NotFound;
- goto Print;
- }
- if((t = alpha(s, "abell"))){
- n = strtoul(t, &u, 10);
- if(n<=0 || n>NAbell)
- goto NotFound;
- if(doreset)
- reset();
- if(!loadabell(n))
- goto NotFound;
- goto Print;
- }
- if((t = alpha(s, "m"))){
- n = strtoul(t, &u, 10);
- if(n<=0 || n>NM)
- goto NotFound;
- mopen();
- for(j=n-1; mindex[j].m<n; j++)
- ;
- if(doreset)
- reset();
- while(mindex[j].m == n){
- if(mindex[j].ngc){
- grow();
- cur->type = NGCN;
- cur->index = mindex[j].ngc;
- }
- j++;
- }
- goto Print;
- }
- for(i=1; i<=Ncon; i++)
- if((t = alpha(s, constel[i]))){
- if(*t){
- fprint(2, "syntax error in constellation\n");
- return;
- }
- constelopen();
- seek(condb, 4L*conindex[i-1], 0);
- j = conindex[i]-conindex[i-1];
- Eread(condb, "con", con, 4*j);
- if(doreset)
- reset();
- for(k=0; k<j; k++){
- grow();
- cur->type = PatchC;
- cur->index = Long(&con[k]);
- }
- goto Print;
- }
- if((t = alpha(s, "expand"))){
- n = 0;
- if(*t){
- if(*t<'0' && '9'<*t){
- Expanderr:
- fprint(2, "syntax error in expand\n");
- return;
- }
- n = strtoul(t, &u, 10);
- t = skipbl(u);
- if(*t)
- goto Expanderr;
- }
- coords(n);
- goto Print;
- }
- if((t = alpha(s, "plot"))){
- if(nrec == 0){
- Bprint(&bout, "empty\n");
- return;
- }
- plot(t);
- return;
- }
- if((t = alpha(s, "astro"))){
- astro(t, 0);
- return;
- }
- if((t = alpha(s, "plate"))){
- pplate(t);
- return;
- }
- if((t = alpha(s, "gamma"))){
- while(*t==' ')
- t++;
- u = t;
- x = strtod(t, &u);
- if(u > t)
- gam.gamma = x;
- Bprint(&bout, "%.2f\n", gam.gamma);
- return;
- }
- if((t = alpha(s, "keep"))){
- if(!cull(t, 1, 0))
- return;
- goto Print;
- }
- if((t = alpha(s, "drop"))){
- if(!cull(t, 0, 0))
- return;
- goto Print;
- }
- for(i=0; planet[i].name[0]; i++){
- if((t = alpha(s, planet[i].name))){
- if(doreset)
- reset();
- loadplanet(i, nil);
- goto Print;
- }
- }
- for(i=0; names[i].name; i++){
- if((t = alpha(s, names[i].name))){
- if(*t){
- fprint(2, "syntax error in type\n");
- return;
- }
- if(doreset)
- reset();
- loadtype(names[i].type);
- goto Print;
- }
- }
- switch(s[0]){
- case '"':
- starts = ++s;
- while(*s != '"')
- if(*s++ == 0){
- fprint(2, "bad star name\n");
- return;
- }
- *s = 0;
- if(doreset)
- reset();
- j = nrec;
- saoopen();
- starts = fromgreek(starts);
- for(i=0; i<NName; i++)
- if(equal(starts, name[i].name)){
- grow();
- if(name[i].sao){
- rec[j].type = NamedSAO;
- rec[j].index = name[i].sao;
- }
- if(name[i].ngc){
- rec[j].type = NamedNGC;
- rec[j].index = name[i].ngc;
- }
- if(name[i].abell){
- rec[j].type = NamedAbell;
- rec[j].index = name[i].abell;
- }
- strcpy(rec[j].named.name, name[i].name);
- j++;
- }
- if(parsename(starts))
- for(i=0; i<NBayer; i++)
- if(bayer[i].name[0]==parsed[0] &&
- (bayer[i].name[1]==parsed[1] || parsed[1]==0) &&
- bayer[i].name[2]==parsed[2]){
- grow();
- rec[j].type = NamedSAO;
- rec[j].index = bayer[i].sao;
- strncpy(rec[j].named.name, starts, sizeof(rec[j].named.name));
- j++;
- }
- if(j == 0){
- *s = '"';
- goto NotFound;
- }
- break;
- case '0': case '1': case '2': case '3': case '4':
- case '5': case '6': case '7': case '8': case '9':
- strtoul(s, &t, 10);
- if(*t != 'h'){
- BadCoords:
- fprint(2, "bad coordinates %s\n", inputline);
- break;
- }
- ra = DEG(getra(s));
- while(*s && *s!=' ' && *s!='\t')
- s++;
- rah = ra/15;
- ra = ra-rah*15;
- ram = ra*4;
- deg = strtol(s, &t, 10);
- if(t == s)
- goto BadCoords;
- /* degree sign etc. is optional */
- if((uint8_t)*t == L'°')
- deg = DEG(getra(s));
- if(doreset)
- reset();
- if(abs(deg)>=90 || rah>=24)
- goto BadCoords;
- if(!loadpatch(patch(rah, ram, deg)))
- goto NotFound;
- break;
- default:
- fprint(2, "unknown command %s\n", inputline);
- return;
- }
- Print:
- if(nrec == 0)
- Bprint(&bout, "empty\n");
- else if(nrec <= 2)
- for(i=0; i<nrec; i++)
- prrec(rec+i);
- else
- Bprint(&bout, "%ld items\n", nrec);
- return;
- NotFound:
- fprint(2, "%s not found\n", inputline);
- return;
- }
- char *ngctypes[] =
- {
- [Galaxy] = "Gx",
- [PlanetaryN] = "Pl",
- [OpenCl] = "OC",
- [GlobularCl] = "Gb",
- [DiffuseN] = "Nb",
- [NebularCl] = "C+N",
- [Asterism] = "Ast",
- [Knot] = "Kt",
- [Triple] = "***",
- [Double] = "D*",
- [Single] = "*",
- [Uncertain] = "?",
- [Nonexistent] = "-",
- [Unknown] = " ",
- [PlateDefect] = "PD",
- };
- char*
- ngcstring(int d)
- {
- if(d<Galaxy || d>PlateDefect)
- return "can't happen";
- return ngctypes[d];
- }
- int16_t descindex[NINDEX];
- void
- printnames(Record *r)
- {
- int i, ok, done;
- done = 0;
- for(i=0; i<NName; i++){ /* stupid linear search! */
- ok = 0;
- if(r->type==SAO && r->index==name[i].sao)
- ok = 1;
- if(r->type==NGC && r->ngc.ngc==name[i].ngc)
- ok = 1;
- if(r->type==Abell && r->abell.abell==name[i].abell)
- ok = 1;
- if(ok){
- if(done++ == 0)
- Bprint(&bout, "\t");
- Bprint(&bout, " \"%s\"", togreek(name[i].name));
- }
- }
- if(done)
- Bprint(&bout, "\n");
- }
- int
- equal(char *s1, char *s2)
- {
- int c;
- while(*s1){
- if(*s1==' '){
- while(*s1==' ')
- s1++;
- continue;
- }
- while(*s2==' ')
- s2++;
- c=*s2;
- if('A'<=*s2 && *s2<='Z')
- c^=' ';
- if(*s1!=c)
- return 0;
- s1++, s2++;
- }
- return 1;
- }
- int
- parsename(char *s)
- {
- char *blank;
- int i;
- blank = strchr(s, ' ');
- if(blank==0 || strchr(blank+1, ' ') || strlen(blank+1)!=3)
- return 0;
- blank++;
- parsed[0] = parsed[1] = parsed[2] = 0;
- if('0'<=s[0] && s[0]<='9'){
- i = atoi(s);
- parsed[0] = i;
- if(i > 100)
- return 0;
- }else{
- for(i=1; i<=24; i++)
- if(strncmp(greek[i], s, strlen(greek[i]))==0){
- parsed[0]=100+i;
- goto out;
- }
- return 0;
- out:
- if('0'<=s[strlen(greek[i])] && s[strlen(greek[i])]<='9')
- parsed[1]=s[strlen(greek[i])]-'0';
- }
- for(i=1; i<=88; i++)
- if(strcmp(constel[i], blank)==0){
- parsed[2] = i;
- return 1;
- }
- return 0;
- }
- char*
- dist_grp(int dg)
- {
- switch(dg){
- default:
- return "unknown";
- case 1:
- return "13.3-14.0";
- case 2:
- return "14.1-14.8";
- case 3:
- return "14.9-15.6";
- case 4:
- return "15.7-16.4";
- case 5:
- return "16.5-17.2";
- case 6:
- return "17.3-18.0";
- case 7:
- return ">18.0";
- }
- }
- char*
- rich_grp(int dg)
- {
- switch(dg){
- default:
- return "unknown";
- case 0:
- return "30-40";
- case 1:
- return "50-79";
- case 2:
- return "80-129";
- case 3:
- return "130-199";
- case 4:
- return "200-299";
- case 5:
- return ">=300";
- }
- }
- char*
- nameof(Record *r)
- {
- NGCrec *n;
- SAOrec *s;
- Abellrec *a;
- static char buf[128];
- int i;
- switch(r->type){
- default:
- return nil;
- case SAO:
- s = &r->sao;
- if(s->name[0] == 0)
- return nil;
- if(s->name[0] >= 100){
- i = snprint(buf, sizeof buf, "%C", greeklet[s->name[0]-100]);
- if(s->name[1])
- i += snprint(buf+i, sizeof buf-i, "%d", s->name[1]);
- }else
- i = snprint(buf, sizeof buf, " %d", s->name[0]);
- snprint(buf+i, sizeof buf-i, " %s", constel[(uintptr_t)s->name[2]]);
- break;
- case NGC:
- n = &r->ngc;
- if(n->type >= Uncertain)
- return nil;
- if(n->ngc <= NNGC)
- snprint(buf, sizeof buf, "NGC%4d ", n->ngc);
- else
- snprint(buf, sizeof buf, "IC%4d ", n->ngc-NNGC);
- break;
- case Abell:
- a = &r->abell;
- snprint(buf, sizeof buf, "Abell%4d", a->abell);
- break;
- }
- return buf;
- }
- void
- prrec(Record *r)
- {
- NGCrec *n;
- SAOrec *s;
- Abellrec *a;
- Planetrec *p;
- int i, rah, ram, dec, nn;
- int32_t key;
- if(r) switch(r->type){
- default:
- fprint(2, "can't prrec type %d\n", r->type);
- exits("type");
- case Planet:
- p = &r->planet;
- Bprint(&bout, "%s", p->name);
- Bprint(&bout, "\t%s %s",
- hms(angle(p->ra)),
- dms(angle(p->dec)));
- Bprint(&bout, " %3.2f° %3.2f°",
- p->az/(double)MILLIARCSEC, p->alt/(double)MILLIARCSEC);
- Bprint(&bout, " %s",
- ms(angle(p->semidiam)));
- if(r->index <= 1)
- Bprint(&bout, " %g", p->phase);
- Bprint(&bout, "\n");
- break;
- case NGC:
- n = &r->ngc;
- if(n->ngc <= NNGC)
- Bprint(&bout, "NGC%4d ", n->ngc);
- else
- Bprint(&bout, "IC%4d ", n->ngc-NNGC);
- Bprint(&bout, "%s ", ngcstring(n->type));
- if(n->mag == UNKNOWNMAG)
- Bprint(&bout, "----");
- else
- Bprint(&bout, "%.1f%c", n->mag/10.0, n->magtype);
- Bprint(&bout, "\t%s %s\t%c%.1f'\n",
- hm(angle(n->ra)),
- dm(angle(n->dec)),
- n->diamlim,
- DEG(angle(n->diam))*60.);
- prdesc(n->desc, desctab, descindex);
- printnames(r);
- break;
- case Abell:
- a = &r->abell;
- Bprint(&bout, "Abell%4d %.1f %.2f° %dMpc", a->abell, a->mag10/10.0,
- DEG(angle(a->rad)), a->dist);
- Bprint(&bout, "\t%s %s\t%.2f %.2f\n",
- hm(angle(a->ra)),
- dm(angle(a->dec)),
- DEG(angle(a->glat)),
- DEG(angle(a->glong)));
- Bprint(&bout, "\tdist grp: %s rich grp: %s %d galaxies/°²\n",
- dist_grp(a->distgrp),
- rich_grp(a->richgrp),
- a->pop);
- printnames(r);
- break;
- case SAO:
- s = &r->sao;
- Bprint(&bout, "SAO%6ld ", r->index);
- if(s->mag==UNKNOWNMAG)
- Bprint(&bout, "---");
- else
- Bprint(&bout, "%.1f", s->mag/10.0);
- if(s->mpg==UNKNOWNMAG)
- Bprint(&bout, ",---");
- else
- Bprint(&bout, ",%.1f", s->mpg/10.0);
- Bprint(&bout, " %s %s %.4fs %.3f\"",
- hms(angle(s->ra)),
- dms(angle(s->dec)),
- DEG(angle(s->dra))*(4*60),
- DEG(angle(s->ddec))*(60*60));
- Bprint(&bout, " %.3s %c %.2s %ld %d",
- s->spec, s->code, s->compid, s->hd, s->hdcode);
- if(s->name[0])
- Bprint(&bout, " \"%s\"", nameof(r));
- Bprint(&bout, "\n");
- printnames(r);
- break;
- case Patch:
- radec(r->index, &rah, &ram, &dec);
- Bprint(&bout, "%dh%dm %d°", rah, ram, dec);
- key = r->patch.key[0];
- Bprint(&bout, " %s", constel[key&0xFF]);
- if((key>>=8) & 0xFF)
- Bprint(&bout, " %s", constel[key&0xFF]);
- if((key>>=8) & 0xFF)
- Bprint(&bout, " %s", constel[key&0xFF]);
- if((key>>=8) & 0xFF)
- Bprint(&bout, " %s", constel[key&0xFF]);
- for(i=1; i<r->patch.nkey; i++){
- key = r->patch.key[i];
- switch(key&0x3F){
- case SAO:
- Bprint(&bout, " SAO%ld", (key>>8)&0xFFFFFF);
- break;
- case Abell:
- Bprint(&bout, " Abell%ld", (key>>8)&0xFFFFFF);
- break;
- default: /* NGC */
- nn = (key>>16)&0xFFFF;
- if(nn > NNGC)
- Bprint(&bout, " IC%d", nn-NNGC);
- else
- Bprint(&bout, " NGC%d", nn);
- Bprint(&bout, "(%s)", ngcstring(key&0x3F));
- break;
- }
- }
- Bprint(&bout, "\n");
- break;
- case NGCN:
- if(r->index <= NNGC)
- Bprint(&bout, "NGC%ld\n", r->index);
- else
- Bprint(&bout, "IC%ld\n", r->index-NNGC);
- break;
- case NamedSAO:
- Bprint(&bout, "SAO%ld \"%s\"\n", r->index, togreek(r->named.name));
- break;
- case NamedNGC:
- if(r->index <= NNGC)
- Bprint(&bout, "NGC%ld \"%s\"\n", r->index, togreek(r->named.name));
- else
- Bprint(&bout, "IC%ld \"%s\"\n", r->index-NNGC, togreek(r->named.name));
- break;
- case NamedAbell:
- Bprint(&bout, "Abell%ld \"%s\"\n", r->index, togreek(r->named.name));
- break;
- case PatchC:
- radec(r->index, &rah, &ram, &dec);
- Bprint(&bout, "%dh%dm %d\n", rah, ram, dec);
- break;
- }
- }
|