//mandelbrot

#include <iostream>
#include <IFM/window>

int main()
{
  ifm::Wstream w(1000,600,"Mandelbrot");
  int pixelxmin= -500;
  int pixelxmax= 500;
  double grid = 0.01;
  int pixelymin= -300;
  int pixelymax= 300;
  int bound = 3000; //change this if you PC ist slow
  double changex = 0;
  double changey =0;
  double x = 0;
  double y = 0;
  int iter = 0;
  double value = 0;
  double xt;
  double yt;
  register double cy;
  register double cx;
  std::cout << "Die Mandelbrotmenge:\nZum zoomen das Rechteck ziehen\n"
	    << "(rechts-oben nach links-unten um nicht zu spiegeln)";
  for(;;){
    for(int positionx = pixelxmin;
	positionx<pixelxmax;positionx++)
      {
	cx = grid * positionx+changex;
	w << ifm::flush;
	for(int positiony = pixelymin;
	    positiony<pixelymax;positiony++)
	  {
	    cy = grid * positiony+changey;
	    x = 0;
	    y = 0;
	    iter = 0;
	    value = 0;
	    while (value <= 2 && iter < bound)
	      {
		xt = x*x-y*y + cx;
		yt = 2*x*y+cy;
		x = xt;
		y = yt;
		iter++;
		value = x*x +y*y;
	      }
	    if(value <= 2)
	      w<<ifm::black 
	       << ifm::Point(positionx-pixelxmin,
			     positiony-pixelymin);
	       
	  }
      }
    w<<ifm::red<<ifm::Line(100,500,200,500)
     <<ifm::Line(450,300,550,300)<<ifm::Line(500,250,500,350)<<ifm::flush;
    ifm::Rectangle zoom;
    w>>zoom;
    changex = changex -(500-(zoom.x1()+zoom.x2())/2)*grid;
    changey = changey -(300-(zoom.y1()+zoom.y2())/2)*grid;
    grid *= (zoom.x2()-zoom.x1())/1000.0;
    if (grid<0)grid = -grid;
    std::cout << "Eckdaten:\nLange der roten Linie: "<<100*grid<<"\n"
	      <<"Mittelpunkt:    ("<<-changex<<"/"<<-changey<<")\n"
	      <<"";
    if(grid < 10e-22&& grid > -10e-22)grid = 0.01;
    w<<ifm::clear;
  }
  return 0;
}