solar.c

#include <u.h>
#include <libc.h>
#include <bio.h>
#include <thread.h>
#include <draw.h>
#include <memdraw.h>
#include <mouse.h>
#include <keyboard.h>
#include <geometry.h>
#include "libgraphics/graphics.h"
#include "fns.h"

enum {
	K↑,
	K↓,
	K←,
	K→,
	Krise,
	Kfall,
	KR↑,
	KR↓,
	KR←,
	KR→,
	KR↺,
	KR↻,
	Kzoomin,
	Kzoomout,
	Khud,
	Ke
};

enum {
	Sfov,
	Scampos,
	Scambx, Scamby, Scambz,
	Sfps,
	Sframes,
	Se
};

enum {
	Cmdwinht = 50,
	Cmdmargin = 10,
	Cmdpadding = 3,

	Cmdlookat = 0,
	Cmdgoto,
};

typedef struct Planet Planet;
typedef struct Camcfg Camcfg;
typedef struct HReq HReq;
typedef struct Cmdbut Cmdbut;
typedef struct Cmdbox Cmdbox;
typedef struct Infobox Infobox;

struct Planet
{
	int id;		/* Horizons API ID */
	char *name;
	double scale;
	Entity *body;
	Material *mtl;
};

struct Camcfg
{
	Point3 p, lookat, up;
	double fov, clipn, clipf;
	int ptype;
};

struct HReq
{
	int pfd[2];
	int pid;	/* planet id */
	char *t0, *t1;	/* start and end times */
};

struct Cmdbut
{
	char *label;
	Rectangle r;
	void (*handler)(Cmdbut*);
};

struct Cmdbox
{
	Rectangle r;
	Cmdbut *cmds;
	ulong ncmds;
};

struct Infobox
{
	char *title;
	char **items;
	Image *image;
};

Rune keys[Ke] = {
 [K↑]		= Kup,
 [K↓]		= Kdown,
 [K←]		= Kleft,
 [K→]		= Kright,
 [Krise]	= Kpgup,
 [Kfall]	= Kpgdown,
 [KR↑]		= 'w',
 [KR↓]		= 's',
 [KR←]		= 'a',
 [KR→]		= 'd',
 [KR↺]		= 'q',
 [KR↻]		= 'e',
 [Kzoomin]	= 'z',
 [Kzoomout]	= 'x',
 [Khud]		= 'h',
};
char *skyboxpaths[] = {
	"cubemap/solar/left.pic",
	"cubemap/solar/right.pic",
	"cubemap/solar/bottom.pic",
	"cubemap/solar/top.pic",
	"cubemap/solar/front.pic",
	"cubemap/solar/back.pic",
};
Planet planets[] = {
	{ .id = 10,	.name = "Sol",		.scale = 695700 },
	{ .id = 1,	.name = "Mercury",	.scale = 2439.4 },
	{ .id = 2,	.name = "Venus",	.scale = 6051.8 },
	{ .id = 399,	.name = "Earth",	.scale = 6371.0084 },
	{ .id = 301,	.name = "Luna",		.scale = 1737.4 },
	{ .id = 4,	.name = "Mars",		.scale = 3389.50 },
	{ .id = 5,	.name = "Jupiter",	.scale = 69911 },
	{ .id = 6,	.name = "Saturn",	.scale = 58232 },
	{ .id = 7,	.name = "Uranus",	.scale = 25362 },
	{ .id = 8,	.name = "Neptune",	.scale = 24622 },
	{ .id = 9,	.name = "Pluto",	.scale = 1188.3 },
};
Planet *selplanet;
char stats[Se][256];
char datefmt[] = "YYYY-MM-DD";
Rectangle viewr, cmdr;
Cmdbox cmdbox;
Infobox *infobox;
Image *screenb;
Mousectl *mctl;
Keyboardctl *kctl;
Channel *drawc;
Mouse om;
int kdown;
Tm date;
char datestr[16];
Scene *scene;

Camera *camera;
Camcfg cameracfg = {
	0,0,0,1,
	0,0,1,0,
	0,1,0,0,
	80*DEG, 1, 1e12, PERSPECTIVE
};
Point3 center = {0,0,0,1};
double speed = 10;

