// Lutz Tautenhahn (c) 1995, 1999 
// Two simple fractal examples in one, try to understand how it works!
// The main idea is the recursive coordinate transformation by
// translation, rotation, stretching: 
// Parent ObjectSet={o_1, ..., o_N_Ob}
//   ->CoordinateTrafo(OperatorSet={o_1, ..., o_N_Op}) 
//     ->Child ObjectSet_1, ..., Child ObjectSet_N_Op
//       ... ->SubChildren  ... ->SubSubChildren
#include <graph.h>
#include <conio.h>                                           
#include <math.h>
#include <stdio.h>

#define Pi 3.14159265 
#define xmin -1.33333     //real x of screen x_min
#define xmax 1.33333      //real x of screen x_max
#define ymin 0.0 //-1.0   //real y of screen y_min
#define ymax 2.0 //1.0    //real y of screen y_max
#define N_Op 3   //5      //Number of Operations (coordinate transformations)
#define N_Ob 1   //4      //Number of Objects (=lines of the parent shape)                    


int ScreenX(double theRealX)
{
  return(int((theRealX-xmin)/(xmax-xmin)*640));
}  

int ScreenY(double theRealY)
{           
  return(int((ymax-theRealY)/(ymax-ymin)*480));
}  

struct Vector
{
  double x;
  double y;
  double z;
};       
         
double ScPr(Vector theVector1, Vector theVector2)
{ double aDouble;
  Vector aVec1 = theVector1;
  Vector aVec2 = theVector2;
  aDouble = aVec1.x * aVec2.x + aVec1.y * aVec2.y + aVec1.z * aVec2.z;
  return(aDouble);
}  

struct Object
{ 
  Vector V[2];
};

struct ObjectSet
{
  Object o[N_Ob];
};  
    
struct Operator
{
  Vector o[3];
};           

struct OperatorSet
{
  Operator o[N_Op];
};  
                  
void Operation(ObjectSet& aTempObSet, Operator aOp, ObjectSet aObSet)
{ for (int jj=0; jj<N_Ob; jj++)
  { for (int ii=0; ii<2; ii++)
    { aTempObSet.o[jj].V[ii].x=ScPr(aOp.o[0],aObSet.o[jj].V[ii]);
      aTempObSet.o[jj].V[ii].y=ScPr(aOp.o[1],aObSet.o[jj].V[ii]); 
      aTempObSet.o[jj].V[ii].z=1;
    }
  }       
}
                  
class Go      
{ public:     
    OperatorSet OpSet;
    ObjectSet ObSet, TempObSet; 
    short T; 
    Go* aGo;
    Go(OperatorSet aOpSet, ObjectSet aObSet, short aT);
};                              

Go::Go(OperatorSet aOpSet, ObjectSet aObSet, short aT)
{ OpSet=aOpSet;
  ObSet=aObSet;
  T=aT;
  for (int ii=0; ii<N_Ob; ii++)    
  { _moveto(ScreenX(aObSet.o[ii].V[0].x), ScreenY(aObSet.o[ii].V[0].y));
    _lineto(ScreenX(aObSet.o[ii].V[1].x), ScreenY(aObSet.o[ii].V[1].y)); 
//    _getch();
  }
  if (aT>0)
  { for (ii=0; ii<N_Op; ii++)
    { Operation(TempObSet, OpSet.o[ii], ObSet),      
      aGo = new Go(OpSet, TempObSet, T-1);  
      delete(aGo);
    }  
  }
}

void main( void )
{ ObjectSet X;
  OperatorSet A;     

  double Fi[N_Op]={-12*Pi/180,20*Pi/180, 6*Pi/180}; 
  double S[N_Op]={0.35, 0.3, 0.4};
  double DX[N_Op]={0, 0, 0};
  double DY[N_Op]={0.4, 0.7, 1};

  X.o[0].V[0].x=0;
  X.o[0].V[0].y=0;
  X.o[0].V[0].z=1;
  X.o[0].V[1].x=0;
  X.o[0].V[1].y=1;
  X.o[0].V[1].z=1;            
                    
/*
  double Fi[N_Op]={ 0, 0, 0, 0, 0}; 
  double S[N_Op]={ 0.33, 0.33, 0.33, 0.33, 0.33};
  double DX[N_Op]={0.67, 0, -0.67, 0, 0};
  double DY[N_Op]={0, 0.67, 0, -0.67, 0};

  X.o[0].V[0].x=-1;
  X.o[0].V[0].y=0;
  X.o[0].V[0].z=1;
  X.o[0].V[1].x=0;
  X.o[0].V[1].y=-1;
  X.o[0].V[1].z=1;            
  X.o[1].V[0].x=0;
  X.o[1].V[0].y=-1;
  X.o[1].V[0].z=1;
  X.o[1].V[1].x=1;
  X.o[1].V[1].y=0;
  X.o[1].V[1].z=1;            
  X.o[2].V[0].x=1;
  X.o[2].V[0].y=0;
  X.o[2].V[0].z=1;
  X.o[2].V[1].x=0;
  X.o[2].V[1].y=1;
  X.o[2].V[1].z=1;            
  X.o[3].V[0].x=0;
  X.o[3].V[0].y=1;
  X.o[3].V[0].z=1;
  X.o[3].V[1].x=-1;
  X.o[3].V[1].y=0;
  X.o[3].V[1].z=1;            
*/
  for (int ii=0; ii<N_Op; ii++)
  {
    A.o[ii].o[0].x = S[ii]*cos(Fi[ii]);
    A.o[ii].o[0].y = S[ii]*sin(Fi[ii]);
    A.o[ii].o[0].z = DX[ii];
    A.o[ii].o[1].x = S[ii]*(-sin(Fi[ii]));
    A.o[ii].o[1].y = S[ii]*cos(Fi[ii]);
    A.o[ii].o[1].z = DY[ii];
    A.o[ii].o[2].x = 0;
    A.o[ii].o[2].y = 0;
    A.o[ii].o[2].z = 1;
  }
  _setvideomode( _VRES16COLOR ); 
  delete(new Go(A, X, 5));
  _getch();    
  _setvideomode( _DEFAULTMODE );
}