namespace; include rt; // the mandelbrot set in 16.16 fixed point. a word is 64 bits, so the product // fits and the point sits 16 bits down. the jit inlines it as a multiply and a shift. enhance fmul(a, b) micro "a = mul a b; a = sar a 16"; const ONE = 65536, W = 72, H = 28, MAXIT = 40; let shades = " .,:;-=+*%#@"; int escape(int cx, int cy) { int x = 0, y = 0, i = 0; while (i < MAXIT) { int xx = fmul(x, x), yy = fmul(y, y); if (xx + yy > 4 * ONE) return i; y = 2 * fmul(x, y) + cy; x = xx - yy + cx; i++; } return MAXIT; } void main() { int row = 0; while (row < H) { int col = 0, cy = (row * 2400 / (H - 1) - 1200) * ONE / 1000; while (col < W) { int cx = (col * 3200 / (W - 1) - 2200) * ONE / 1000, it = escape(cx, cy); rt.put_c(it == MAXIT ? '@' : shades[it * 11 / MAXIT]); col++; } rt.put_c(10); row++; } return; }