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

int numrr=12; // number of radii (must be even)
int numtt=8; // number of theta angles (in-plane angle)
int numpp=8; // number of phi angles (out-of-plane angle)
int pnumrr=12; // prev. number of radii (must be even)
int pnumtt=8; // prev. number of theta angles (in-plane angle)
int pnumpp=8; // prev. number of phi angles (out-of-plane angle)
float minrr=20; // minimum radius
float maxrr=240; // maximum radius
float mintt=0; // minimum theta angle
float maxtt=TWO_PI; // maximum theta angle
float minpp=PI/15.0; // minimum phi angle
float maxpp=PI-PI/7.5; // maximum phi angle
float noisrr=0.25;
float noistt=.5;
float noispp=.5;
float pnoisrr=0.25;
float pnoistt=.5;
float pnoispp=.5;
float rax=0.0;
float ray=0.0;
float raz=0.0;
float[] rr=new float[18]; // array of radii
float[] tt=new float[26]; // array of theta angles
float[] pp=new float[26]; // array of phi angles
PointSphericalCoords[][][] p = new PointSphericalCoords[18][26][26];
HScrollbar hs1,hs2,hs3,hs4,hs5,hs6,hs7,hs8,hs9;

float viewAngleX=0.0;
float viewAngleY=0.0;
float viewAngleZ=0.0;

void setup(){
  size(800,480);
  background(255);
  framerate(18);
  hs1 = new HScrollbar(660, 40, 120, 6, 1,60);
  hs2 = new HScrollbar(660, 60, 120, 6,1,40);
  hs3 = new HScrollbar(660, 80, 120, 6,1,40);
  hs4 = new HScrollbar(660, 100, 120, 6,1,30);
  hs5 = new HScrollbar(660, 120, 120, 6,1,30);
  hs6 = new HScrollbar(660, 140, 120, 6,1,30);
  hs7 = new HScrollbar(660, 160, 120, 6,1,40);
  hs8 = new HScrollbar(660, 180, 120, 6,1,60);
  hs9 = new HScrollbar(660, 200, 120, 6,1,30);
  reset();
}

