#include <stdio.h>
#include <math.h>
#include <string.h>
#include <fstream.h>
#include "screen.h"
#include <string.h>

float XPoint(float,float, float,float, float,float);
float YPoint(float,float, float,float, float,float);
void fractal(float,float,int,const float);
void fname(char [], const int);

const float PI=3.1415927;

Screen S(600,600);

int main(){
  ofstream scriptfile("makegif");
  ofstream mergegifs("mergegif");
  ofstream removal("gifremoval");
  char name[10];
  S.setrect(-1.5,-1.5,3,3);
  mergegifs << "gifsicle --delay 50 -O2 -l ";
  for (int i=0;i<18;i++){
    fname(name,i);
    fractal(0,0,7,PI*i/36);
    scriptfile << "convert " << name << ".xpm " << name << ".gif"<< endl;
    scriptfile << "rm " << name << ".xpm" << endl;
    mergegifs << name << ".gif ";
    removal << "rm " << name << ".gif" << endl;
    S.matrix2xpm(strcat(name,".xpm"));
    cout << "Done: Image " << i  << endl;
    S.clear();
  }
  mergegifs << "> anim.gif";
  return 0;
}

float XPoint(float xc,float yc, float theta,float scale, float x,float y){
  return ((x*cos(theta)-y*sin(theta))*scale+xc);
}

float YPoint(float xc,float yc, float theta,float scale, float x,float y){
  return ((x*sin(theta)+y*cos(theta))*scale+yc);
}

void fractal(float x, float y, int depth, const float ang){
  if (depth==0){
    S.plotxy(x,y);  
  } else {
    fractal(XPoint(cos(ang),sin(ang),ang,.33333333,x,y),YPoint(cos(ang),sin(ang),ang,.33333333,x,y),depth-1, ang);
    fractal(XPoint(.70710678*cos(ang+PI/4),.70710678*sin(ang+PI/4),ang,.33333333,x,y),YPoint(.70710678*cos(ang+PI/4),.70710678*sin(ang+PI/4),ang,.33333333,x,y),depth-1, ang);
    fractal(XPoint(cos(ang+PI/2),sin(ang+PI/2),ang,.33333333,x,y),YPoint(cos(ang+PI/2),sin(ang+PI/2),ang,.33333333,x,y),depth-1, ang);
    fractal(XPoint(.70710678*cos(ang+3*PI/4),.70710678*sin(ang+3*PI/4),ang,.33333333,x,y),YPoint(.70710678*cos(ang+3*PI/4),.70710678*sin(ang+3*PI/4),ang,.33333333,x,y),depth-1, ang);
    fractal(XPoint(cos(ang+PI),sin(ang+PI),ang,.33333333,x,y),YPoint(cos(ang+PI),sin(ang+PI),ang,.33333333,x,y),depth-1, ang);
    fractal(XPoint(.70710678*cos(ang+5*PI/4),.70710678*sin(ang+5*PI/4),ang,.33333333,x,y),YPoint(.70710678*cos(ang+5*PI/4),.70710678*sin(ang+5*PI/4),ang,.33333333,x,y),depth-1, ang);
    fractal(XPoint(cos(ang+3*PI/2),sin(ang+3*PI/2),ang,.33333333,x,y),YPoint(cos(ang+3*PI/2),sin(ang+3*PI/2),ang,.33333333,x,y),depth-1, ang);
    fractal(XPoint(.70710678*cos(ang+7*PI/4),.70710678*sin(ang+7*PI/4),ang,.33333333,x,y),YPoint(.70710678*cos(ang+7*PI/4),.70710678*sin(ang+7*PI/4),ang,.33333333,x,y),depth-1, ang);
  }
}

void fname(char name[], const int n){
  strcpy(name,"aa");;
  name[0]+=(n/26)%26;
  name[1]+=n%26;
}