/* Sequential Mandelbrot program */


#include <X11/Xlib.h>
#include <X11/Xutil.h>
#include <X11/Xos.h>
#include <stdio.h>
#include <string.h>
#include <math.h>

#define		X_RESN	800      /* x resolution */
#define		Y_RESN	800       /* y resolution */
#define       MAX_COLORS 256
/* #define       MAX_COLORS 100 */

typedef struct complextype
	{
        float real, imag;
	} Compl;


int main ()
{
	Window		win;                            /* initialization for a window */
	unsigned
	int             width, height,                  /* window size */
                        x, y,                           /* window position */
                        border_width,                   /*border width in pixels */
                        display_width, display_height,  /* size of screen */
                        screen;                         /* which screen */

	char            *window_name = "Mandelbrot Set", *display_name = NULL;
	GC              gc;
	unsigned
	long		valuemask = 0;
	XGCValues	values;
	Display		*display;
	XSizeHints	size_hints;
	Pixmap		bitmap;
	XPoint		points[800];
	FILE		*fp, *fopen ();
	char		str[100];
	
	Colormap colormap;
	XColor colors[MAX_COLORS];
	
	XSetWindowAttributes attr[1];

       /* Mandlebrot variables */
        int i, j, k;
        Compl	z, c;
        float	lengthsq, temp;
       
	/* connect to Xserver */

	if (  (display = XOpenDisplay (display_name)) == NULL ) {
	   fprintf (stderr, "drawon: cannot connect to X server %s\n",
				XDisplayName (display_name) );
	exit (-1);
	}
	
	/* get screen size */

	screen = DefaultScreen (display);
	display_width = DisplayWidth (display, screen);
	display_height = DisplayHeight (display, screen);

	/* set window size */

	width = X_RESN;
	height = Y_RESN;

	/* set window position */

	x = 0;
	y = 0;

        /* create opaque window */

	border_width = 4;
	win = XCreateSimpleWindow (display, RootWindow (display, screen),
				x, y, width, height, border_width, 
				BlackPixel (display, screen), WhitePixel (display, screen));

	size_hints.flags = USPosition|USSize;
	size_hints.x = x;
	size_hints.y = y;
	size_hints.width = width;
	size_hints.height = height;
	size_hints.min_width = 300;
	size_hints.min_height = 300;
	
	XSetNormalHints (display, win, &size_hints);
	XStoreName(display, win, window_name);
	
	colormap = DefaultColormap(display, DefaultScreen(display));
	
	/*if you want regist colors by value, please do this
	 * for(i=0; i<MAX_COLORS; i++)
	 * {
	 *	colors[i].red = 200*i;   //value from 0~65535
	 *	colors[i].green = 3000+i*200; //value from 0~65535
	*	colors[i].blue = 0; //value from 0~65535
	*	XAllocColor(dispay, colormap, &colors[i]);
	* }
	*/
	
	/*register colors*/
	for(i=0; i<MAX_COLORS; i++)
	{
	    if(i<10)
	    XAllocNamedColor(display, colormap, "red",&colors[i], &colors[i]);
	    else if(i<20)
	    XAllocNamedColor(display, colormap, "yellow",&colors[i], &colors[i]);
	    else if(i<30)
	    XAllocNamedColor(display, colormap, "green",&colors[i], &colors[i]);
	    else if(i<50)
	    XAllocNamedColor(display, colormap, "brown",&colors[i], &colors[i]);
	    else if(i<60)
	    XAllocNamedColor(display, colormap, "blue",&colors[i], &colors[i]);
	    else if(i<70)
	    XAllocNamedColor(display, colormap, "red",&colors[i], &colors[i]);
	    else if(i<80)
	    XAllocNamedColor(display, colormap, "blue",&colors[i], &colors[i]);
	    else if(i<90)
	    XAllocNamedColor(display, colormap, "green",&colors[i], &colors[i]);
	    else 
	    XAllocNamedColor(display, colormap, "red",&colors[i], &colors[i]);
	    
	    
	} 

        /* create graphics context */

	gc = XCreateGC (display, win, valuemask, &values);

	XSetBackground (display, gc, WhitePixel (display, screen));
	//XSetForeground (display, gc, BlackPixel (display, screen));
	XSetLineAttributes (display, gc, 1, LineSolid, CapRound, JoinRound);

	attr[0].backing_store = Always;
	attr[0].backing_planes = 1;
	attr[0].backing_pixel = BlackPixel(display, screen);

	XChangeWindowAttributes(display, win, CWBackingStore | CWBackingPlanes | CWBackingPixel, attr);

	XMapWindow (display, win);
	XSync(display, 0);
      	 
        /* Calculate and draw points.  The beginning of the actual
	   mandelbrot code */

        for(i=0; i < X_RESN; i++) 
        for(j=0; j < Y_RESN; j++) {

          z.real = z.imag = 0.0;
          c.real = ((float) j - 400.0)/200.0;               /* scale factors for 800 x 800 window */
	  c.imag = ((float) i - 400.0)/200.0;
          k = 0;

          do  {                                             /* iterate for pixel color */

            temp = z.real*z.real - z.imag*z.imag + c.real;
            z.imag = 2.0*z.real*z.imag + c.imag;
            z.real = temp;
            lengthsq = z.real*z.real+z.imag*z.imag;
            k++;

          } while (lengthsq < 4.0 && k < MAX_COLORS);


	XSetForeground (display, gc, colors[k].pixel);
	XDrawPoint (display, win, gc, j, i);
	

        }
	 
	XFlush (display);
	sleep (10); 

	/* Program Finished */

}