/* 
   "Stand By":- pre beta - by Alfredo Calosci, based on:
   Structure 3 (work in progress) 
   
   A surface filled with one hundred medium to small sized circles. 
   Each circle has a different size and direction, but moves at the same slow rate. 
   Display: 
   A. The instantaneous intersections of the circles 
   B. The aggregate intersections of the circles 
 
   Implemented by Casey Reas <http://groupc.net> 
   Uses circle intersection code from William Ngan <http://metaphorical.net> 
   Processing v.1.0.9  <http://processing.org> 
 
*/ 
 
 
Circle circleA, circleB, circleC; 
 
void setup() 
{ 
  size( 768, 628 ); 
  frameRate( 30 ); 
 
  circleA = new Circle( 150, 300, 90 ); 
  circleB = new Circle( 200, 300, 90 );
 circleC = new Circle( 300, 400, 120 ); 
 
  ellipseMode(CENTER_DIAMETER); 
  //rectMode(CENTER_DIAMETER); 
  noStroke();
  smooth();
 background(255); 
} 
 
 
void draw() { 
  
  circleA.update(); 
  circleB.update();
  circleC.update();
   
  //circleA.display();
  //circleB.display(); 
  //circleC.display(); 
 
 
  intersect( circleA, circleB ); 
  intersect( circleA, circleC ); 
  intersect( circleB, circleC ); 
  
  // tapa
  fill(128);
  noStroke();
  rect(0,0,(72/4)*16,(16/4)*16);
  rect((120/4)*16,0,(72/4)*16,(16/4)*16);
  
  rect(0,(16/4)*16,(36/4)*16,(16/4)*16);
  rect((156/4)*16,(16/4)*16,(36/4)*16,(16/4)*16);
  
} 
 
 
class Circle 
{ 
  float x, y, r, r2; 
  float xspeed, xdir; 
  float yspeed, ydir; 
 
  Circle( float px, float py, float pr ) { 
    x = px; 
    y = py; 
    r = pr; 
    r2 = r*r; 
    xspeed = random(-2.0, 2.0); 
    yspeed = random(-2.0, 2.0); 
    xdir = 1; 
    ydir = -1; 
  } 
  
  void update() { 
  
    x += xspeed * xdir; 
    if(x > width - r || x < r) { 
      xdir = xdir * -1; 
    } 
    
    y += yspeed * ydir; 
    //if(y > height - r || y < r) { 
    if(y > height  || y < r) { 
      ydir = ydir * -1; 
    } 
    
  } 
  
  void display(){
    stroke(0, 0, 0, 26); 
    noFill();
    ellipse(x,y,r*2,r*2);
  }
} 
 
void intersect( Circle cA, Circle cB ) { 
 
  float dx = cA.x - cB.x; 
  float dy = cA.y - cB.y; 
  float d2 = dx*dx + dy*dy; 
  float d = sqrt( d2 ); 
 
  if ( d>cA.r+cB.r || d<abs(cA.r-cB.r) ) return; // no solution 
  
  float a = (cA.r2 - cB.r2 + d2) / (2*d); 
  float h = sqrt( cA.r2 - a*a ); 
  float x2 = cA.x + a*(cB.x - cA.x)/d; 
  float y2 = cA.y + a*(cB.y - cA.y)/d; 
 
  stroke(0, 0, 0, 26); 
 
  float paX = x2 + h*(cB.y - cA.y)/d; 
  float paY = y2 - h*(cB.x - cA.x)/d; 
  float pbX = x2 - h*(cB.y - cA.y)/d; 
  float pbY = y2 + h*(cB.x - cA.x)/d; 
 
 
  //line(paX, paY, pbX, pbY); 
  dotLine(paX, paY, pbX, pbY);
} 

void dotLine(float x1, float y1, float x2, float y2){
  if(x2== x1) return;
  float M = (y2-y1)/(x2-x1);
  float Q = ((x2*y1)-(x1*y2))/(x2-x1);
  
  if(y2==y1){
    // y = y1 
  }
  
  // x normalizadas
  int x1n = floor(x1/16)*16;
  int x2n = floor(x2/16)*16;
  
  for(int n = x1n; n<= x2n; n+=16){
    int yy = int(M*n + Q);
    int yn = floor(yy/16)*16;
    fill(0,128,0,12);
    noStroke();
    //rect(n,yn,16,16);
    ellipse(n+8,yn+8,12,12);
  }
  
}

void mousePressed(){
  //saveFrame("standBy##.png");
}