#!/usr/bin/nawk -f $0
BEGIN {
  if(ARGC < 3) {
    frame_delay = 0.02;
    num_frames = 100
    num_pyrms = 1
    period = 2
  }
  else {
    num_frames  = ARGV[1];
    num_pyrms   = ARGV[2];
    frame_delay = ARGV[3];
  }
  pi2= 2*atan2(0,-1);

  verts = 19
  num_verts = num_pyrms * verts

  g[0] = 0; g[1] = 1; g[2] = 2; g[3] = 3; g[4] = 4; g[5] = 5;
  g[6] = 6; g[7] = 7; g[8] = 8; g[9] = 9; g[10] = 10; g[11] = 11; 
  g[12] = 12; g[13] = 13; g[14] = 14; g[15] = 15; g[16] = 16; g[17] = 17; 
  g[18] = 0;

  v[0,0] = -3; v[0,1] = -3; v[0,2] = -4;
  v[1,0] = 4; v[1,1] = -3; v[1,2] = -4;
  v[2,0] = 4; v[2,1] = 3; v[2,2] = -4;
  v[3,0] = 4; v[3,1] = 3; v[3,2] = 3;
  v[4,0] = 4; v[4,1] = -4; v[4,2] = 3;
  v[5,0] = 4; v[5,1] = -4; v[5,2] = -3;
  v[6,0] = -3; v[6,1] = -4; v[6,2] = -3;
  v[7,0] = -3; v[7,1] = -4; v[7,2] = 4;
  v[8,0] = 3; v[8,1] = -4; v[8,2] = 4;
  v[9,0] = 3; v[9,1] = 3; v[9,2] = 4;
  v[10,0] = -4; v[10,1] = 3; v[10,2] = 4;
  v[11,0] = -4; v[11,1] = -3; v[11,2] = 4;
  v[12,0] = -4; v[12,1] = -3; v[12,2] = -3;
  v[13,0] = -4; v[13,1] = 4; v[13,2] = -3;
  v[14,0] = -4; v[14,1] = 4; v[14,2] = 3;
  v[15,0] = 3; v[15,1] = 4; v[15,2] = 3;
  v[16,0] = 3; v[16,1] = 4; v[16,2] = -4;
  v[17,0] = -3; v[17,1] = 4; v[17,2] = -4;

  printf("%d %d %f\n",num_frames,num_verts,frame_delay); 

  for(k=0;k<num_frames;k++) {
    t1 = (k/num_frames)*pi2;
    t3 = (k/(num_frames-1))*pi2;

    for(j=0;j<num_verts;j++) {
      t2 = (j/num_verts)*pi2;
      t4 = (j/(num_verts-1))*pi2;

      x = v[g[(j%verts)],0];
      y = v[g[(j%verts)],1];
      z = v[g[(j%verts)],2];

      scale((int(j/verts)+1)*(verts/num_verts));
      scale(0.25);
      rotz(pi2/0.3)
      rotx(1.06)
      roty(t3)
      red   = cos(t4) * 127 + 128
      green = sin(t4) * 127 + 128
      blue  = int((j/(num_verts-1)) * 255)

      printf("%d %f %f %f %02x%02x%02x\n",flags,x,y,z,blue,green,red);
    }
  }
  exit(0);
}

function rotx(angle) {
  y1 = y
  y =  cos(angle) * y1 + sin(angle) * z
  z = -sin(angle) * y1 + cos(angle) * z
}

function roty(angle) {
  z1 = z
  z =  cos(angle)*z1 + sin(angle) * x
  x = -sin(angle)* z1 + cos(angle) * x
}

function rotz(angle) {
  x1 = x
  x =  cos(angle) * x1 + sin(angle) * y
  y = -sin(angle) * x1 + cos(angle) * y
}

function scale(value) {
  x = x * value;
  y = y * value;
  z = z * value;
}
