//competitive plotter attractor net
//by Steed
neuron [] n = new neuron [256];
plot one;
mouse f1,f2;
BImage output;
float E = 2.71828182845904523536;
void setup(){
  size (640,320);
  background(255);
  output = loadImage("output.gif");
  for (int i=0; i<n.length; i++){
    n[i] = new neuron();
    one = new plot(500,200);
    f1 = new mouse(100,100);
    f2 = new mouse(200,200);
  }
}
void loop(){
fill(0);
stroke(0);
rect(0,0,320,320);
  net2image();
 loadn();
}

void net2image (){
  for (int i =0; i < n.length; i++){
    output.pixels[i] = color (n[i].b*300,0,0);
  }
  push();
  scale(20);
  image(output,0,0);
  pop();
}

class neuron{
  float w,a,b,d;
  neuron (){
    this.a = 1.0;
    this.w = 0.0;
    this.d = 1.0;
  }
  void fire(){
    d = (d+0.05)%6.0;
    w = normalD (d, 1.0, 3.0);
    b = a * w;
  }
}//class;

void loadn(){
  int nM = 0;
  int nF1 = 0;
  int nF2 = 0;
  for (int i1 = 0; i1 < 16; i1++){
    for (int i2 = 0; i2 < 16; i2++){
      if (overRect(mouseX,mouseY,i1*20,i2*20,(i1+1)*20,(i2+1)*20)){
        n[i1+(i2*16)].fire();
        nM = i1+(i2*16);
      }
      if (overRect(f1.x,f1.y,i1*20,i2*20,(i1+1)*20,(i2+1)*20)){
        n[i1+(i2*16)].fire();
        nF1 = i1+(i2*16);
      }
      if (overRect(f2.x,f2.y,i1*20,i2*20,(i1+1)*20,(i2+1)*20)){
        n[i1+(i2*16)].fire();
        nF2 = i1+(i2*16);
      }
    }
  }
  float x2 = 0.0;
  float y2 = 0.0;
float c[] = { n[nM].b, n[nF1].b, n[nF2].b };
  sort(c);
  int col = 3;
  if (c[c.length-1] == n[nM].b){
    x2 = mouseX+320;
    y2 = mouseY;
    col = 0;
  }
  if (c[c.length-1] == n[nF1].b){
    x2 = f1.x+320;
    y2 = f1.y;
    col = 1;
  }
  if (c[c.length-1] == n[nF2].b){
    x2 = f2.x+320;
    y2 = f2.y;
    col = 2;
  }
  println(col);
  cShift(col);
  one.move(1.0,x2,y2);
  f1.move();
  f2.move();
}
class mouse{
  float x,y,th;
  mouse(float x, float y){
    this.x = x;
    this.y = y;
    this.th = 0.0;
  }
  void move(){
    th += random (1.0)-0.5;
    float m = 1.0;
    x = fwrap(x + cos(th) * m,315);
    y = fwrap(y + sin(th) * m,315);
    stroke(200);
    fill(100);
    ellipse(x-5,y-5,10,10);
 }
}
class plot{
  float x,y,th,cr,cg,cb;
  plot(float x, float y){
    this.x = x;
    this.y = y;
    cr = 93;
    cg = 93;
    cb = 76;
  }
  void move(float speed,float x2,float y2){
    float t = atan2 (y2-y,x2-x);
    if (t > th){
    th += PI/10;
    }
    if (t < th){
    th -= PI/10;
    }
    float m = speed;
    x += cos(th) * m;
    y += sin(th) * m;
    stroke(cr,cg,cb);
    println(cr+" "+cg+" "+cb);
    point(x,y);
  }
}
float normalD (float xx, float sigma, float mu){
  return  pow(E,-(sq((xx-mu))) / (2 * sq(sigma))) / (sigma*sqrt(TWO_PI));
}//where sigma is the scale of the graph & mu shifts the location of the graph along x

boolean overRect(float xo, float yo, float x1, float y1, float x2, float y2){
  if (xo >= x1 && xo <= x2 &&
  yo >= y1 && yo <= y2) {
    return true;
  } else {
    return false;
  }
}
float fwrap (float fvalu, float flimi){
  float fcar=fvalu;
if (fcar>flimi){ fcar=(fcar)%flimi; }
  if (fcar<0){
    fcar=flimi-fcar;
    fcar=(fcar)%flimi;
    fcar=flimi-fcar;
  }
  return fcar;
}//fwrap;
void cShift(int c){
if (c == 0){
one.cr = shift(one.cr, 10, 15);
one.cb = shift(one.cb, 10, 15);
one.cg = shift(one.cg, 10, 15);
}
if (c == 1){
one.cr = shift(one.cr, 200, 15);
one.cb = shift(one.cb, 10, 10);
one.cg = shift(one.cg, 10, 5);
}
if (c == 2){
one.cr = shift(one.cr, 200, 15);
one.cb = shift(one.cb, 200, 10);
one.cg = shift(one.cg, 200, 5);
}
}

float shift (float a, float b, float v){
float r = 0.0;
if (a >= b){
r = a - v;
}
if (a < b){
r = a + v;
}
return r;
}