#include "stdio.h"
#include "stdlib.h"
#include "math.h"
#include <sys/types.h>
#include <time.h>
#include "paulslib.h"
#include "bitmaplib.h"

#define NX 1000
#define NY 1000
#define N 10000000
#define SCALE (NX / 2)

int main(int argc,char **argv)
{
   int c,i,ix,iy;
   double sq15 = sqrt(1.5);
   double x=0.5,y=0.75,xa,xb,xc,ya,yb,yc,a,b,xr,yr,x2,y2;
   double xmin=1e32,ymin=1e32,xmax=-1e32,ymax=-1e32;
   BITMAP *image,white={255,255,255},black={0,0,0};
   BITMAP colour;
   char fname[64];
   FILE *fptr;

   if (argc < 1) {
      fprintf(stderr,"Usage: %s \n",argv[0]);
      exit(-1);
   }

   /* Create the image */
   image = CreateBitmap(NX,NY);
   EraseBitmap(image,NX,NY,white);
   srand(time(NULL));

   for (i=0;i<N;i++) {
      xr = cos(TWOPI/3)*x - sin(TWOPI/3)*y;
      yr = cos(TWOPI/3)*y + sin(TWOPI/3)*x;
      x2 = xr*xr - yr*yr;
      y2 = 2 * xr * yr;
      switch (rand() % 3) {
      case 0:
         xa = (SQRT3-1)*x+1;
         ya = (SQRT3-1)*y;
         xb = -x+(SQRT3+1);
         yb = -y;
         xc = xb/(xb*xb+yb*yb);
         yc = -yb/(xb*xb+yb*yb);
         break;
      case 1:
         xa = (SQRT3-1)*xr+1;
         ya = (SQRT3-1)*yr;
         xb = -xr+(SQRT3+1);
         yb = -yr;
         xc = xb/(xb*xb+yb*yb);
         yc = -yb/(xb*xb+yb*yb);
         break;
      case 2:
         xa = (SQRT3-1)*x2+1;
         ya = (SQRT3-1)*y2;
         xb = -x2+(SQRT3+1);
         yb = -y2;
         xc = xb/(xb*xb+yb*yb);
         yc = -yb/(xb*xb+yb*yb);
         break;
      }
      x = (xa*xc - ya*yc) * sq15;
      y = (xa*yc + xc*ya) * sq15;
      xmin = MIN(x,xmin);
      ymin = MIN(y,ymin);
      xmax = MAX(x,xmax);
      ymax = MAX(y,ymax);
      if ((i % (N/100)) == 0)
         fprintf(stderr,"%d (%g,%g) - (%g,%g)\n",i,xmin,ymin,xmax,ymax);
      if (i < 100)
         continue;
      ix = x * SCALE + NX/2;
      iy = y * SCALE + NY/2;
      DrawPixel(image,NX,NY,ix,iy,black);
   }

   /* Write the image */
   sprintf(fname,"mauldin.tga");
   if ((fptr = fopen(fname,"w")) == NULL) {
      fprintf(stderr,"Unable to open bitmap file\n");
      exit(0);
   }
   WriteBitmap(fptr,image,NX,NY,12);
   fclose(fptr);
   DestroyBitmap(image);
}