#include <stdio.h>
#include <stdlib.h>
#include <omp.h>

int escapes(double cr, double ci, int it) {
  double zr = 0;
  double zi = 0;
  double zrtmp;
  int i;
  
  for(i=0; i<it; i++) {
    // z <- z^2 + c
    zrtmp = zr*zr - zi*zi + cr;
    zi = 2*zr*zi + ci;
    zr = zrtmp;
    if (zr*zr + zi*zi > 4) {
      return 1;
    }
  }
  return 0;
}

void mandel(double xmin, double xmax, int xstep, double ymin, double ymax, int ystep, int iters) {
  int yc;
  
  // array of string to store result
  char *m = (char *) malloc(ystep * (xstep + 1) * sizeof(char));

  #pragma omp parallel for  
  for(yc=0; yc<ystep; yc++) {
    double y = yc*(ymax-ymin)/ystep + ymin;
	int xc;
    for(xc=0; xc<xstep; xc++) {
      double x = xc*(xmax-xmin)/xstep + xmin;
      escapes(x, y, iters);
      if (escapes(x, y, iters)) {
        m[yc * (xstep + 1) + xc] = ' ';
      } else {
        m[yc * (xstep + 1) + xc] = 'X';   
      }
    }
    // add end of string
    m[yc * (xstep+1) + xstep] = '\0';  
    
  }
  
  for(yc=0; yc<ystep; yc++) {
	  printf("%s\n", &m[yc * (xstep+1)]);
  }
  
  free(m);
}

##</code><code>##

// mandelbrot.h
void mandel(double xmin, double xmax, int xstep, double ymin, double ymax, int ystep, int iters);

##</code><code>##

use warnings;
use Mandelbrot;
 
Mandelbrot::mandel(-2.0, 1.0, 256, -1.0, 1.0, 256, 100000); 

##</code><code>##

	LIBS      => ['-lgomp'], # e.g., '-lm'
	CCFLAGS   => '-fopenmp',

##</code><code>##

export OMP_NUM_THREADS=2