#macro KWR_Level(Step,Transform)
  #if(Step>0)
    #local Counter=0;
    #while (Counter<N)
      KWR_Level(Step-1,transform{transform{MT[Counter]}transform{Transform}})
      #local Counter=Counter+1;
    #end
    Branch()
  #else
    Leaf()
  #end
#end

// ============================================
// #macro KWR_3D(Entries3D,NSteps)
//
// Entries3D - array of parameters of 3D matrix
// Entries3D[N][4][3] - for N matrices
//
// x' = a x + b y + c z + d
//       1     1     1     1
//
// y' = a x + b y + c z + d
//       2     2     2     2
//
// z' = a x + b y + c z + d
//       3     3     3     3
//
// ============================================
#macro KWR_3D(Entries3D,NSteps)
  #local MT=array[dimension_size(Entries3D,1)]
  #local N=dimension_size(Entries3D,1);
  #local Counter=0;
  #while (Counter<N)
    #local MT[Counter]=transform{matrix<
      Entries3D[Counter][0][0],Entries3D[Counter][0][1],Entries3D[Counter][0][2]
      Entries3D[Counter][1][0],Entries3D[Counter][1][1],Entries3D[Counter][1][2]
      Entries3D[Counter][2][0],Entries3D[Counter][2][1],Entries3D[Counter][2][2]
      Entries3D[Counter][3][0],Entries3D[Counter][3][1],Entries3D[Counter][3][2]
    >}
    #local Counter=Counter+1;
  #end
  KWR_Level(NSteps,transform{translate 0})
#end

// ============================================
// #macro KWR_2D(Entries2D,NSteps)
//
// Entries2D - array of parameters of 2D matrix
// Entries2D[N][3][2] - for N matrices
//
//     x'=ax+by+c               y'=dx+ey+f
// a=Entries2D[n][0][0]    d=Entries2D[n][0][1]
// b=Entries2D[n][1][0]    e=Entries2D[n][1][1]
// c=Entries2D[n][2][0]    f=Entries2D[n][2][1]
//
// ============================================
#macro KWR_2D(Entries2D,NSteps)
  #local Entries3D=array[dimension_size(Entries2D,1)][4][3]
  #local Counter=0;
  #while (Counter<dimension_size(Entries2D,1))
    #local Entries3D[Counter][0][0]=Entries2D[Counter][0][0];
    #local Entries3D[Counter][1][0]=Entries2D[Counter][1][0];
    #local Entries3D[Counter][2][0]=0;
    #local Entries3D[Counter][3][0]=Entries2D[Counter][2][0];
    #local Entries3D[Counter][0][1]=Entries2D[Counter][0][1];
    #local Entries3D[Counter][1][1]=Entries2D[Counter][1][1];
    #local Entries3D[Counter][2][1]=0;
    #local Entries3D[Counter][3][1]=Entries2D[Counter][2][1];
    #local Entries3D[Counter][0][2]=0;
    #local Entries3D[Counter][1][2]=0;
    #local Entries3D[Counter][2][2]=1;
    #local Entries3D[Counter][3][2]=0;
    #local Counter=Counter+1;
  #end
  KWR_3D(Entries3D,NSteps)
#end

#declare P1=array[3][3][2]{
  {{0.500,0.000},{0.000,0.500},{0.0000,0.0000}},
  {{0.500,0.000},{0.000,0.500},{0.5000,0.0000}},
  {{0.500,0.000},{0.000,0.500},{0.0000,0.5000}},
}

#declare P2=array[3][3][2]{
  {{0.000,0.500},{-.500,0.000},{0.5000,0.0000}},
  {{0.000,-.500},{0.500,0.000},{0.5000,0.5000}},
  {{0.500,0.000},{0.000,0.500},{0.2500,0.5000}},
}

#declare P3=array[3][3][2]{
  {{0.000,-.577},{0.577,0.000},{0.0951,0.5893}},
  {{0.000,-.577},{0.577,0.000},{0.4413,0.7893}},
  {{0.000,-.577},{0.577,0.000},{0.0952,0.9893}},
}

#declare P4=array[3][3][2]{
  {{0.000,1.000},{-.333,0.000},{0.3333,0.0000}},
  {{0.333,0.000},{0.000,0.333},{0.3333,0.6666}},
  {{0.000,-1.00},{0.333,0.000},{0.6666,0.9999}},
}

#declare P5=array[3][3][2]{
  {{0.387,0.430},{0.430,-.387},{0.2560,0.5220}},
  {{0.441,-.009},{-.091,-.322},{0.4219,0.5059}},
  {{-.468,-.113},{0.020,0.015},{0.4000,0.4000}},
}

#declare P6=array[4][3][2]{
  {{0.255,0.000},{0.000,0.255},{0.3726,0.6714}},
  {{0.255,0.000},{0.000,0.255},{0.1146,0.2232}},
  {{0.255,0.000},{0.000,0.255},{0.6306,0.2232}},
  {{0.370,0.642},{-.642,0.370},{0.6356,-.0061}},
}

#declare P7=array[5][3][2]{
  {{0.382,0.000},{0.000,0.382},{0.3072,0.6190}},
  {{0.382,0.000},{0.000,0.382},{0.6033,0.4044}},
  {{0.382,0.000},{0.000,0.382},{0.4920,0.0595}},
  {{0.382,0.000},{0.000,0.382},{0.1253,0.0595}},
  {{0.382,0.000},{0.000,0.382},{0.0139,0.4044}},
}

#declare P8=array[5][3][2]{
  {{0.195,0.344},{-.488,0.443},{0.4431,0.2452}},
  {{0.462,-.252},{0.414,0.361},{0.2511,0.5692}},
  {{-.058,0.453},{-.070,-.111},{0.5976,0.0969}},
  {{-.035,-.469},{0.070,-.022},{0.4884,0.5069}},
  {{-.637,0.000},{0.000,0.501},{0.8562,0.2513}},
}

#declare P9=array[4][3][2]{
  {{0.849,-.037},{0.037,0.849},{0.0750,0.1830}},
  {{0.197,0.226},{-.226,0.197},{0.4000,0.0400}},
  {{-.150,0.260},{0.283,0.237},{0.5750,-.0840}},
  {{0.001,0.000},{0.000,0.160},{0.5000,0.0000}},
}