void loop(){
  background(255);
  hs1.update();
  hs2.update();
  hs3.update();
  hs4.update();
  hs5.update();
  hs6.update();
  hs7.update();
  hs8.update();
  hs9.update();
  hs1.draw();
  hs2.draw();
  hs3.draw();
  hs4.draw();
  hs5.draw();
  hs6.draw();
  hs7.draw();
  hs8.draw();
  hs9.draw();
  barDraw();
  if (updateAndCheckChanges()) reset();
  
  for(int ii=0;ii<numrr;ii+=2){
    fill(135+120.0/numrr*ii,135+120.0/numrr*ii,135+120.0/numrr*ii,200.0-180.0/numrr*ii);
    stroke(0,0,0,150-150.0/numrr*ii);
    if (ii%6==2){
      for(int jj=0;jj<numtt;jj++){
        beginShape(TRIANGLE_STRIP);
        for(int kk=0;kk<numpp;kk++){
          vertex(screenX(p[ii][jj][kk].x,p[ii][jj][kk].y,p[ii][jj][kk].z),screenY(p[ii][jj][kk].x,p[ii][jj][kk].y,p[ii][jj][kk].z),-300);
          vertex(screenX(p[ii+1][jj][kk].x,p[ii+1][jj][kk].y,p[ii+1][jj][kk].z),screenY(p[ii+1][jj][kk].x,p[ii+1][jj][kk].y,p[ii+1][jj][kk].z),-300);
        }
        endShape();
      }
    }
    else if (ii%6==4){
      for(int jj=0;jj<numtt;jj++){
        beginShape(TRIANGLE_STRIP);
        for(int kk=0;kk<numpp;kk++){
          vertex(screenX(p[ii][jj][kk].x,p[ii][jj][kk].y,p[ii][jj][kk].z),screenY(p[ii][jj][kk].x,p[ii][jj][kk].y,p[ii][jj][kk].z),-300);
          vertex(screenX(p[ii+1][jj][kk].x,p[ii+1][jj][kk].y,p[ii+1][jj][kk].z),screenY(p[ii+1][jj][kk].x,p[ii+1][jj][kk].y,p[ii+1][jj][kk].z),-300);
        }
        endShape();
      }
    }
    else {
      for(int jj=0;jj<numtt;jj++){
        beginShape(TRIANGLE_STRIP);
        for(int kk=0;kk<numpp;kk++){
          vertex(screenX(p[ii][jj][kk].x,p[ii][jj][kk].y,p[ii][jj][kk].z),screenY(p[ii][jj][kk].x,p[ii][jj][kk].y,p[ii][jj][kk].z),-300);
          vertex(screenX(p[ii+1][jj][kk].x,p[ii+1][jj][kk].y,p[ii+1][jj][kk].z),screenY(p[ii+1][jj][kk].x,p[ii+1][jj][kk].y,p[ii+1][jj][kk].z),-300);
        }
        endShape();
      }
    }
  }
  translate(300,240,-200);
  rotateX(viewAngleX+=rax);
  rotateY(viewAngleY+=ray);
  rotateZ(viewAngleZ+=raz);
  for(int ii=0;ii<numrr;ii+=2){
    fill(255.0/numrr*ii,255.0/numrr*ii,255.0/numrr*ii,200.0-180.0/numrr*ii);
    stroke(0,0,0,250-250.0/numrr*ii);
    if (ii%6==2){
      for(int jj=0;jj<numtt;jj++){
        beginShape(TRIANGLE_STRIP);
        for(int kk=0;kk<numpp;kk++){
          vertex(p[ii][jj][kk].x,p[ii][jj][kk].y,p[ii][jj][kk].z);
          vertex(p[ii+1][jj][kk].x,p[ii+1][jj][kk].y,p[ii+1][jj][kk].z);
        }
        endShape();
      }
    }
    else if (ii%6==4){
      for(int jj=0;jj<numtt;jj++){
        beginShape(TRIANGLE_STRIP);
        for(int kk=0;kk<numpp;kk++){
          vertex(p[ii][jj][kk].y,p[ii][jj][kk].x,p[ii][jj][kk].z);
          vertex(p[ii+1][jj][kk].y,p[ii+1][jj][kk].x,p[ii+1][jj][kk].z);
        }
        endShape();
      }
    }
    else {
      for(int jj=0;jj<numtt;jj++){
        beginShape(TRIANGLE_STRIP);
        for(int kk=0;kk<numpp;kk++){
          vertex(p[ii][jj][kk].z,p[ii][jj][kk].y,p[ii][jj][kk].x);
          vertex(p[ii+1][jj][kk].z,p[ii+1][jj][kk].y,p[ii+1][jj][kk].x);
        }
        endShape();
      }
    }
  }
  
}

void barDraw(){
  fill(200);
  rect(660, 17, 120, 6);
  rect(660, 217, 120, 24);
  rect(660, 257, 120, 6);
  rect(660, 277, 120, 6);
  rect(660, 297, 120, 6);
  rect(660, 317, 120, 6);
  rect(660, 337, 120, 6);
  fill(210);
  rect(660, 357, 120, 6);
  fill(220);
  rect(660, 377, 120, 6);
  fill(230);
  rect(660, 397, 120, 6);
  fill(240);
  rect(660, 417, 120, 6);
  fill(250);
  rect(660, 437, 120, 6);
  fill(255);
  rect(660, 457, 120, 6);
  

}

void keyPressed()
{
  reset();
}

void reset(){
  for(int ii=1;ii<numrr+1;ii++){
    rr[ii-1]=minrr+ii*(maxrr-minrr)/float(numrr);
  }
  for(int ii=1;ii<numtt+1;ii++){
    tt[ii-1]=mintt+ii*(maxtt-mintt)/float(numtt);
  }
  for(int ii=1;ii<numpp+1;ii++){
    pp[ii-1]=minpp+ii*(maxpp-minpp)/float(numpp);
  }
  for(int ii=0;ii<numrr;ii++){
    for(int jj=0;jj<numtt;jj++){
      for(int kk=0;kk<numpp;kk++){
        p[ii][jj][kk]=new PointSphericalCoords(rr[ii],tt[jj],pp[kk]);
      }
    }
  }

}