static int museummode;
static int showskybox;
static int doprof;
static int showhud;

static void
refreshinfobox(Infobox *info)
{
	static Point pad = {2, 2};
	static Image *tbg, *tfg, *bg[2], *fg;
	Rectangle tr, cr, rr;	/* title, content and row rects */
	char **s;
	int i;

	assert(info != nil && info->image != nil && info->items != nil);

	if(tbg == nil){
		tbg = eallocimage(display, UR, RGBA32, 1, 0x2222447F);
		tfg = display->white;
		bg[0] = eallocimage(display, UR, RGBA32, 1, 0xEEEEEEEE);
		bg[1] = eallocimage(display, UR, RGBA32, 1, 0xAAAAAAAA);
		fg = display->black;
	}

	tr = info->image->r;
	tr.max.y = tr.min.y + 2*font->height;

	draw(info->image, tr, tbg, nil, ZP);
	string(info->image, addpt(tr.min, Pt(Dx(tr)/2 - stringwidth(font, info->title)/2,font->height/2)),
		tfg, ZP, font, info->title);

	cr = info->image->r;
	cr.min.y = tr.max.y;

	draw(info->image, cr, bg[1], nil, ZP);

	rr.min = addpt(cr.min, pad);
	rr.max = subpt(cr.max, pad);
	rr.max.y = rr.min.y + font->height+4;

	for(s = info->items, i = 0; *s != nil; s++, i ^= 1){
		draw(info->image, rr, bg[i], nil, ZP);
		string(info->image, rr.min, fg, ZP, font, *s);
		rr = rectaddpt(rr, Pt(0,Dy(rr)));
	}

	border(info->image, info->image->r, 1, tbg, ZP);
}

static Infobox *
mkplanetinfobox(Planet *p, Rectangle r)
{
	enum { ID, POS, RADIUS, NITEMS };
	Infobox *info;
	char **items, buf[256];
	int i;

	if(p == nil || badrect(r))
		return nil;

	items = emalloc((NITEMS + 1)*sizeof(*items));
	for(i = 0; i < NITEMS; i++){
		switch(i){
		case ID: snprint(buf, sizeof buf, "id: %d", p->id); break;
		case POS: snprint(buf, sizeof buf, "position (in km): %V", p->body->p); break;
		case RADIUS: snprint(buf, sizeof buf, "radius (in km): %g", p->scale); break;
		}
		items[i] = strdup(buf);
		if(items[i] == nil)
			sysfatal("strdup: %r");
	}
	items[i] = nil;

	info = emalloc(sizeof *info);
	memset(info, 0, sizeof *info);
	info->title = p->name;
	info->items = items;
	info->image = eallocimage(display, r, RGBA32, 0, DTransparent);
	refreshinfobox(info);
	return info;
}

static void
freeinfobox(Infobox *info)
{
	char **s;

	if(info == nil)
		return;

	freeimage(info->image);
	for(s = info->items; *s != nil; s++)
		free(*s);
	free(info);
}

static void
drawinfobox(Image *dst, Infobox *info)
{
	if(info == nil)
		return;

	draw(dst, rectaddpt(info->image->r, dst->r.min), info->image, nil, info->image->r.min);
}

