//W:Mute 2004 - Abstract
//www.wmute.org

float[][] points={{0,0,1},
{0.894427,0,0.447214},
{0.276393,0.850651,0.447214},
{-0.723607,0.525731,0.447214},
{-0.723607,-0.525731,0.447214},
{0.276393,-0.850651,0.447214},
{0.723607,0.525731,-0.447214},
{-0.276393,0.850651,-0.447214},
{-0.894427,0,-0.447214},
{-0.276393,-0.850651,-0.447214},
{0.723607,-0.525731,-0.447214},
{0,0,-1},
{0.360729,0,0.932671},
{0.672883,0,0.739749},
{0.111471,0.343074,0.932671},
{0.207932,0.63995,0.739749},
{-0.291836,0.212031,0.932671},
{-0.544374,0.395511,0.739749},
{-0.291836,-0.212031,0.932671},
{-0.544374,-0.395511,0.739749},
{0.111471,-0.343074,0.932671},
{0.207932,-0.63995,0.739749},
{0.784354,0.343074,0.516806},
{0.568661,0.63995,0.516806},
{-0.0839038,0.851981,0.516806},
{-0.432902,0.738584,0.516806},
{-0.83621,0.183479,0.516806},
{-0.83621,-0.183479,0.516806},
{-0.432902,-0.738585,0.516806},
{-0.0839036,-0.851981,0.516806},
{0.568661,-0.63995,0.516806},
{0.784354,-0.343074,0.516806},
{0.964719,0.212031,0.156077},
{0.905103,0.395511,-0.156077},
{0.0964608,0.983023,0.156077},
{-0.0964609,0.983023,-0.156077},
{-0.905103,0.395511,0.156077},
{-0.964719,0.212031,-0.156077},
{-0.655845,-0.738585,0.156077},
{-0.499768,-0.851981,-0.156077},
{0.499768,-0.851981,0.156077},
{0.655845,-0.738584,-0.156077},
{0.964719,-0.212031,0.156077},
{0.905103,-0.39551,-0.156077},
{0.499768,0.851981,0.156077},
{0.655845,0.738584,-0.156077},
{-0.655845,0.738584,0.156077},
{-0.499768,0.851981,-0.156077},
{-0.905103,-0.395511,0.156077},
{-0.964719,-0.212031,-0.156077},
{0.0964609,-0.983023,0.156077},
{-0.0964608,-0.983023,-0.156077},
{0.432902,0.738584,-0.516806},
{0.0839037,0.851981,-0.516806},
{-0.568661,0.63995,-0.516806},
{-0.784354,0.343074,-0.516806},
{-0.784354,-0.343074,-0.516806},
{-0.568661,-0.63995,-0.516806},
{0.0839038,-0.851981,-0.516806},
{0.432902,-0.738584,-0.516806},
{0.83621,-0.183479,-0.516806},
{0.83621,0.183479,-0.516806},
{0.291836,0.212031,-0.932671},
{0.544374,0.395511,-0.739749},
{-0.111471,0.343074,-0.932671},
{-0.207932,0.63995,-0.739749},
{-0.360729,0,-0.932671},
{-0.672883,0,-0.739749},
{-0.111471,-0.343074,-0.932671},
{-0.207932,-0.63995,-0.739749},
{0.291836,-0.212031,-0.932671},
{0.544374,-0.39551,-0.739749},
{0.479506,0.348381,0.805422},
{-0.183155,0.563693,0.805422},
{-0.592702,0,0.805422},
{-0.183155,-0.563693,0.805422},
{0.479506,-0.348381,0.805422},
{0.985456,0,-0.169933},
{0.304523,0.937224,-0.169933},
{-0.79725,0.579236,-0.169933},
{-0.79725,-0.579236,-0.169933},
{0.304523,-0.937224,-0.169933},
{0.79725,0.579236,0.169933},
{-0.304523,0.937224,0.169933},
{-0.985456,0,0.169933},
{-0.304522,-0.937224,0.169933},
{0.79725,-0.579236,0.169933},
{0.183155,0.563693,-0.805422},
{-0.479506,0.348381,-0.805422},
{-0.479506,-0.348381,-0.805422},
{0.183155,-0.563693,-0.805422},
{0.592702,1.2143e-007,-0.805422}};

int num=92;
int alpha=150;
float scale;
float ax=0;
float ay=0;
float x,y,z,tx,ty,tz,mx,my,sf;
BImage bkg;

Triplets Ts;

