import peasy.test.*; import peasy.org.apache.commons.math.*; import peasy.*; import peasy.org.apache.commons.math.geometry.*; import peasy.test.*; import peasy.org.apache.commons.math.*; import peasy.*; import peasy.org.apache.commons.math.geometry.*; import peasy.*; int num = 500; float RADIUS = 300; float SCALE = 500.0; int globKonexID = 0; PeasyCam cam; Neuron n[]; void setup(){ size(1440,900,P3D); n = new Neuron[num]; cam = new PeasyCam(this, 800); cam.setMinimumDistance(300); cam.setMaximumDistance(3000); cam = new PeasyCam(this, 800); cam.setMinimumDistance(300); cam.setMaximumDistance(3000); for(int i = 0;i<n.length;i++){ n[i] = new Neuron(i); } for(int i = 0;i<n.length;i++){ n[i].makeConnections(); } noSmooth(); stroke(1.0); for(int i = 0;i<n[0].k.length;i++) n[0].makeSignal(i); } void mousePressed(){ if(mouseButton==RIGHT){ for(int i = 0;i<n.length;i++){ n[i].x = random(-width/2,width/2); n[i].y = random(-width/2,width/2); n[i].z = random(-width/2,width/2); } for(int i = 0;i<n.length;i++) n[i].makeConnections(); for(int i = 0;i<n[0].sig.length;i++) n[0].makeSignal(i); } } void draw(){ background(0); for(int i = 0;i<n.length;i++){ n[i].draw(); } } ////////////////////////// >> class Neuron{ float x,y,z,radius,xx,yy,zz,val; int id; boolean imobile = false; Konexe k[] = new Konexe[0]; Signal sig[] = new Signal[0]; Neuron(int _id){ xx = x = random(-width/2,width/2); yy = y = random(-width/2,width/2); zz = z = random(-width/2,width/2); id= _id; radius = RADIUS; val = 0; if(id==1)imobile = true; } void makeConnections(){ k = new Konexe[0]; sig = new Signal[0]; for(int i = 0 ;i<n.length;i++){ if(i!=id && dist(x,y,z,n[i].x,n[i].y,n[i].z) < radius){ k = (Konexe[])expand(k,k.length+1); k[k.length-1] = new Konexe(id,i); sig = (Signal[])expand(sig,sig.length+1); sig[sig.length-1] = new Signal(k[sig.length-1]); } } } void makeSignal(int which){ int i = which; sig[i].x = xx; sig[i].y = yy; sig[i].z = zz; sig[i].running = true; } void drawConnections(){ for(int i = 0 ;i<k.length;i++){ line(xx,yy,zz,n[k[i].B].xx,n[k[i].B].yy,n[k[i].B].zz); } } void draw(){ xx += (x-xx) / 100.0; yy += (y-yy) / 100.0; zz += (z-zz) / 100.0; if(sig.length>0){ for(int i = 0;i<sig.length;i+=1){ if(sig[i].running){ if(sig[i].moving){ pushMatrix(); pushStyle(); //strokeWeight(3); stroke(255,20); //line(sig[i].x,sig[i].y,sig[i].z,sig[i].lx,sig[i].ly,sig[i].lz); translate(sig[i].x,sig[i].y,sig[i].z); box(5); popStyle(); popMatrix(); } sig[i].step(); } } } stroke(lerpColor(#FFFFFF,#33CCFF,norm(val,10,200)),200); drawConnections(); } } ////////////////////////// >> class Konexe{ int A,B,id; float weight = 1.5; Konexe(int _A,int _B){ A = _A; B = _B; id = globKonexID++; weight = random(101,200)/SCALE; } } class Signal{ Konexe base; int cyc = 0; float x,y,z,lx,ly,lz; float speed = 10.1; boolean running = false; boolean visible = true; boolean moving = false; int deadnum = 2000; int deadcount = 0; Signal(Konexe _base){ deadnum = (int)random(2,40); base = _base; lx = x = n[base.A].x; ly = y = n[base.A].y; lz = z = n[base.A].z; speed *= base.weight; } void step(){ running = true; if(abs(lx-x)>.1||abs(ly-y)>.1||abs(lz-z)>.1) moving = true; else moving = false; lx = x; ly = y; lz = z; if(!n[base.A].imobile){ x += (n[base.B].xx-x) / speed;//(speed+(dist(n[base.A].x,n[base.A].y,n[base.B].x,n[base.B].y)+1)/100.0); y += (n[base.B].yy-y) / speed;//(speed+(dist(n[base.A].x,n[base.A].y,n[base.B].x,n[base.B].y)+1)/100.0); z += (n[base.B].zz-z) / speed; } n[base.A].val+=(0-n[base.A].val)/5000.0; if(dist(x,y,z,n[base.B].xx,n[base.B].yy,n[base.B].zz)<1.0){ if(deadcount<0){ deadcount = deadnum; //deadnum += (int)random(-1,1); //println("run "+base.A+" : "+base.B); running = false; for(int i = 0; i < n[base.B].k.length;i++){ if(!n[base.B].sig[i].running && base.A!=n[base.B].sig[i].base.B){ n[base.B].makeSignal(i); n[base.B].sig[i].base.weight += (base.weight-n[base.B].sig[i].base.weight)/((dist(x,y,z,n[base.A].xx,n[base.A].yy,n[base.A].zz)+1.0)/200.0); } } //base.weight = random(1001,3000) / 1000.0; if(n[base.A].id!=1){ n[base.A].xx+=((n[base.B].x-n[base.A].x)/1.1)*noise((frameCount+n[base.A].id)/18.0);; n[base.A].yy+=((n[base.B].y-n[base.A].y)/1.1)*noise((frameCount+n[base.A].id)/17.0); n[base.A].zz+=((n[base.B].z-n[base.A].z)/1.1)*noise((frameCount+n[base.A].id)/16.0); n[base.A].xx-=((n[base.B].x-n[base.A].x)/1.1)*noise((frameCount+n[base.B].id)/18.2);; n[base.A].yy-=((n[base.B].y-n[base.A].y)/1.1)*noise((frameCount+n[base.B].id)/17.2); n[base.A].zz-=((n[base.B].z-n[base.A].z)/1.1)*noise((frameCount+n[base.B].id)/16.2); } lx = n[base.A].x; ly = n[base.A].y; lz = n[base.A].z; n[base.A].val+=(255-n[base.A].val)/10.0; }else{ deadcount--; } } } }