*b = pow(iterations, color_mandelbrot_parms[12]);
}
+struct
+{
+ int n;
+ double darkness;
+ double power;
+ double density;
+}
+color_starfield_parms;
+void color_starfield(double x, double y, double z, double *r, double *g, double *b)
+{
+ static struct
+ {
+ double x, y, z, e;
+ double R, G, B, A;
+ } *starfield = NULL;
+ int i;
+ double f;
+
+ if(!starfield)
+ {
+ fprintf(stderr, "Initializing starfield...\n");
+ starfield = malloc(sizeof(*starfield) * color_starfield_parms.n);
+ for(i = 0; i < color_starfield_parms.n; ++i)
+ {
+ double r;
+ do
+ {
+ starfield[i].x = rnd() * 2 - 1;
+ starfield[i].y = rnd() * 2 - 1;
+ starfield[i].z = rnd() * 2 - 1;
+ r = starfield[i].x * starfield[i].x
+ + starfield[i].y * starfield[i].y
+ + starfield[i].z * starfield[i].z;
+ }
+ while(r > 1);
+ r = sqrt(r);
+ starfield[i].x /= r;
+ starfield[i].y /= r;
+ starfield[i].z /= r;
+
+ starfield[i].e = color_starfield_parms.density * pow(rnd(), -color_starfield_parms.power);
+
+ starfield[i].R = rnd();
+ starfield[i].G = rnd();
+ starfield[i].B = rnd();
+ f = starfield[i].R * 0.299 + starfield[i].G * 0.587 + starfield[i].B * 0.114;
+ starfield[i].R /= f;
+ starfield[i].G /= f;
+ starfield[i].B /= f;
+ starfield[i].A = rnd();
+ }
+ fprintf(stderr, "Done.\n");
+ }
+
+ *r = *g = *b = 0;
+ for(i = 0; i < color_starfield_parms.n; ++i)
+ {
+ double dot = x * starfield[i].x + y * starfield[i].y + z * starfield[i].z;
+ if(dot <= 0)
+ continue;
+ double f = pow(dot, starfield[i].e) * starfield[i].A;
+ *r += starfield[i].R * f;
+ *g += starfield[i].G * f;
+ *b += starfield[i].B * f;
+ }
+ // make fit in 0..1
+ *r = *r / (color_starfield_parms.darkness + *r);
+ *g = *g / (color_starfield_parms.darkness + *g);
+ *b = *b / (color_starfield_parms.darkness + *b);
+}
+
void map_back(double x_in, double y_in, double *x_out, double *y_out, double *z_out)
{
*x_out = 2 * x_in - 1;
{
int x, y;
uint8_t tga[18];
+ int percent, p;
FILE *file = fopen(fn, "wb");
if(!file)
tga[16] = 24; // pixel size
fwrite(&tga, sizeof(tga), 1, file);
+ percent = 0;
for(y = height-1; y >= 0; --y)
+ {
for(x = 0; x < width; ++x)
{
uint8_t rgb[3];
rgb[0] = floor(rnd() + bb * 255);
fwrite(rgb, sizeof(rgb), 1, file);
}
+ p = (100 * (height - y)) / height;
+ if(p != percent)
+ {
+ percent = p;
+ fprintf(stderr, "%d%%\r", percent);
+ }
+ }
+ fprintf(stderr, "\n");
fclose(file);
}
void map_all(const char *fn, colorfunc_t f, int width, int height)
{
char buf[1024];
+ fprintf(stderr, "%s_bk.tga\n", fn);
snprintf(buf, sizeof(buf), "%s_bk.tga", fn); buf[sizeof(buf) - 1] = 0; writepic(f, map_back, buf, width, height);
+ fprintf(stderr, "%s_ft.tga\n", fn);
snprintf(buf, sizeof(buf), "%s_ft.tga", fn); buf[sizeof(buf) - 1] = 0; writepic(f, map_front, buf, width, height);
+ fprintf(stderr, "%s_rt.tga\n", fn);
snprintf(buf, sizeof(buf), "%s_rt.tga", fn); buf[sizeof(buf) - 1] = 0; writepic(f, map_right, buf, width, height);
+ fprintf(stderr, "%s_lf.tga\n", fn);
snprintf(buf, sizeof(buf), "%s_lf.tga", fn); buf[sizeof(buf) - 1] = 0; writepic(f, map_left, buf, width, height);
+ fprintf(stderr, "%s_up.tga\n", fn);
snprintf(buf, sizeof(buf), "%s_up.tga", fn); buf[sizeof(buf) - 1] = 0; writepic(f, map_up, buf, width, height);
+ fprintf(stderr, "%s_dn.tga\n", fn);
snprintf(buf, sizeof(buf), "%s_dn.tga", fn); buf[sizeof(buf) - 1] = 0; writepic(f, map_down, buf, width, height);
}
color_mandelbrot_parms[11] = argc<=15 ? 0.5 : atof(argv[15]);
color_mandelbrot_parms[12] = argc<=16 ? 0.2 : atof(argv[16]);
}
+ else if(!strcmp(argv[3], "starfield"))
+ {
+ f = color_starfield;
+ color_starfield_parms.n = argc<= 4 ? 8192 : atoi(argv[4]);
+ color_starfield_parms.darkness = argc<= 5 ? 0.4 : atof(argv[5]);
+ color_starfield_parms.power = argc<= 6 ? 2.5 : atof(argv[6]);
+ color_starfield_parms.density = argc<= 7 ? 60000 : atof(argv[7]);
+ }
else
{
f = color_test;
}
- map_all(argv[1], color_mandelbrot, res, res);
+ map_all(argv[1], f, res, res);
return 0;
}