static void
selectplanet(Planet *p)
{
	static Planet *oldp;
	struct { Point3 min, max; } aabb;
	Entity *e, *esel;
	Model *msel;
	Primitive l;
	int i, j;

	if(p == oldp)
		return;

	oldp = selplanet = p;
	esel = scene->getent(scene, "selection");
	if(esel != nil)
		scene->delent(scene, esel);
	lockdisplay(display);
	freeinfobox(infobox);
	infobox = nil;
	unlockdisplay(display);
	if(p == nil)
		return;

	e = p->body;
	msel = newmodel();
	esel = newentity("selection", msel);
	esel->RFrame3 = e->RFrame3;

	lockdisplay(display);
	infobox = mkplanetinfobox(p, Rpt(subpt(viewr.max, Pt(400,200)), viewr.max));
	unlockdisplay(display);

	memset(&aabb, 0, sizeof aabb);
	for(i = 0; i < e->mdl->nprims; i++)
		for(j = 0; j < e->mdl->prims[i].type+1; j++){
			aabb.min = minpt3(aabb.min, e->mdl->prims[i].v[j].p);
			aabb.max = maxpt3(aabb.max, e->mdl->prims[i].v[j].p);
		}
	aabb.min = mulpt3(aabb.min, p->scale*0.8);
	aabb.max = mulpt3(aabb.max, p->scale*0.8);
	aabb.min.w = aabb.max.w = 1;

	memset(&l, 0, sizeof l);
	l.type = PLine;
	l.v[0].c = l.v[1].c = Pt3(0.2666, 0.5333, 0.2666, 1);
	/* bottom */
	l.v[0].p = aabb.min; l.v[1].p = qrotate(aabb.min, Vec3(0,1,0), PI/2);
	msel->addprim(msel, l);
	for(i = 0; i < 3; i++){
		l.v[0].p = l.v[1].p; l.v[1].p = qrotate(l.v[1].p, Vec3(0,1,0), PI/2);
		msel->addprim(msel, l);
	}
	/* top */
	l.v[0].p = aabb.max; l.v[1].p = qrotate(aabb.max, Vec3(0,1,0), PI/2);
	msel->addprim(msel, l);
	for(i = 0; i < 3; i++){
		l.v[0].p = l.v[1].p; l.v[1].p = qrotate(l.v[1].p, Vec3(0,1,0), PI/2);
		msel->addprim(msel, l);
	}
	/* struts */
	l.v[0].p = aabb.min; l.v[1].p = qrotate(aabb.max, Vec3(0,1,0), PI);
	msel->addprim(msel, l);
	for(i = 0; i < 3; i++){
		l.v[0].p = qrotate(l.v[0].p, Vec3(0,1,0), PI/2);
		l.v[1].p = qrotate(l.v[1].p, Vec3(0,1,0), PI/2);
		msel->addprim(msel, l);
	}

	scene->addent(scene, esel);
}

static void
sailor(void *arg)
{
	char buf[128], pidstr[8];
	HReq *r;

	r = arg;

	close(r->pfd[0]);
	dup(r->pfd[1], 1);
	close(r->pfd[1]);

	getwd(buf, sizeof(buf)-1);
	snprint(buf+strlen(buf), sizeof(buf)-strlen(buf), "/tools/horizonget");
	snprint(pidstr, sizeof pidstr, "%d", r->pid);

	execl(buf, "horizonget", pidstr, r->t0, r->t1, nil);
	sysfatal("execl: %r");
}

static char *
getplanetstate(int id, Tm *t)
{
	Biobuf *bin;
	char *line, *lastline, t0[16], t1[16];
	HReq r;

	lastline = nil;
	r.pid = id;
	r.t0 = t0;
	r.t1 = t1;

	snprint(t1, sizeof t1, "%τ", tmfmt(t, datefmt));
	t->mday--;
	snprint(t0, sizeof t0, "%τ", tmfmt(t, datefmt));
	t->mday++;

	if(pipe(r.pfd) < 0)
		sysfatal("pipe: %r");
	switch(fork()){
	case -1:
		sysfatal("fork: %r");
	case 0:
		sailor(&r);
	default:
		close(r.pfd[1]);
		bin = Bfdopen(r.pfd[0], OREAD);
		if(bin == nil)
			sysfatal("Bfdopen: %r");
		while((line = Brdline(bin, '\n')) != nil){
			line[Blinelen(bin)-1] = 0;
			lastline = line;
		}
		if(lastline != nil)
			lastline = strdup(lastline);
		Bterm(bin);
		close(r.pfd[0]);
	}

	return lastline;
}

void
updateplanets(void)
{
	char *s, *p;
	int i;

	fprint(2, "loading planet states...\n");
	for(i = 1; i < nelem(planets); i++){
		s = getplanetstate(planets[i].id, &date);
		if(s == nil){
			fprint(2, "couldn't load planet: %s", planets[i].name);
			continue;
		}
		p = strchr(s, '=');
		planets[i].body->p.x = strtod(++p, nil);
		p = strchr(p, '=');
		planets[i].body->p.y = strtod(++p, nil);
		p = strchr(p, '=');
		planets[i].body->p.z = strtod(++p, nil);
		planets[i].body->p.w = 1;
		free(s);
		fprint(2, "%s ready\n", planets[i].name);
	}
}

