Skip to content

Instantly share code, notes, and snippets.

@DavidBuchanan314
Last active June 28, 2023 15:12
Show Gist options
  • Star 106 You must be signed in to star a gist
  • Fork 8 You must be signed in to fork a gist
  • Save DavidBuchanan314/b9230fe7d335a1caf90483dbb00a5375 to your computer and use it in GitHub Desktop.
Save DavidBuchanan314/b9230fe7d335a1caf90483dbb00a5375 to your computer and use it in GitHub Desktop.
Compile-time mandelbrot in pure C. Outputs a PGM image file to stdout. Output can be seen at https://twitter.com/David3141593/status/1062468528115200001
#include <stdio.h>
#define SQ(x) (x)*(x)
#define M0(x,y) SQ(x)+SQ(y)<4?0:0xe0
#define M1(x,y,x0,y0) (SQ(x)+SQ(y)<4)?M0(SQ(x)-SQ(y)+(x0),2*(x)*(y)+(y0)):0xc0
#define M2(x,y,x0,y0) (SQ(x)+SQ(y)<4)?M1(SQ(x)-SQ(y)+(x0),2*(x)*(y)+(y0),x0,y0):0xa0
#define M3(x,y,x0,y0) (SQ(x)+SQ(y)<4)?M2(SQ(x)-SQ(y)+(x0),2*(x)*(y)+(y0),x0,y0):0x80
#define M4(x,y,x0,y0) (SQ(x)+SQ(y)<4)?M3(SQ(x)-SQ(y)+(x0),2*(x)*(y)+(y0),x0,y0):0x60
#define M5(x,y,x0,y0) (SQ(x)+SQ(y)<4)?M4(SQ(x)-SQ(y)+(x0),2*(x)*(y)+(y0),x0,y0):0x40
#define M6(x,y,x0,y0) (SQ(x)+SQ(y)<4)?M5(SQ(x)-SQ(y)+(x0),2*(x)*(y)+(y0),x0,y0):0x20
#define XF(x) (x)/20-2.2
#define YF(y) (y)/20-1.6
#define C1(x,y) M4(XF(x), YF(y), XF(x), YF(y)),
// ^- change this to 5 or 6 to increase detail
#define C2(x,y) C1(x,y) C1(x+1,y)
#define C4(x,y) C2(x,y) C2(x+2,y)
#define C8(x,y) C4(x,y) C4(x+4,y)
#define C16(x,y) C8(x,y) C8(x+8,y)
#define C32(x,y) C16(x,y) C16(x+16,y)
#define C64(x,y) C32(x,y) C32(x+32,y)
#define R2(y) C64(0.0,y) C64(0.0,y+1)
#define R4(y) R2(y) R2(y+2)
#define R8(y) R4(y) R4(y+4)
#define R16(y) R8(y) R8(y+8)
#define R32(y) R16(y) R16(y+16)
#define R64(y) R32(y) R32(y+32)
static const unsigned char pixels[] = {
R64(0.0)
};
int main() {
fputs("P5\n64 64\n255\n", stdout);
fwrite(pixels, sizeof(pixels), 1, stdout);
return 0;
}
@klotzambein
Copy link

Very nice, killed my laptop though.

@Nodraak
Copy link

Nodraak commented Nov 14, 2018

That was close. The compilation finished after 60 sec (cpu at 3.2 GHz) and 7 GB of ram.

@george-chakhidze
Copy link

george-chakhidze commented Nov 14, 2018

MSVC 2017 x64 (v19.15.26730): Error C1060: compiler is out of heap space.

@PhilipTrauner
Copy link

PhilipTrauner commented Nov 14, 2018

26GBs of RAM for the M6 variant 😅

test

@oherrala
Copy link

👍 Nice work! However, my definition of "pure C" doesn't include // comments. Those are something from the land of C++. My "pure C" include only /* comments */. But you can substitute your own definition of "pure C" in your own code any time. 😸

@ror3d
Copy link

ror3d commented Nov 14, 2018

@oherrala it is if you use C99. Maybe what you mean by "pure C" is really "old C".

@oherrala
Copy link

@rorlork Nothing wrong with old and good K&R C. 😄

@seebs
Copy link

seebs commented Nov 14, 2018

sure there is, it hasn't got prototypes.

@cptnapalm
Copy link

M6: 45.8 GB of RAM and got me a load average of 64.08. Impressive!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment