12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394959697989910010110210310410510610710810911011111211311411511611711811912012112212312412512612712812913013113213313413513613713813914014114214314414514614714814915015115215315415515615715815916016116216316416516616716816917017117217317417517617717817918018118218318418518618718818919019119219319419519619719819920020120220320420520620720820921021121221321421521621721821922022122222322422522622722822923023123223323423523623723823924024124224324424524624724824925025125225325425525625725825926026126226326426526626726826927027127227327427527627727827928028128228328428528628728828929029129229329429529629729829930030130230330430530630730830931031131231331431531631731831932032132232332432532632732832933033133233333433533633733833934034134234334434534634734834935035135235335435535635735835936036136236336436536636736836937037137237337437537637737837938038138238338438538638738838939039139239339439539639739839940040140240340440540640740840941041141241341441541641741841942042142242342442542642742842943043143243343443543643743843944044144244344444544644744844945045145245345445545645745845946046146246346446546646746846947047147247347447547647747847948048148248348448548648748848949049149249349449549649749849950050150250350450550650750850951051151251351451551651751851952052152252352452552652752852953053153253353453553653753853954054154254354454554654754854955055155255355455555655755855956056156256356456556656756856957057157257357457557657757857958058158258358458558658758858959059159259359459559659759859960060160260360460560660760860961061161261361461561661761861962062162262362462562662762862963063163263363463563663763863964064164264364464564664764864965065165265365465565665765865966066166266366466566666766866967067167267367467567667767867968068168268368468568668768868969069169269369469569669769869970070170270370470570670770870971071171271371471571671771871972072172272372472572672772872973073173273373473573673773873974074174274374474574674774874975075175275375475575675775875976076176276376476576676776876977077177277377477577677777877978078178278378478578678778878979079179279379479579679779879980080180280380480580680780880981081181281381481581681781881982082182282382482582682782882983083183283383483583683783883984084184284384484584684784884985085185285385485585685785885986086186286386486586686786886987087187287387487587687787887988088188288388488588688788888989089189289389489589689789889990090190290390490590690790890991091191291391491591691791891992092192292392492592692792892993093193293393493593693793893994094194294394494594694794894995095195295395495595695795895996096196296396496596696796896997097197297397497597697797897998098198298398498598698798898999099199299399499599699799899910001001100210031004100510061007100810091010101110121013101410151016101710181019102010211022102310241025102610271028102910301031103210331034103510361037103810391040104110421043104410451046104710481049105010511052105310541055105610571058105910601061106210631064106510661067106810691070107110721073107410751076107710781079108010811082108310841085108610871088108910901091109210931094109510961097109810991100110111021103110411051106110711081109111011111112111311141115111611171118111911201121112211231124112511261127112811291130113111321133113411351136113711381139114011411142114311441145114611471148114911501151115211531154115511561157115811591160116111621163116411651166116711681169117011711172117311741175117611771178117911801181118211831184118511861187118811891190119111921193119411951196119711981199120012011202120312041205120612071208120912101211121212131214121512161217121812191220122112221223122412251226122712281229123012311232123312341235123612371238123912401241124212431244124512461247124812491250125112521253125412551256125712581259126012611262126312641265126612671268126912701271127212731274127512761277127812791280128112821283128412851286128712881289129012911292129312941295129612971298129913001301130213031304130513061307130813091310131113121313131413151316131713181319132013211322132313241325132613271328132913301331133213331334133513361337133813391340134113421343134413451346134713481349135013511352135313541355135613571358135913601361136213631364136513661367136813691370137113721373137413751376137713781379138013811382138313841385138613871388138913901391139213931394139513961397139813991400140114021403140414051406140714081409141014111412141314141415141614171418141914201421142214231424142514261427142814291430143114321433143414351436143714381439144014411442144314441445144614471448144914501451145214531454145514561457145814591460146114621463146414651466146714681469147014711472147314741475147614771478147914801481148214831484148514861487148814891490149114921493149414951496149714981499150015011502150315041505150615071508150915101511151215131514151515161517151815191520152115221523152415251526152715281529153015311532153315341535153615371538153915401541154215431544154515461547154815491550155115521553155415551556155715581559156015611562156315641565156615671568156915701571157215731574157515761577157815791580158115821583158415851586158715881589159015911592159315941595159615971598159916001601160216031604160516061607160816091610161116121613161416151616161716181619162016211622162316241625162616271628162916301631163216331634163516361637163816391640164116421643164416451646164716481649165016511652165316541655165616571658165916601661166216631664166516661667166816691670167116721673167416751676 |
- /*
- * 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;
- 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 *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[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;
- }
- }
|