static Planet *
getplanet(char *name)
{
	int i;

	for(i = 0; i < nelem(planets); i++)
		if(strcmp(planets[i].name, name) == 0)
			return &planets[i];
	return nil;
}

static void
gotoplanet(Planet *p)
{
	movecamera(camera, addpt3(p->body->p, Vec3(0,0,1.5*p->scale)));
	aimcamera(camera, p->body->p);
}

Point3
identvshader(Shaderparams *sp)
{
	Planet *p;
	Point3 pos;

	p = getplanet(sp->su->entity->name);

	if(p != nil){
		Matrix3 S = {
			p->scale, 0, 0, 0,
			0, p->scale, 0, 0,
			0, 0, p->scale, 0,
			0, 0, 0, 1,
		};
		pos = xform3(sp->v->p, S);
		sp->v->mtl = p->mtl;
		sp->v->c = p->mtl->diffuse;
	}else
		pos = sp->v->p;

	return world2clip(sp->su->camera, model2world(sp->su->entity, pos));
}

Color
identshader(Shaderparams *sp)
{
	if(sp->v->mtl != nil && sp->v->mtl->diffusemap != nil && sp->v->uv.w != 0)
		return sampletexture(sp->v->mtl->diffusemap, sp->v->uv, neartexsampler);
	return sp->v->c;
}

Shadertab shader = { "ident", identvshader, identshader };

void
zoomin(void)
{
	camera->fov = fclamp(camera->fov - 1*DEG, 1*DEG, 180*DEG);
	reloadcamera(camera);
}

void
zoomout(void)
{
	camera->fov = fclamp(camera->fov + 1*DEG, 1*DEG, 180*DEG);
	reloadcamera(camera);
}

void
drawstats(void)
{
	int i;

	snprint(stats[Sfov], sizeof(stats[Sfov]), "FOV %g°", camera->fov/DEG);
	snprint(stats[Scampos], sizeof(stats[Scampos]), "%V", camera->p);
	snprint(stats[Scambx], sizeof(stats[Scambx]), "bx %V", camera->bx);
	snprint(stats[Scamby], sizeof(stats[Scamby]), "by %V", camera->by);
	snprint(stats[Scambz], sizeof(stats[Scambz]), "bz %V", camera->bz);
	snprint(stats[Sfps], sizeof(stats[Sfps]), "FPS %.0f/%.0f/%.0f/%.0f",
		!camera->stats.max? 0: 1e9/camera->stats.max,
		!camera->stats.avg? 0: 1e9/camera->stats.avg,
		!camera->stats.min? 0: 1e9/camera->stats.min,
		!camera->stats.v? 0: 1e9/camera->stats.v);
	snprint(stats[Sframes], sizeof(stats[Sframes]), "frame %llud", camera->stats.nframes);
	for(i = 0; i < Se; i++)
		stringbg(screen, addpt(screen->r.min, Pt(10,10 + i*font->height)), display->black, ZP, font, stats[i], display->white, ZP);
}

void
redraw(void)
{
	int i;

	lockdisplay(display);
	draw(screen, rectaddpt(viewr, screen->r.min), screenb, nil, ZP);
	draw(screen, rectaddpt(cmdbox.r, screen->r.min), display->white, nil, ZP);
	drawinfobox(screen, infobox);
	for(i = 0; i < cmdbox.ncmds; i++){
		border(screen, rectaddpt(cmdbox.cmds[i].r, screen->r.min), 1, display->black, ZP);
		string(screen, addpt(screen->r.min, addpt(cmdbox.cmds[i].r.min, Pt(Cmdpadding,Cmdpadding))), display->black, ZP, font, cmdbox.cmds[i].label);
	}
	if(showhud)
		drawstats();
	flushimage(display, 1);
	unlockdisplay(display);
}

void
renderproc(void *)
{
	static Image *bg;
	uvlong t0, Δt;

	threadsetname("renderproc");

	bg = eallocimage(display, UR, XRGB32, 1, DBlack);

	t0 = nsec();
	for(;;){
		shootcamera(camera, &shader);
		Δt = nsec() - t0;
		if(Δt > HZ2MS(60)*1000000ULL){
			lockdisplay(display);
			draw(screenb, screenb->r, bg, nil, ZP);
			camera->view->draw(camera->view, screenb, nil);
			unlockdisplay(display);
			nbsend(drawc, nil);
			t0 += Δt;
		}
	}
}