boolean updateAndCheckChanges(){
  boolean somethingChanged=false;
  numrr=2*int(hs1.getPos()/120.0*7.0+1.0);
  numtt=int(hs2.getPos()/120.0*22.0+2.0);
  numpp=int(hs3.getPos()/120.0*22.0+2.0);
  noisrr=hs4.getPos()/120.0;
  noistt=hs5.getPos()/60.0;
  noispp=hs6.getPos()/60.0;
  somethingChanged|= (numrr!=pnumrr) || (numtt!=pnumtt) || (numpp!=pnumpp) || (noisrr!=pnoisrr) || (noistt!=pnoistt) || (noispp!=pnoispp);
  pnumrr=numrr;
  pnumtt=numtt;
  pnumpp=numpp;
  pnoisrr=noisrr;
  pnoistt=noistt;
  pnoispp=noispp;
  rax=-0.05+hs7.getPos()/1200.0;
  ray=-0.05+hs8.getPos()/1200.0;
  raz=-0.05+hs9.getPos()/1200.0;

  return somethingChanged ;
}

class PointSphericalCoords{
  float r,theta,phi,x,y,z;

  PointSphericalCoords(float rr, float tt, float pp) {
    r=rr*(1.0+random(-noisrr,noisrr));
    theta=tt+noistt*random(-PI/6.0,PI/6.0);
    phi=pp+noispp*random(-PI/6.0,PI/6.0);
    x=r*sin(phi)*cos(theta);
    y=r*sin(phi)*sin(theta);
    z=r*cos(phi);
  }

  void update(){
    x=r*sin(phi)*cos(theta);
    y=r*sin(phi)*sin(theta);
    z=r*cos(phi);
  }

  void draw(){
    point(x,y,z);
  }

}

class HScrollbar
{
  int swidth, sheight;    // width and height of bar
  int xpos, ypos;         // x and y position of bar
  float spos, newspos;    // x position of slider
  int sposMin, sposMax;   // max and min values of slider
  int loose;              // how loose/heavy
  boolean over;           // is the mouse over the slider?
  boolean locked;
  float ratio;

  HScrollbar (int xp, int yp, int sw, int sh, int l,int sp) {
    swidth = sw;
    sheight = sh;
    int widthtoheight = sw - sh;
    ratio = (float)sw / (float)widthtoheight;
    xpos = xp;
    ypos = yp-sheight/2;
    spos = xpos + sp;
    newspos = spos;
    sposMin = xpos;
    sposMax = xpos + swidth - sheight;
    loose = l;
  }

  void update() {
    if(over()) {
      over = true;
    } else {
      over = false;
    }
    if(mousePressed && over) {
      locked = true;
    }
    if(!mousePressed) {
      locked = false;
    }
    if(locked) {
      newspos = constrain(mouseX-sheight/2, sposMin, sposMax);
    }
    if(abs(newspos - spos) > 1) {
      spos = spos + (newspos-spos)/loose;
    }
  }

  int constrain(int val, int minv, int maxv) {
    return min(max(val, minv), maxv);
  }

  boolean over() {
    if(mouseX > xpos && mouseX < xpos+swidth &&
    mouseY > ypos && mouseY < ypos+sheight) {
      return true;
    } else {
      return false;
    }
  }

  void draw() {
    stroke(102);
    fill(102,102,102,50);
    rect(xpos, ypos, swidth, sheight);
    if(over || locked) {
      fill(200,200,200);
    } else {
      fill(102,102, 102);
    }
    rect(spos, ypos, sheight, sheight);
  }

  float getPos() {
    // convert spos to be values between
    // 0 and the total width of the scrollbar
    return (spos-xpos) * ratio;
  }
}