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


int numanchors=5;
int numchains=15;
int numbeads=5;
float dotSiz=2.0;
float minLuminance=0.2;
float maxLuminance=1.0;
float weightFactor=1.0;
int weightToLevel=150;
chain2d[][] chains = new chain2d[numanchors][numchains];
chain2d[][] roamers = new chain2d[numanchors][numchains];
int[] level = new int[numanchors];
float[] weight= new float[numanchors];
int[] output = new int[1];
boolean first=true;
float[] luminance = new float[numanchors];
float lum;
int clr;
BFont font;
boolean thinking=true;
boolean redraw=true;
int counter=0;
BImage bkg;

void setup(){
  size(640,360);
  background(0);
  framerate(1);
  noStroke();
  ellipseMode(CENTER_RADIUS);
  output=new int[width*height];
  bkg=loadImage("D7AP_bkg.jpg");
  setall();
  reset();
}

void setall(){
  int xi, yi;
  float d;
  for(int a=0;a<numanchors;a++){
    level[a]=0;
    weight[a]=0.0;
    xi=int(random(0,width));
    yi=int(random(0,height));
    for(int i=0;i<numchains;i++) chains[a][i] = new chain2d(xi,yi,random(0,width),random(0,height),numbeads);
    for(int i=0;i<numchains;i++) roamers[a][i] = new chain2d(xi,yi,random(0,width),random(0,height),numbeads);
    luminance[a]=random(minLuminance,maxLuminance);
  }
}

void reset() {
  int clr;
  float d;
  background(0);
  counter=0;
  thinking=true;
  setall();
  for (int u = 0; u < width; u++) {
    for (int v = 0; v < height; v++) {
      lum=0;
      for (int a = 0; a<numanchors; a++) {
        weight[a]=0;
        for (int c = 0; c <numchains; c++) {
          for (int b = 0; b<numbeads; b++) {
            d=sqrt(sq(u - chains[a][c].beads[b].xi) + sq(v - chains[a][c].beads[b].yi));
            weight[a] +=chains[a][c].beads[b].sign*weightFactor*sqrt(numbeads-b)*dotSiz/d;

          }
        }
        level[a]=(int)(weight[a]*weightToLevel);
        lum+=level[a]*luminance[a];
      }
      clr=color(lum/numanchors,lum/numanchors,lum/numanchors);
      output[v*width+u] =clr;

    }
  }

  thinking=false;
}

void loop(){

  if(thinking){
    redraw=true;
    counter++;
    background(bkg);
    for(int a=0;a<numanchors;a++) {
      for(int c=0;c<numchains;c++){
        chains[a][c].drawarcs();
      }
    }
    for(int a=0;a<numanchors;a++) {
      for(int c=0;c<numchains;c++){
        chains[a][c].drawcrosslines();
      }
    }
  }
  else
  {
    if(redraw){
      for (int u = 0; u < width; u++) {
        for (int v = 0; v < height; v++) {
          pixels[v*width+u]=output[v*width+u];
        }
      }
      for(int a=0;a<numanchors;a++) {
        for(int c=0;c<numchains;c++){
          chains[a][c].drawarcs();
        }
      }
      for(int a=0;a<numanchors;a++) {
        for(int c=0;c<numchains;c++){
          chains[a][c].drawcrosslines();
        }
      }
      redraw=false;

    }
  }
}
void mousePressed(){
  if (!redraw) reset();

}

class v2d{
  float xi,yi,xf,yf,a,a0,da;
  boolean cl,annot;
  int sign=1;

  v2d(float xxi, float yyi, float xxf, float yyf, float aa0, float aa, boolean ccl, boolean ann) {
    xi=xxi;
    yi=yyi;
    xf=xxf;
    yf=yyf;
    a=aa;
    a0=aa0;
    da=random(-0.05,0.05);
    cl=ccl;
    annot=ann;
    sign=(random(100)<0)?-1:1;
  }
  void drawdots(){
    noStroke();
    fill(255,120);
    ellipse(xi,yi,4,4);
  }
  void drawlines(){
    stroke(255,20);
    line(xi,yi,xf,yf);
  }

}