void
drawproc(void *)
{
	threadsetname("drawproc");

	for(;;){
		recv(drawc, nil);
		redraw();
	}
}

static char *
genplanetmenu(int idx)
{
	if(idx < nelem(planets))
		return planets[idx].name;
	return nil;
}

void
lookat_cmd(Cmdbut *)
{
	static Menu menu = { .gen = genplanetmenu };
	Planet *p;
	int idx;

	lockdisplay(display);
	idx = menuhit(1, mctl, &menu, _screen);
	if(idx >= 0){
		p = &planets[idx];
		aimcamera(camera, p->body->p);
	}
	unlockdisplay(display);
	nbsend(drawc, nil);
}

void
goto_cmd(Cmdbut *)
{
	static Menu menu = { .gen = genplanetmenu };
	int idx;

	lockdisplay(display);
	idx = menuhit(1, mctl, &menu, _screen);
	if(idx >= 0)
		gotoplanet(&planets[idx]);
	unlockdisplay(display);
	nbsend(drawc, nil);
}

void
date_cmd(Cmdbut *)
{
	Tm t;
	char buf[16];

	if(museummode)
		return;

	memmove(buf, datestr, sizeof buf);
	lockdisplay(display);
	if(enter("new date", buf, sizeof buf, mctl, kctl, nil) <= 0)
		goto nodate;
	if(tmparse(&t, datefmt, buf, nil, nil) == nil)
		goto nodate;
	date = t;
	snprint(datestr, sizeof datestr, "%τ", tmfmt(&date, datefmt));
	updateplanets();
nodate:
	unlockdisplay(display);
	nbsend(drawc, nil);
}

Cmdbut cmds[] = {
	{ .label = "look at", .handler = lookat_cmd },
	{ .label = "go to", .handler = goto_cmd },
	{ .label = datestr, .handler = date_cmd },
};

void
lmb(void)
{
	Point3 p0;
	Point mp;
	Planet *p;
	Cmdbut *cmd;
	double lastz, z;
	int i;

	if((om.buttons ^ mctl->buttons) == 0)
		return;

	mp = subpt(mctl->xy, screen->r.min);
	if(ptinrect(mp, viewr)){
		p0 = viewport2world(camera, Pt3(mp.x,mp.y,1,1));
		p = nil;
		lastz = Inf(1);

		for(i = 0; i < nelem(planets); i++)
			if(lineXsphere(nil, camera->p, p0, planets[i].body->p, planets[i].scale, 1) > 0){
				z = vec3len(subpt3(planets[i].body->p, camera->p));
				/* select the closest one */
				if(z < lastz){
					lastz = z;
					p = &planets[i];
				}
			}
		selectplanet(p);
		return;
	}

	cmd = nil;
	for(i = 0; i < cmdbox.ncmds; i++)
		if(ptinrect(mp, cmdbox.cmds[i].r))
			cmd = &cmdbox.cmds[i];

	if(cmd == nil)
		return;
	cmd->handler(cmd);
}

void
mmb(void)
{
	enum {
		CHGSPEED,
		QUIT,
	};
	static char *items[] = {
	 [CHGSPEED]	"change speed",
	 [QUIT]	"quit",
		nil,
	};
	static Menu menu = { .item = items };
	char buf[128];

	if((om.buttons ^ mctl->buttons) == 0)
		return;

	lockdisplay(display);
	switch(menuhit(2, mctl, &menu, _screen)){
	case CHGSPEED:
		snprint(buf, sizeof buf, "%g", speed);
		if(enter("speed (km)", buf, sizeof buf, mctl, kctl, nil) > 0)
			speed = fabs(strtod(buf, nil));
		break;
	case QUIT:
		threadexitsall(nil);
	}
	unlockdisplay(display);
	nbsend(drawc, nil);
}

void
mouse(void)
{
	if((mctl->buttons & 1) != 0)
		lmb();
	if((mctl->buttons & 2) != 0)
		mmb();
	if((mctl->buttons & 8) != 0)
		zoomin();
	if((mctl->buttons & 16) != 0)
		zoomout();

	om = mctl->Mouse;
}

