gambas-source-code/benchmark/mandelbrot.py

#!/usr/bin/python

import sys

MAXITER = 50
LIMIT = 4
out = sys.stdout;

def mandelbrot(w, h) :

  xmin = -1.5
  ymin = -1
  invN = 2.0 / w

  checknext = True

  for y in range(h) :
    Ci = y * invN + ymin

    for x in range(w) :
      Zr = 0.0
      Zi = 0.0
      Tr = 0.0
      Ti = 0.0
      Cr = x * invN + xmin
      if (checknext) :
	for k in range(MAXITER) :
	  Zi = 2 * Zr * Zi + Ci
	  Zr = Tr - Ti + Cr
	  Ti = Zi * Zi
	  Tr = Zr * Zr
	  if (Tr + Ti) > LIMIT :
	    break
	if k == MAXITER :
	  out.write("1")
	else :
	  out.write("0")
	  checknext = False
      else :
	for k in range(MAXITER) :
	  Zi = 2 * Zr * Zi + Ci
	  Zr = Tr - Ti + Cr
	  Ti = Zi * Zi
	  Tr = Zr * Zr
	if (Tr + Ti) < LIMIT :
	  out.write("1")
	else :
	  out.write("0")
	  checknext = True

    out.write("\n")

for i in range(20) :
  mandelbrot(200, 200)