import javax.media.opengl.*;
import processing.opengl.*;

PFont font;
PImage spot;
SphericalDistanceCalculator SDC;
PointOnUnitSphere pous[];


int nPous;
int maxNPous;
float mx = 0;
float my = 0;

void setup(){
  size(400,440,OPENGL);
  hint(ENABLE_OPENGL_4X_SMOOTH);

  spot = loadImage("spot32inv.jpg"); 
  font = loadFont("SnellRoundhand-Bold-16.vlw"); 

  nPous = 1;
  maxNPous = 256;
  pous = new PointOnUnitSphere[maxNPous];
  for (int i=0; i<maxNPous; i++){
    float lat = acos(random(-1,1)) - HALF_PI;
    float lon = random(0, TWO_PI); 
    pous[i] = new PointOnUnitSphere(i);
    pous[i].setPositionFromCorrectLatLon(lat, lon);
  }

  SDC = new SphericalDistanceCalculator();
}




long lastFrameTime = 0;
long lastAddTime = 0;
float mspf = 16.6666;
//--------------------------------------------
void draw(){
  
  long now = millis();
  updateEu();

  background(0,15,45);
  noStroke();

  fill(26,46,106);
  textFont(font);
  textAlign(CENTER);
  text("Wishing you interesting new connections in 2009\nFrom Golan, Andrea and Zenon", width/2,height-44); 


  // Hold the framerate constant by adjusting the number of points we see.
  float A = 0.95;
  float B = 1.0-A;
  mspf = A*mspf + B*(now - lastFrameTime);
  lastFrameTime = now;
  if (mspf > 60){
    nPous = max(20, nPous-4);
  } 
  if (frameCount % 10 == 0){
    if ((mspf < 16.66666) || (millis()<5000)) {
      nPous = min(maxNPous, nPous+1);
      lastAddTime = now;
    } 
    else if (mspf > 17.66666){
      if (millis()>5000){
        nPous = max(20, nPous-1);
      }
    }
  }


  // Choose the right blend function.
  fill(255);
  int a,b;
  PGraphicsOpenGL pgl = (PGraphicsOpenGL) g;  // g may change
  GL gl = pgl.beginGL(); 
  gl.glDisable(gl.GL_DEPTH_TEST);
  gl.glEnable(gl.GL_BLEND);
  int sfact[] = {
    gl.GL_ZERO, gl.GL_ONE, gl.GL_SRC_COLOR, gl.GL_ONE_MINUS_SRC_COLOR, gl.GL_DST_COLOR, 
    gl.GL_ONE_MINUS_DST_COLOR, gl.GL_SRC_ALPHA, gl.GL_ONE_MINUS_SRC_ALPHA, gl.GL_DST_ALPHA, 
    gl.GL_ONE_MINUS_DST_ALPHA, gl.GL_CONSTANT_COLOR, gl.GL_ONE_MINUS_CONSTANT_COLOR,
    gl.GL_CONSTANT_ALPHA, gl.GL_ONE_MINUS_CONSTANT_ALPHA, gl.GL_SRC_ALPHA_SATURATE                                                                                };
  int dfact[] = {
    gl.GL_ZERO, gl.GL_ONE, gl.GL_SRC_COLOR, gl.GL_ONE_MINUS_SRC_COLOR, gl.GL_DST_COLOR, 
    gl.GL_ONE_MINUS_DST_COLOR, gl.GL_SRC_ALPHA, gl.GL_ONE_MINUS_SRC_ALPHA, gl.GL_DST_ALPHA, 
    gl.GL_ONE_MINUS_DST_ALPHA, gl.GL_CONSTANT_COLOR, gl.GL_ONE_MINUS_CONSTANT_COLOR, 
    gl.GL_CONSTANT_ALPHA, gl.GL_ONE_MINUS_CONSTANT_ALPHA                                                                                };
  a = 2;
  b = 8;
  gl.glBlendFunc (sfact[a], dfact[b]);
 
 
  // Rotate and draw the points.
  float R = 300;
  mx = 0.99*mx + 0.01*(mouseX-width/2);
  my = 0.99*my + 0.01*(mouseY);
  
  pushMatrix();
  translate(width/1, 400, 0); 
  rotateY(millis()/5000.0);
  rotateY(mx/50.0);
  rotateX(HALF_PI); 
  rotateX(0-my/50.0); 

  for (int i=0; i<nPous; i++){
    float d = 140 + 30*noise(i+millis()/(i+147.0));
    float r = 120 + 60*noise((i+millis())/1003.0);
    float g = 150 + 50*noise(i+millis()/ 777.0);

    tint(r,g,255);
    pous[i].draw(R,d);

    tint(220,255,255);
    pous[i].draw(R+3,d*0.50);

  }

  // Draw the connections between the points.
  a = 4;//(mouseX/10)%sfact.length; //4
  b = 5;//(mouseY/10)%dfact.length; //14;//(millis()/600)%dfact.length; //5
  gl.glBlendFunc (sfact[a], dfact[b]);
  smooth();
  stroke(255,255,255,10);
  float M = 0.45;
  for (int i=0; i<nPous; i++){
    PointOnUnitSphere S = pous[i];
    for (int j=0; j<i; j++){
      PointOnUnitSphere F = pous[j];
      float distance = SDC.getSphericalAngleArccos (S, F);
      if (distance < M){
        SDC.drawArcBetweenTwoPointsOnUnitSphere (S, F, R, distance);
      }
    }
  }


  popMatrix();

  // Draw a highlight that sweeps across the text
  tint(180,220,255);
  a = 14;
  b = 13;
  gl.glBlendFunc (sfact[a], dfact[b]);
  float t = ((millis()-1500)/3)%6000;
  image(spot,t-200, 810-50,400,100);
  
  pgl.endGL();
}




//--------------------------------------------
// Update the forces on the points, 
// so that they mutually repel or attract.
//
void updateEu(){

  float sub = 0.47;
  float C =  0.00002 * (noise(millis()/200.0)-sub);
  if (mousePressed){
    C = -0.00001;
  }

  for (int i=0; i<nPous; i++){
    PointOnUnitSphere S = pous[i];
    for (int j=0; j<i; j++){
      PointOnUnitSphere F = pous[j];

      float dx = S.x - F.x;
      float dy = S.y - F.y; 
      float dz = S.z - F.z;
      float dh2 = (dx*dx + dy*dy + dz*dz); 

      float fx = C * dx/dh2;
      float fy = C * dy/dh2;
      float fz = C * dz/dh2;

      S.accumulateEuclideanForce ( fx, fy, fz);
      F.accumulateEuclideanForce (-fx,-fy,-fz);

    }
  }


  for (int i=0; i<nPous; i++){
    pous[i].updateEuclideanForces();
  }
}