void
kbdproc(void *)
{
	Rune r, *a;
	char buf[128], *s;
	int fd, n;

	threadsetname("kbdproc");

	if((fd = open("/dev/kbd", OREAD)) < 0)
		sysfatal("kbdproc: %r");
	memset(buf, 0, sizeof buf);

	for(;;){
		if(buf[0] != 0){
			n = strlen(buf)+1;
			memmove(buf, buf+n, sizeof(buf)-n);
		}
		if(buf[0] == 0){
			if((n = read(fd, buf, sizeof(buf)-1)) <= 0)
				break;
			buf[n-1] = 0;
			buf[n] = 0;
		}
		if(buf[0] == 'c'){
			chartorune(&r, buf+1);
			if(r == Kdel){
				close(fd);
				threadexitsall(nil);
			}else
				nbsend(kctl->c, &r);
		}
		if(buf[0] != 'k' && buf[0] != 'K')
			continue;
		s = buf+1;
		kdown = 0;
		while(*s){
			s += chartorune(&r, s);
			for(a = keys; a < keys+Ke; a++)
				if(r == *a){
					kdown |= 1 << a-keys;
					break;
				}
		}
	}
}

void
keyproc(void *c)
{
	threadsetname("keyproc");

	for(;;){
		nbsend(c, nil);
		sleep(HZ2MS(100));	/* key poll rate */
	}
}

void
handlekeys(void)
{
	static int okdown;

	if(kdown & 1<<K↑)
		movecamera(camera, mulpt3(camera->bz, -speed));
	if(kdown & 1<<K↓)
		movecamera(camera, mulpt3(camera->bz, speed));
	if(kdown & 1<<K←)
		movecamera(camera, mulpt3(camera->bx, -speed));
	if(kdown & 1<<K→)
		movecamera(camera, mulpt3(camera->bx, speed));
	if(kdown & 1<<Krise)
		movecamera(camera, mulpt3(camera->by, speed));
	if(kdown & 1<<Kfall)
		movecamera(camera, mulpt3(camera->by, -speed));
	if(kdown & 1<<KR↑)
		rotatecamera(camera, camera->bx, 1*DEG);
	if(kdown & 1<<KR↓)
		rotatecamera(camera, camera->bx, -1*DEG);
	if(kdown & 1<<KR←)
		rotatecamera(camera, camera->by, 1*DEG);
	if(kdown & 1<<KR→)
		rotatecamera(camera, camera->by, -1*DEG);
	if(kdown & 1<<KR↺)
		rotatecamera(camera, camera->bz, 1*DEG);
	if(kdown & 1<<KR↻)
		rotatecamera(camera, camera->bz, -1*DEG);
	if(kdown & 1<<Kzoomin)
		zoomin();
	if(kdown & 1<<Kzoomout)
		zoomout();

	if((okdown & 1<<Khud) == 0 && (kdown & 1<<Khud) != 0)
		showhud ^= 1;

	okdown = kdown;
}

void
resize(void)
{
	lockdisplay(display);
	if(getwindow(display, Refnone) < 0)
		fprint(2, "can't reattach to window\n");
	unlockdisplay(display);
	nbsend(drawc, nil);
}

static void
confproc(void)
{
	char buf[64];
	int fd;

	snprint(buf, sizeof buf, "/proc/%d/ctl", getpid());
	fd = open(buf, OWRITE);
	if(fd < 0)
		sysfatal("open: %r");

	if(doprof)
		fprint(fd, "profile\n");
	close(fd);
}

void
usage(void)
{
	fprint(2, "usage: %s [-ms]\n", argv0);
	exits("usage");
}

