// L'auteur original de ce fichier est Jonathan Godin.
// jonathan.godin@umontreal.ca
// Il est permis de réutiliser ou de modifier 
// le code à condition de laisser cet entête.
//
// Itération de p(z) = z + z^2 
//
// Il faut compiler le programme comme suit
//
//$ gcc -std=gnu99 -O3 dendrite.c -o dendrite -lm 
//
// Tout dépendant de la taille de W et H, il est possible
// qu'il faille avoir accès à plus de mémoire pour la taille
// des tableaux. Dans ce cas, on peut utiliser
//
//$ ulimit -s unlimited
//
#define _GNU_C_SOURCE
#include <stdio.h>
#include <math.h>
#include <complex.h>
#define D(z1,z2) cabs((z1)-(z2))
#define MIN(a,b) ((a<b)?(a):(b))
#define W 1200 //Résolution
#define H 1200 //Résolution
#define w W/2
#define h H/2
#define C 1.5 //Changement d'échelle; on voit le carré
  //[-C/2,C/2]x[-C/2,C/2]
/* gcc -std=gnu99 -O3 z2para.c -o z2para -lm */

void inline couleur(int iter, char rgb[W][H][3], int lig, int col );

int main ()
{
    char area[W][H][3];
    complex double z,
        r = 0 + 0*I;
 
    for ( int im=-h; im<h; im++ )
        for ( int re=-w, iter; re<w; iter=0, re++ )
        {
            z = (C*re/((float) w)) + (C*im/((float) h))*I;
            while( D(z,r) < 1000 && iter++ < 100) 
              z = z + z*z;
            couleur(iter, area, im+h, re+ w);
        }

    /* Produire l'image dans un ppm */
    FILE *out = fopen ( "z2para.ppm", "w" );
    fprintf ( out, "P6 %d %d 255\n",W,H );
    fwrite ( area, 3, W*H, out );
    fclose ( out );
}

void inline couleur(int iter, char rgb[W][H][3], int lig, int col ){
  int r,g,b;
  if(iter < 2){
    r = 250;
    g = 30;
    b = 100;
  }
  else if(iter < 3){
    r = 230;
    g = 30;
    b = 100;
  }
  else if(iter < 4){
    r = 210;
    g = 40;
    b = 100;
  }
  else if(iter < 5){
    r = 255;
    g = 110;
    b = 0;
  }
  else if(iter < 6){
    r = 230;
    g = 255;
    b = 0;
  }
  else if(iter < 7){
    r = 130;
    g = 240;
    b = 110;
  }
  else if(iter < 8){
    r = 0;
    g = 255;
    b = 20;
  }
  else if(iter < 9){
    r = 0;
    g = 200;
    b = 255;
  }
  else if(iter < 10){
    r = 10;
    g = 125;
    b = 175;
  }
  else if(iter < 20){
    r = 50;
    g = 0;
    b = 250;
  }
  else if(iter < 40){
    r = 140;
    g = 80;
    b = 180;
  }
  else if(iter < 70){
    r = 250;
    g = 0;
    b = 240;
  }
  else if(iter < 90){
    r = 250;
    g = 0;
    b = 0;
  }
  else if(iter >=90){
    r = 0;
    g = 0;
    b = 0;
  }
  rgb[lig][col][0] = r; rgb[lig][col][1] = g; rgb[lig][col][2] = b;
};