#include <unistd.h>
#include <stdio.h>
#include <time.h>
+#include <getopt.h>
#include "random.h"
#include "planet.h"
#define MAX_FPS 60
-static void fade_buf(SDL_Surface *screen, int amount)
+static void fade_buf(SDL_Surface *screen, double fade)
{
int i;
unsigned int *buf = screen->pixels;
unsigned char *b;
+ unsigned char amount;
+ static double threshold;
+ threshold += fade;
+
+ if (threshold > 1) {
+ amount = threshold;
+ threshold -= amount;
+ } else {
+ return;
+ }
- for (i = 0; i < screen->pitch * screen->h / sizeof(*buf); i ++) {
+ for (i = 0; i < screen->pitch * screen->h / sizeof(*buf); i++) {
if (!buf[i])
continue;
b = (unsigned char *)&buf[i];
{
struct timespec tp;
unsigned long ret;
- static unsigned long start = 0;
+ static unsigned long start;
clock_gettime(CLOCK_MONOTONIC, &tp);
case SDLK_a:
time_scale_rate = 1 / 1.5;
break;
+ case SDLK_1:
+ status->time_scale = 1;
+ break;
+ case SDLK_2:
+ status->time_scale = 10;
+ break;
+ case SDLK_3:
+ status->time_scale = 20;
+ break;
+ case SDLK_4:
+ status->time_scale = 30;
+ break;
+ case SDLK_0:
+ status->time_scale = 0;
+ break;
case SDLK_t:
status->tracers_enabled =
!status->tracers_enabled;
}
break;
case SDL_VIDEORESIZE:
- status->screen = SDL_SetVideoMode(event.resize.w,
- event.resize.h,
- 32,
- status->screen->flags);
+ status->screen =
+ SDL_SetVideoMode(event.resize.w,
+ event.resize.h,
+ 32,
+ status->screen->flags);
break;
case SDL_QUIT:
goto quit;
double range)
{
struct sim_status status;
- struct planet *planet, *pl1, *pl2;
+ struct planet *planet, *pl1, *pl2, *planet_root;
struct camera camera;
+ struct planet_search_iterator itr;
+ struct vector vect;
+
int planets;
int framecount = 0, last_fps_time = 0;
int last_framecount = 0, last_step_count = 0;
unsigned long step_count = 0;
double last_fps = 0, last_sps = 0, fade_amount;
double step_time = 0, true_time = 0;
+ int visible_planets;
init_camera(&camera);
status.screen = screen;
status.tracers_enabled = 0;
+ itr.screen = screen;
+ itr.cam = &camera;
+ itr.qt_iterator.direction = planet_search_rectangular;
+ itr.qt_iterator.callback = planet_draw_iterator;
+
planet = malloc(sizeof(*planet));
init_planet(planet);
create_planets(planet, num_of_planets, total_mass, range);
while (1) {
planets = 0;
- list_for_each_entry(pl1, &planet->list, list) {
- pl2 = list_to_planet(pl1->list.next);
- list_for_each_entry_from(pl2, &planet->list, list) {
- struct planet *ptmp;
- if (!gravitize_planets(pl1, pl2,
- step_time))
- continue;
-
- ptmp = list_to_planet(pl2->list.prev);
- merge_planets(pl1, pl2);
- pl2 = ptmp;
+ if (status.time_scale > 0) {
+ list_for_each_entry(pl1, &planet->list, list) {
+ pl2 = list_to_planet(pl1->list.next);
+ list_for_each_entry_from(pl2, &planet->list,
+ list) {
+ struct planet *ptmp;
+ if (!gravitize_planets(pl1, pl2,
+ step_time))
+ continue;
+
+ ptmp = list_to_planet(pl2->list.prev);
+ merge_planets(pl1, pl2);
+ pl2 = ptmp;
+ }
+
+ planet_root = move_planet(pl1, step_time);
+ planets++;
}
-
- move_planet(pl1, step_time);
}
-
move_camera(&camera, true_time);
if (poll_events(&status, true_time))
step_count++;
/*
- * Do not draw to the screen more often than MAS_FPS
+ * Do not draw to the screen more often than MAX_FPS
* per second
*/
if (last_frame_tick + (1000 * 1000) / MAX_FPS > ticks)
clear_buf(screen);
}
- list_for_each_entry(pl1, &planet->list, list) {
- draw_planet(screen, pl1, &camera);
- planets++;
- }
+ itr.limit[0] = camera.pos;
+ vect.x = screen->w / 2;
+ vect.y = screen->h / 2;
+ vector_scale(&vect, 1 / camera.zoom, &vect);
+ vector_add(&itr.limit[0], &vect, &itr.limit[1]);
+ vector_sub(&itr.limit[0], &vect, &itr.limit[0]);
+
+ itr.qt_iterator.head = &planet_root->tree;
+
+ visible_planets = walk_quadtree(&itr.qt_iterator);
SDL_UnlockSurface(screen);
last_fps_time = ticks;
}
-
- printf(" \rFrames/s: %.2f, steps/s: %.2f, planets: %d"
- ", scale %.2f, zoom %.2f, step %ld",
+ printf(" \rfps: %.2f, steps/s: %.2f, planets: %d"
+ ", scale %.2f, zoom %.2f, step %ld, visible %d,"
+ " depth %ld, c:%ld",
last_fps, last_sps, planets, status.time_scale,
- camera.zoom, step_count);
+ camera.zoom, step_count, visible_planets,
+ planet_root->tree.depth, planet_root->tree.children);
fflush(stdout);
int planets = 100, xres = 800, yres = 600;
double total_mass = 50000;
double range = 500;
+ int optind = 0, c;
+ static struct option long_options[] = {
+ { .val = 'x', .name = "xres", .has_arg = 1, },
+ { .val = 'y', .name = "yres", .has_arg = 1, },
+ { .val = 'p', .name = "planets", .has_arg = 1 },
+ { .val = 'm', .name = "total-mass", .has_arg = 1 },
+ { .val = 'r', .name = "range", .has_arg = 1 },
+ {},
+ };
+ char short_options[] = "x:y:p:m:r:";
if (SDL_Init(SDL_INIT_VIDEO) != 0) {
fprintf(stderr, "Unable to initialize SDL: %s\n",
}
atexit(SDL_Quit);
- if (argc >= 2)
- planets = atoi(argv[1]);
+ while (1) {
+ c = getopt_long(argc, argv, short_options, long_options,
+ &optind);
+
+ if (c == -1)
+ break;
- if (argc >= 3)
- total_mass = atof(argv[2]);
+ printf("%c: %s\n", c, optarg);
- if (argc >= 4)
- range = atof(argv[3]);
+ switch (c) {
+ case 'x':
+ xres = atoi(optarg);
+ break;
+
+ case 'y':
+ yres = atoi(optarg);
+ break;
+
+ case 'p':
+ planets = atoi(optarg);
+ break;
+
+ case 'r':
+ range = atoi(optarg);
+ break;
+ }
+ }
screen = SDL_SetVideoMode(xres, yres, 32, flags);
if (screen == NULL) {