void
threadmain(int argc, char *argv[])
{
	Renderer *rctl;
	Channel *keyc;
	Entity *subject;
	Model *model;
	Point lblsiz;
	int fd, i, j;

	tmfmtinstall();
	GEOMfmtinstall();
	ARGBEGIN{
	case 'm': museummode++; break;
	case 's': showskybox++; break;
	case 'p': doprof++; break;
	default: usage();
	}ARGEND;
	if(argc != 0)
		usage();

	confproc();

	fd = open("mdl/planet.mdl", OREAD);
	if(fd < 0)
		sysfatal("open: %r");
	model = readmodel(fd);
	if(model == nil)
		sysfatal("readmodel: %r");
	close(fd);
	/*
	 * normalize the vertices so that we can scale
	 * each planet based on its radius
	 */
	for(i = 0; i < model->nprims; i++)
		for(j = 0; j < model->prims[i].type+1; j++){
			model->prims[i].v[j].p = normvec3(model->prims[i].v[j].p);
			model->prims[i].v[j].p.w = 1;
		}
	scene = newscene(nil);
	for(i = 0; i < nelem(planets); i++){
		subject = newentity(planets[i].name, model);
		scene->addent(scene, subject);
		planets[i].body = subject;
		for(j = 0; j < model->nmaterials; j++)
			if(strcmp(planets[i].name, model->materials[j].name) == 0)
				planets[i].mtl = &model->materials[j];
		if(i == 0){
			subject->p = Pt3(0,0,0,1);
			continue;
		}else if(museummode)
			subject->p.x = planets[i-1].body->p.x + 1.5*planets[i-1].scale + planets[i].scale;
	}
	tmnow(&date, nil);
	snprint(datestr, sizeof datestr, "%τ", tmfmt(&date, datefmt));
	if(!museummode)
		updateplanets();
	if(showskybox)
		scene->skybox = readcubemap(skyboxpaths);

	if(memimageinit() != 0)
		sysfatal("memimageinit: %r");
	if((rctl = initgraphics()) == nil)
		sysfatal("initgraphics: %r");
	if(initdraw(nil, nil, "solar") < 0)
		sysfatal("initdraw: %r");
	if((mctl = initmouse(nil, screen)) == nil)
		sysfatal("initmouse: %r");

	rctl->doprof = doprof;

	viewr = rectsubpt(Rpt(screen->r.min, subpt(screen->r.max, Pt(0,Cmdwinht))), screen->r.min);
	cmdbox.r = Rect(viewr.min.x, viewr.max.y, Dx(viewr), Dy(screen->r));
	cmdbox.cmds = cmds;
	cmdbox.ncmds = nelem(cmds);
	for(i = 0; i < nelem(cmds); i++){
		lblsiz = stringsize(font, cmds[i].label);
		cmds[i].r = Rect(0,0,Cmdpadding+lblsiz.x+Cmdpadding,Cmdpadding+lblsiz.y+Cmdpadding);
		if(i == 0)
			cmds[i].r = rectaddpt(cmds[i].r, addpt(cmdbox.r.min, Pt(Cmdmargin,Cmdmargin)));
		else
			cmds[i].r = rectaddpt(cmds[i].r, Pt(cmds[i-1].r.max.x+Cmdmargin,cmds[i-1].r.min.y));
	}

	screenb = eallocimage(display, viewr, XRGB32, 0, DNofill);
	camera = Cam(screenb->r, rctl, cameracfg.ptype, cameracfg.fov, cameracfg.clipn, cameracfg.clipf);
	placecamera(camera, scene, cameracfg.p, cameracfg.lookat, cameracfg.up);
	camera->cullmode = CullBack;
	gotoplanet(getplanet("Sol"));

	kctl = emalloc(sizeof *kctl);
	kctl->c = chancreate(sizeof(Rune), 16);
	keyc = chancreate(sizeof(void*), 1);
	drawc = chancreate(sizeof(void*), 1);
	display->locking = 1;
	unlockdisplay(display);

	proccreate(kbdproc, nil, mainstacksize);
	proccreate(keyproc, keyc, mainstacksize);
	proccreate(renderproc, nil, mainstacksize);
	proccreate(drawproc, nil, mainstacksize);

	for(;;){
		enum {MOUSE, RESIZE, KEY};
		Alt a[] = {
			{mctl->c, &mctl->Mouse, CHANRCV},
			{mctl->resizec, nil, CHANRCV},
			{keyc, nil, CHANRCV},
			{nil, nil, CHANEND}
		};
		switch(alt(a)){
		case MOUSE: mouse(); break;
		case RESIZE: resize(); break;
		case KEY: handlekeys(); break;
		}
	}
}