/*  * 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 "sky.h" double	PI_180	= 0.0174532925199432957692369;double	TWOPI	= 6.2831853071795864769252867665590057683943387987502;double	LN2	= 0.69314718055994530941723212145817656807550013436025;static double angledangle=(180./PI)*MILLIARCSEC; intrint(char *p, int n){	int i=0; 	while(*p==' ' && n)		p++, --n;	while(n--)		i=i*10+*p++-'0';	return i;} DAngledangle(Angle angle){	return angle*angledangle;} Angleangle(DAngle dangle){	return dangle/angledangle;} doublerfloat(char *p, int n){	double i, d=0; 	while(*p==' ' && n)		p++, --n;	if(*p == '+')		return rfloat(p+1, n-1);	if(*p == '-')		return -rfloat(p+1, n-1);	while(*p == ' ' && n)		p++, --n;	if(n == 0)		return 0.0;	while(n-- && *p!='.')		d = d*10+*p++-'0';	if(n <= 0)		return d;	p++;	i = 1;	while(n--)		d+=(*p++-'0')/(i*=10.);	return d;} intsign(int c){	if(c=='-')		return -1;	return 1;} char*hms(Angle a){	static char buf[20];	double x;	int h, m, s, ts; 	x=DEG(a)/15;	x += 0.5/36000.;	/* round up half of 0.1 sec */	h = floor(x);	x -= h;	x *= 60;	m = floor(x);	x -= m;	x *= 60;	s = floor(x);	x -= s;	ts = 10*x;	sprint(buf, "%dh%.2dm%.2d.%ds", h, m, s, ts);	return buf;} char*dms(Angle a){	static char buf[20];	double x;	int sign, d, m, s, ts; 	x = DEG(a);	sign='+';	if(a<0){		sign='-';		x=-x;	}	x += 0.5/36000.;	/* round up half of 0.1 arcsecond */	d = floor(x);	x -= d;	x *= 60;	m = floor(x);	x -= m;	x *= 60;	s = floor(x);	x -= s;	ts = floor(10*x);	sprint(buf, "%c%d°%.2d'%.2d.%d\"", sign, d, m, s, ts);	return buf;} char*ms(Angle a){	static char buf[20];	double x;	int d, m, s, ts; 	x = DEG(a);	x += 0.5/36000.;	/* round up half of 0.1 arcsecond */	d = floor(x);	x -= d;	x *= 60;	m = floor(x);	x -= m;	x *= 60;	s = floor(x);	x -= s;	ts = floor(10*x);	if(d != 0)		sprint(buf, "%d°%.2d'%.2d.%d\"", d, m, s, ts);	else		sprint(buf, "%.2d'%.2d.%d\"", m, s, ts);	return buf;} char*hm(Angle a){	static char buf[20];	double x;	int h, m, n; 	x = DEG(a)/15;	x += 0.5/600.;	/* round up half of tenth of minute */	h = floor(x);	x -= h;	x *= 60;	m = floor(x);	x -= m;	x *= 10;	n = floor(x);	sprint(buf, "%dh%.2d.%1dm", h, m, n);	return buf;} char*hm5(Angle a){	static char buf[20];	double x;	int h, m; 	x = DEG(a)/15;	x += 2.5/60.;	/* round up 2.5m */	h = floor(x);	x -= h;	x *= 60;	m = floor(x);	m -= m % 5;	sprint(buf, "%dh%.2dm", h, m);	return buf;} char*dm(Angle a){	static char buf[20];	double x;	int sign, d, m, n; 	x = DEG(a);	sign='+';	if(a<0){		sign='-';		x=-x;	}	x += 0.5/600.;	/* round up half of tenth of arcminute */	d = floor(x);	x -= d;	x *= 60;	m = floor(x);	x -= m;	x *= 10;	n = floor(x);	sprint(buf, "%c%d°%.2d.%.1d'", sign, d, m, n);	return buf;} char*deg(Angle a){	static char buf[20];	double x;	int sign, d; 	x = DEG(a);	sign='+';	if(a<0){		sign='-';		x=-x;	}	x += 0.5;	/* round up half degree */	d = floor(x);	sprint(buf, "%c%d°", sign, d);	return buf;} char*getword(char *ou, char *in){	int c; 	for(;;) {		c = *in++;		if(c == ' ' || c == '\t')			continue;		if(c == 0)			return 0;		break;	} 	if(c == '\'')		for(;;) {			if(c >= 'A' && c <= 'Z')				c += 'a' - 'A';			*ou++ = c;			c = *in++;			if(c == 0)				return 0;			if(c == '\'') {				*ou = 0;				return in-1;			}		}	for(;;) {		if(c >= 'A' && c <= 'Z')			c += 'a' - 'A';		*ou++ = c;		c = *in++;		if(c == ' ' || c == '\t' || c == 0) {			*ou = 0;			return in-1;		}	}} /* * Read formatted angle.  Must contain no embedded blanks */Anglegetra(char *p){	Rune r;	char *q;	Angle f, d;	int neg; 	neg = 0;	d = 0;	while(*p == ' ')		p++;	for(;;) {		if(*p == ' ' || *p=='\0')			goto Return;		if(*p == '-') {			neg = 1;			p++;		}		if(*p == '+') {			neg = 0;			p++;		}		q = p;		f = strtod(p, &q);		if(q > p) {			p = q;		}		p += chartorune(&r, p);		switch(r) {		default:		Return:			if(neg)				d = -d;			return RAD(d);		case 'h':			d += f*15;			break;		case 'm':			d += f/4;			break;		case 's':			d += f/240;			break;		case L'°':			d += f;			break;		case '\'':			d += f/60;			break;		case '\"':			d += f/3600;			break;		}	}} doublexsqrt(double a){ 	if(a < 0)		return 0;	return sqrt(a);} Angledist(Angle ra1, Angle dec1, Angle ra2, Angle dec2){	double a; 	a = sin(dec1) * sin(dec2) +		cos(dec1) * cos(dec2) *		cos(ra1 - ra2);	a = atan2(xsqrt(1 - a*a), a);	if(a < 0)		a = -a;	return a;} intdogamma(Pix c){	float f; 	f = c - gam.min;	if(f < 1)		f = 1; 	if(gam.absgamma == 1)		c = f * gam.mult2;	else		c = exp(log(f*gam.mult1) * gam.absgamma) * 255;	if(c > 255)		c = 255;	if(gam.neg)		c = 255-c;	return c;}