You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
81 lines
1.6 KiB
C
81 lines
1.6 KiB
C
#include "SDL.h"
|
|
#include <complex.h>
|
|
#include <stdbool.h>
|
|
|
|
/* Compute the out-coloring based on the iteration counter. */
|
|
Uint32 outcolor(int it) {
|
|
return 0x00010001 * it;
|
|
}
|
|
|
|
/* Compute the in-coloring. */
|
|
Uint32 incolor() {
|
|
return 0x00000000; /* black */
|
|
}
|
|
|
|
void drawmandelbrot(SDL_Surface *surface) {
|
|
if (SDL_MUSTLOCK(surface)) {
|
|
SDL_LockSurface(surface);
|
|
}
|
|
Uint32 *pixels = (Uint32 *) surface->pixels;
|
|
|
|
for (int i = 0; i < surface->w * surface->h; i++) {
|
|
int y = i / surface->w;
|
|
int x = i % surface->w;
|
|
|
|
float complex c = ((3.0 * x / surface->w) - 2.0)
|
|
+ I * ((2.0 * y / surface->h) - 1.0);
|
|
|
|
bool diverges = false;
|
|
float complex z = 0;
|
|
int it;
|
|
for (it = 1; it < 256; it++) {
|
|
z = cpowf(z, 2) + c;
|
|
if (cabs(z) > 2) {
|
|
diverges = true;
|
|
break;
|
|
}
|
|
}
|
|
|
|
Uint32 color;
|
|
if (diverges) {
|
|
color = outcolor(it);
|
|
} else {
|
|
color = incolor();
|
|
}
|
|
|
|
pixels[i] = color;
|
|
|
|
/* Update the screen every 10 lines. */
|
|
if (y % 10 == 0) {
|
|
SDL_Flip(surface);
|
|
}
|
|
}
|
|
|
|
if (SDL_MUSTLOCK(surface)) {
|
|
SDL_UnlockSurface(surface);
|
|
}
|
|
SDL_Flip(surface);
|
|
}
|
|
|
|
int main(int argc, char *argv[]) {
|
|
/* Set up SDL. */
|
|
if (SDL_Init(SDL_INIT_EVERYTHING) < 0) {
|
|
fprintf(stderr, "Error: Could not initialize SDL: %s.\n", SDL_GetError());
|
|
exit(1);
|
|
}
|
|
SDL_Surface *screen = SDL_SetVideoMode(320, 240, 32, SDL_SWSURFACE);
|
|
|
|
/* Do the mandelbrot. */
|
|
drawmandelbrot(screen);
|
|
|
|
/* Save BMP */
|
|
char *file = "mandelbrot.bmp";
|
|
if (SDL_SaveBMP(screen, file) != 0) {
|
|
fprintf(stderr, "Could not write %s!\n", file);
|
|
}
|
|
|
|
/* Quit. */
|
|
SDL_Delay(20000);
|
|
SDL_Quit();
|
|
}
|