class chain2d{
  v2d[] beads = new v2d[1];
  int numbeads = 1;

  chain2d(float xxi, float yyi, float xxf, float yyf,int nn) {
    numbeads=nn;
    if (numbeads<1) numbeads=1;
    beads=new v2d[numbeads];

    float dx=random(-100.0,100.0);
    float dy=random(-100.0,100.0);
    for (int i=0; i<numbeads;i++) beads[i]= new v2d(xxi+(xxf-xxi)/numbeads*i,yyi+(yyf-yyi+dy/numbeads*i)/numbeads*i,xxi+(xxf-xxi)/numbeads*(i+1),yyi+(yyf-yyi)/numbeads*(i+1),random(-0.02,0.02),random(-HALF_PI,HALF_PI),(random(100.0)<20.0),(random(100.0)<10.0));
    beads[0].sign=1;
  }
  void drawdots(){
    for(int i=0;i<numbeads;i++) beads[i].drawdots();
  }
  void drawlines(){
    for(int i=0;i<numbeads;i++) beads[i].drawlines();
  }
  void drawarcs(){
    noFill();
    for(int i=1;i<numbeads;i++){
      stroke(122+beads[i].sign*122,10);
      arc(beads[0].xi,beads[0].yi,beads[i].xi,beads[i].yi,beads[i].a0,beads[i].a,beads[i].da);
    }
  }
  void arc(float xxc, float yyc, float xxi, float yyi, float a0, float a,float da){
    float d=sqrt(sq(xxc-xxi)+sq(yyc-yyi));
    float angle;
    beginShape(LINE_STRIP);
    for(int i=0; i<17;i++){
      angle=i*a/16+a0+counter*da;
      curveVertex(xxc+(xxi-xxc)*cos(angle)+(yyi-yyc)*sin(angle),yyc+(xxi-xxc)*(-sin(angle))+(yyi-yyc)*cos(angle));
    }
    endShape();
    //ellipse(xxc,yyc,d,d);
  }

  void drawcrosslines(){
    float di, dii,r,angle,xxi,xxf,yyi,yyf;
    noFill();

    for(int i=2;i<numbeads;i++){
      if(beads[i].cl){
        stroke(122+beads[i].sign*122,50);
        di=sqrt(sq(beads[0].xi-beads[i].xi)+sq(beads[0].yi-beads[i].yi));
        dii=sqrt(sq(beads[0].xi-beads[i-1].xi)+sq(beads[0].yi-beads[i-1].yi));
        r=dii/di;
        angle=i*beads[i].a/16+beads[i].a0+counter*beads[i].da;
        xxi=beads[0].xi+(beads[i].xi-beads[0].xi)*cos(angle)+(beads[i].yi-beads[0].yi)*sin(angle);
        yyi=beads[0].yi+(beads[i].xi-beads[0].xi)*(-sin(angle))+(beads[i].yi-beads[0].yi)*cos(angle);
        xxf=beads[0].xi+r*(beads[i].xi-beads[0].xi)*cos(angle)+r*(beads[i].yi-beads[0].yi)*sin(angle);
        yyf=beads[0].yi+r*(beads[i].xi-beads[0].xi)*(-sin(angle))+r*(beads[i].yi-beads[0].yi)*cos(angle);
        line(xxi,yyi,xxf,yyf);
      }

    }
  }

  void drawtext(int a, int c){
    font=loadFont("Futura-Medium.vlw.gz");
    textFont(font,12);
    float xxi,yyi;
    fill(255);
    noStroke();
    for(int i=0;i<numbeads;i++){
      if(beads[i].annot){
        text(a + "." + c + "." + i,beads[i].xi+5,beads[i].yi+15);

      }
    }

  }
}