void setup(){
  smooth();
  size(640,360);
  framerate(18);
  stroke(0,0,0,alpha);
  scale=120.0;
  Ts=new Triplets();
  bkg = loadImage("D4MD_bkg.jpg");
  background(bkg);
}

void loop(){
  background(bkg);
  translate(200,height/2.0,0);
  my=min(1.0,1.0*float(mouseX-width/2)/width);
  mx=min(1.0,1.0*float(mouseY-width/2)/height);
  ax+=mx*TWO_PI/500;
  ay+=my*TWO_PI/500;

  for (int i=0; i<num;i++){
    x=scale*points[i][0];
    y=scale*points[i][1];
    z=scale*points[i][2];
    tx=x*cos(ay)+y*sin(ax)*sin(ay)+z*cos(ax)*sin(ay);
    ty=y*cos(ax)-z*sin(ax);
    tz=-x*sin(ay)+y*sin(ax)*cos(ay)+z*cos(ax)*cos(ay);
    Ts.update(i,100,0,0,0,20,tx,ty,tz);

  }
  Ts.draw();
}

class Triplets {
  float[] cx = new float[num];
  float[] cy = new float[num];
  float[] lead = new float[num];
  float[] offset = new float[num];
  float[] wid = new float[num];
  float[] px = new float[num];
  float[] py = new float[num];
  float[] pz = new float[num];

  Triplets(){
    for (int i=0; i<num;i++){
      cx[i]=0;
      cy[i]=0;
      lead[i]=0;
      offset[i]=0;
      wid[i]=0;
      px[i]=0;
      py[i]=0;
      pz[i]=0;
    }
  }

  void update(int i, float cxx, float cyy, float ll, float oo, float ww, float pxx, float pyy, float pzz){
    cx[i]=cxx;
    cy[i]=cyy;
    lead[i]=ll;
    offset[i]=oo;
    wid[i]=ww;

    px[i]=pxx;
    py[i]=pyy;
    pz[i]=pzz;

  }

  void draw(){
    for (int i=0; i<num;i++){
      float sx=1.0;
      float sy=1.0;
      float sz=1.0;
      float sxx=cx[i]+(lead[i]+offset[i]);
      float syx=px[i];
      float sxy=cx[i]+lead[i]+offset[i];
      float syy=py[i];
      float sxz=cx[i]+lead[i]+offset[i];
      float syz=pz[i];

      if(abs(px[i]-py[i])<10){
        stroke(15*abs(px[i]-py[i]),0,0,alpha);
        bezier(sxx,syx,sxx-50,syx,-0.5*(sxy-50)+0.866*syy,-0.5*syy-0.866*(sxy-50),-0.5*sxy+0.866*syy,-0.5*syy-0.866*sxy);
        sx+=0.5;
        sy+=0.5;
      }
      if(abs(px[i]-pz[i])<10){
        stroke(15*abs(px[i]-pz[i]),0,0,alpha);
        bezier(sxx,syx,sxx-50,syx,-0.5*(sxz-50)-0.866*syz,-0.5*syz+0.866*(sxz-50),-0.5*sxz-0.866*syz,-0.5*syz+0.866*sxz);
        sx+=0.5;
        sz+=0.5;
      }
      if(abs(py[i]-pz[i])<10){
        stroke(15*abs(py[i]-pz[i]),0,0,alpha);
        bezier(-0.5*sxy+0.866*syy,-0.5*syy-0.866*sxy,-0.5*(sxy-50)+0.866*syy,-0.5*syy-0.866*(sxy-50),-0.5*(sxz-50)-0.866*syz,-0.5*syz+0.866*(sxz-50),-0.5*sxz-0.866*syz,-0.5*syz+0.866*sxz);
        sy+=0.5;
        sz+=0.5;
      }
      float exx=cx[i]+lead[i]+offset[i]+sx*wid[i];
      float eyx=px[i];
      float exy=cx[i]+lead[i]+offset[i]+sy*wid[i];
      float eyy=py[i];
      float exz=cx[i]+lead[i]+offset[i]+sz*wid[i];
      float eyz=pz[i];

      line(sxx,syx,exx,eyx);

      line(-0.5*sxy+0.866*syy,-0.5*syy-0.866*sxy,-0.5*exy+0.866*eyy,-0.5*eyy-0.866*exy);

      line(-0.5*sxz-0.866*syz,-0.5*syz+0.866*sxz,-0.5*exz-0.866*eyz,-0.5*eyz+0.866*exz);

    }
  }

}