/*
 *
 *
 */

#include <stdlib.h>
#include <stdio.h>
#include <assert.h>
#include <math.h>
#include "mandelbrot.h"

#define MAX_ESCAPESTEPS 256

int main (int argc, char *argv[]) {

   double xcentre = 0;
   double ycentre = 0;
   int z = 8;

   double pixeldistance = pow(2, -z);

   // start from the top left of the region
   double firstX = (xcentre - 256); //* pixeldistance;
   double lastX = (xcentre + 256);
   double firstY = (ycentre + 256);
   double lastY = (ycentre - 256);
   double middleX = (firstX + lastX);
   double middleY = (firstY + lastY);
 
   printf("midpoint: (%lf,%lf)\n", middleX, middleY);  
   printf("firstX = %lf\n", firstX);
   printf("lastX = %lf\n", lastX);
   printf("firstY = %lf\n", firstY);
   printf("lastY = %lf\n", lastY);
   
   int x;
   for (x = -256; x <= 256; x++) {
      double y;

    //  printf("%i\n", escapeSteps(x, y));


      for (y = -256*pixeldistance; y <= 256*pixeldistance; y+=pixeldistance) {
         printf("(%f,%f)", x, y);

//         printf("%i\n", escapeSteps(x, y));

         if (escapeSteps(x, y) < MAX_ESCAPESTEPS) {
            printf("*");
         } else {
            printf(" ");                   
         }
      }
      printf("\n");
   }

   return EXIT_SUCCESS;
}

int escapeSteps (double x, double y) {
   int steps = 0;
   double real = 0;
   double imaginary = 0;
   double modSquared = ((real * real) + (imaginary * imaginary));
   double newImaginary;
   double newReal;
   while (modSquared < 4 && steps < MAX_ESCAPESTEPS) {
      newImaginary = (2 * imaginary * real)+y;
      newReal = ((real * real) - (imaginary * imaginary))+x;
      imaginary = newImaginary;
      real = newReal;
      modSquared = ((newReal * newReal) + (newImaginary * newImaginary));
      steps++;
   }
 
   return steps;
}