/*----------------------------------------------------------*/
/*                                                          */
/*         AUTHOR : Eric VIOLARD                            */
/*         E-MAIL : violard@icps.u-strasbg.fr               */
/*         ORGANISM : Université Louis Pasteur (Strasbourg) */
/*         CREATION : 23/11/02                              */
/*                                                          */
/* ---------------------------------------------------------*/
  
/*  

     Programme permettant de dessiner l'ensemble de Mandelbrot.
 
     Le programme demande à l'utilisateur les dimensions de la fenêtre graphique.

     (version qui utilise les bibliothèques complexe et intervalle)
*/

#include <stdio.h> /* pour utiliser les fonctions printf et scanf */
#include <graph.h> /* pour utiliser les fonctions graphiques */ 
#include <complexe.h> /* pour utiliser les fonctions sur les nombres complexes */
#include <intervalle.h> /* pour utiliser les fonctions sur les intervalles */

/* --- fonctions utiles */

float min(float a, float b)
     /* minimum de 2 réels */
{
  float r;

  r=a;
  if(b<a) r=b;
  return r;
}

float max(float a, float b)    
    /* maximum de 2 réels */
{
  float r;

  r=b;
  if(a>b) r=a;
  return r;
}

#define N 200 /* borne utilisée pour déterminer la convergence des suites */


main()
{
  

  /* déclarations */
  
  int largeur,hauteur; /* dimensions de la fenêtre */
  complexe zmin=cons_complexe(-2.05,-1.3),
           zmax=cons_complexe(0.55,1.3); /* défini l'intervalle des nombres complexes considérés */
  int X,Y; /* les coordonnées d'un pixel de la fenêtre graphique */
  complexe c; /* le nombre complexe dans l'intervalle [zmin,zmax] et qui correpond au pixel (X,Y)  */
  complexe z;  /* termes de la suite (zn) */
  int n; /* indice des termes de la suite (zn) */


  /* --- initialisation de la fenêtre graphique */

  printf("Entrer les dimensions de la fenêtre graphique (largeur x hauteur) en nombre de pixels : "); 
  scanf("%d %d",&largeur,&hauteur); 
  initgraph(largeur,hauteur);
  cleargraph();


  /* --- tracé de l'ensemble --- */ 

  setcolor(0.,1.,0.); /* en vert */


  /* on parcours tous les pixels */

  for(X=0;X<largeur;X++)
  for(Y=0;Y<hauteur;Y++) 
    {

      /* --- construction de c, le nombre complexe correspondant au pixel (X,Y) --- */
      { 
      float x,y; /* partie réelle et partie imaginaire de c (variables temporaires) */

      x = thales((float)X,cons_intervalle(0.,(float)largeur),
                   cons_intervalle(partie_reelle(zmin),partie_reelle(zmax)));
      y = thales((float)Y,cons_intervalle(0.,(float)hauteur),
                   cons_intervalle(partie_imaginaire(zmin),partie_imaginaire(zmax)));
      c = cons_complexe(x,y);
      }
 
      /* --- calcul de termes de la suite (zn) --- */
      
      z = cons_complexe(0.,0.);      
      n = 0;

      while(module_complexe(z)<2 && n<=N)
        {
	  z = somme_complexe(produit_complexe(z,z),c);   
	  n = n+1;
        }

      
      /* --- illumination éventuelle du pixel (X,Y) --- */

      if(module_complexe(z)<2) putpixel(X,Y);
    }

  flushgraph(); /* pour forcer l'affichage */


  /* --- terminaisons --- */

  waitgraph(); /* attente d'un frappe au clavier */
  closegraph(); /* fermeture de la fenêtre graphique */

}