int NN = 762;
int N = NN;
PSound flush;
long len;
long time;
boolean doFlush;
boolean doIntro = true;
float[] px = new float[N];
float[] py = new float[N];
float[] pz = new float[N];

float[] dx = new float[N];
float[] dy = new float[N];
float[] dz = new float[N];

float ss = 0.1;
float s = (ss *600)/N;
void setup() {
  framerate(20);
  size(300,200,P3D);

  for(int i = 0; i<N; i++) {
    px[i] = random(200)-100;
    pz[i] = random(200)-100;
    py[i] = random(60)-30;

    dx[i] = .0*pz[i];
    dz[i] = -.0*px[i];
    dy[i] = 0;
  }

  flush = loadSound("flush.wav");
  len = (long)(1000*flush.duration());

  //sphereDetail(3);
}

void draw() {
  perspective(PI/4, 3./2, -100, 100);

  translate(150, 100, 0);
  rotateX(alpha);
  rotateY(beta);

  step();
  background(0);
  stroke(255);
  beginShape(POINTS);
  for(int i = 0; i<N; i++) {
    vertex(px[i],py[i],pz[i]);
    //    translate(px[i],py[i],pz[i]);
    //sphere(.5);
    //translate(-px[i],-py[i],-pz[i]);
  }
  endShape();
}

void step() {


  for(int i = 0; i<N; i++) {
    for(int j = i+1; j<N; j++) {
      float x = px[i] -px[j];
      float y = px[i] - py[j];
      float z = px[i] - pz[j];
      float d = -s/(x*x +y*y + z*z);

      dx[i] += x*d;
      dx[j] += -x*d;
      dy[i] += y*d;
      dy[j] += -y*d;
      dz[i] += z*d;
      dz[j] += -z*d;
    }
  }
  if(doIntro) {

    if(inCnt>300) {
      inCnt = 0;
      doIntro = false;
      N = NN;
    } 
    else {
      intro(inCnt);
      inCnt++;
    }
    return;
  } 
  if(doFlush) {
    for(int i = 0; i<NN; i++) {

      float d = sqrt(px[i]*px[i] +py[i]*py[i] + pz[i]*pz[i] );
      d *=d;

      dx[i] = 0.9*dx[i];
      dy[i] = 0.9*dy[i];
      dz[i] = 0.9*dz[i];

      px[i] += dx[i]- .5*px[i]/d;
      py[i] += dy[i]- .5*py[i]/d;
      pz[i] += dz[i]- .5*pz[i]/d;

      if(System.currentTimeMillis()-time >len) {
        doFlush = false;
        flush.stop();
        N = NN;
        inCnt = 0;
        doIntro = true;
      }

    }
    N-= 5;
  } 
  else
    for(int i = 0; i<N; i++) {

      float d = sqrt(px[i]*px[i] +py[i]*py[i] + pz[i]*pz[i] );

      d *=d;
      dx[i] = 0.98*dx[i] - ss*px[i]/d;
      dy[i] = 0.98*dy[i] - ss*py[i]/d;
      dz[i] = 0.98*dz[i] - ss*pz[i]/d;

      px[i] += dx[i];
      py[i] += dy[i];
      pz[i] += dz[i];


    }
}

int oldX;
int oldY;
float alpha= 0*PI/6.;
;
float beta = 0*PI/3.;
void mousePressed() {
  oldX = mouseX;
  oldY = mouseY;
  if(//mouseButton == RIGHT & 
  !doIntro) {
    flush.play();
    time = System.currentTimeMillis();
    doFlush = true;
  }
}
//void mouseDragged() {
//  alpha += (mouseX-oldX)/200.;
//  beta += (mouseY-oldY)/200.;
//  oldX = mouseX;
//  oldY = mouseY;
//}

int inCnt=0;
void intro(int i) {
  if(i<150)   N = 7;
  int x = i<150?150-i:0;
  int y = 80;
  px[0] = 10 +x;
  py[0] = 14;
  pz[0] = 0;

  px[1] = 14 +x;
  py[1] = 42;
  pz[1] = 0;

  px[2] = 48 +x;
  py[2] = 42;
  pz[2] = 0;

  px[3] = 52 +x;
  py[3] = 17;
  pz[3] = 0;

  px[4] = 64 +x;
  py[4] = 14;
  pz[4] = 0;

  px[5] = 78 +x;
  py[5] = 8;
  pz[5] = 0;

  px[6] = 98 +x;
  py[6] = 14;
  pz[6] = 0;


  if(i>=150) {
    i -= 150;
    if(i>150)  i = 150;

    float si = sin(PI*i/300.);
    float cs = cos(PI*i/300.);
    for(int j = 0; j<7;j++) {
      float xx = cs*px[j] + si*py[j];
      py[j] = -si*px[j] + cs*py[j];
      px[j] = xx; 
    }

    for(int j = N; j<N+5;j++) {
      px[j] = px[0] ;
      py[j] = py[0];
      pz[j] = 0;
      dx[j] = random(2.)-1.9;
      dy[j] = 2.5;
      dz[j] = random(2.)-1.;

    }
    N = N + 5;
    for(int j = 7; j<N;j++) { 
      px[j] += dx[j];
      py[j] += dy[j];
      pz[j] += dz[j];
      dy[j] += .1;
      if(py[j]>100) {
        px[j] = px[0] ;
        py[j] = py[0];
        pz[j] = 0;
        dy[j] = 0;
      }
    }
  }
}