\/*
\n* Copyright (c) 2012 Gwyllim Jahn
\n* http:\/\/scripts.crida.net\/gh<\/p>\n
* Developed By Gwyllim Jahn, Mark Kowalyov and ZHU Teng Fei
\n* during the 2012 RMIT nCertainties studio.
\n*
\n* This code makes extensive use of several freely distributed
\n* processing libraries \u2013 thankyou especially to toxi.<\/p>\n
* This code is free software; you can redistribute it and\/or
\n* modify it under the terms of the GNU Lesser General Public
\n* License as published by the Free Software Foundation; either
\n* version 2.1 of the License, or (at your option) any later version.
\n*
\n* http:\/\/creativecommons.org\/licenses\/LGPL\/2.1\/
\n*
\n* This code is distributed in the hope that it will be useful,
\n* but WITHOUT ANY WARRANTY; without even the implied warranty of
\n* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
\n* Lesser General Public License for more details.
\n*
\n*\/<\/p>\n
\/\/\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013strata \u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013\/\/
\nimport processing.video.*;<\/p>\n
import processing.dxf.*;<\/p>\n
import processing.opengl.*;
\nimport javax.media.opengl.*;
\nimport peasy.org.apache.commons.math.*;
\nimport peasy.*;
\nimport peasy.org.apache.commons.math.geometry.*;<\/p>\n
import toxi.math.conversion.*;
\nimport toxi.geom.*;
\nimport toxi.geom.mesh2d.*;
\nimport toxi.util.datatypes.*;
\nimport toxi.util.events.*;
\nimport toxi.geom.mesh.subdiv.*;
\nimport toxi.geom.mesh.*;
\nimport toxi.util.*;
\nimport toxi.math.noise.*;
\nimport toxi.volume.*;<\/p>\n
import toxi.math.*;
\nimport controlP5.*;
\nimport toxi.color.*;<\/p>\n
\/\/——————————Volume Params<\/p>\n
int DX=300;
\nint DY=300;
\nint DZ=300;
\nOscillatingVolume volumeV;
\nVolumetricBrush brushV;
\nVec3D SCALE=new Vec3D(DX,DY,DZ);
\nint sX = int(SCALE.x\/DX);
\nint sY = int(SCALE.y\/DY);
\nint sZ = int(SCALE.z\/DZ);<\/p>\n
\/\/——————————Agent params<\/p>\n
\/\/==============================
\n\/\/Trail length and agent speed – obviously if you want
\n\/\/to run the agents for ages, you will need longer trails.
\n\/\/agent speed also has a significant effect on the root geom.
\nint trailLength = 200;
\nfloat agentSpeed = 1;
\n\/\/ammount to confine to a particular threshold
\nfloat confineToThreshold = 0.1;
\n\/\/distance between trail points
\nint dropSpacing = 2;
\n\/\/number of agents
\nint numAgents = 500;<\/p>\n
\/\/if the agent cant find a value in its neighbouring
\n\/\/voxels that is less than this number from the ideal
\n\/\/threshold (by default 0.5) =dead. LOWER=MORE DEATH
\nfloat cutoffToDie = 0.01;
\n\/\/if the agent finds a value in its neighbouring
\n\/\/voxels that is less than this number from the ideal
\n\/\/threshold (by default 0.5), then it starts counting
\n\/\/towards branching LOWER=LESS BRANCHING
\nfloat cutoffToThinkAboutBranching = 0.00015;
\n\/\/The agent needs to have been within its branching cutoff
\n\/\/for this number of steps, then it branches.
\nint numEasyMovesBeforeBranch =10;
\n\/\/number of new agents to make.
\nint splitIntoAgents=3;<\/p>\n
\/\/==============================
\n\/\/parameters for how the agents modify the field…
\n\/\/how big the brush is that modifies the diffusion f parameter.
\nfloat rootRadius = 2;
\nint searchRad = 3;
\n\/\/density value for the agents – a negative
\n\/\/number removes material from the diffusion,
\n\/\/numbers closer to 0 mean that more agents need to fly
\n\/\/down a similar path to actually remove the material …
\nfloat density = -0.02;
\nfloat noiseScale = 2;<\/p>\n
Boolean drawAgents = true;
\nBoolean saveRaw = false;
\nMovieMaker mm;
\nBoolean filming =true;<\/p>\n
\/\/——————————Diffusion params
\nint NUM_ITERATIONS = 10;
\nControlP5 controlP5;
\nControlWindow controlWindow;<\/p>\n
PeasyCam cam;
\nWorld world;<\/p>\n
NoiseInterpolator nI = new NoiseInterpolator();<\/p>\n
void setup() {
\nsize(500,500,OPENGL);
\ncam = new PeasyCam(this,180);<\/p>\n
\/\/setup voxels
\nvolumeV=new OscillatingVolume(SCALE,DX,DY,DZ);<\/p>\n
brushV=new RoundBrush(volumeV,5);<\/p>\n
volumeV.setLandscape();<\/p>\n
setupP5();<\/p>\n
world = new World();
\nworld.spawn();
\n}<\/p>\n
void draw() {
\nbackground(255);<\/p>\n
world.update();<\/p>\n
volumeV.drawPoints();<\/p>\n
}<\/p>\n
void keyPressed(){<\/p>\n
if(key == ‘q’){
\nGsaveData(volumeV,”E:\/”+day()+”_”+hour()+”_”+minute()+”_”+second()+”voxels.raw”);
\n}<\/p>\n
if(key == ‘w’){
\nsavePts();
\n}<\/p>\n
if(key == ‘e’){
\nsaveParams();
\n}<\/p>\n
if(key == ‘ ‘){
\nmm.finish();
\n}
\n}<\/p>\n
\/\/\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013Agent \u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013\/\/
\n\/\/simple agent class to interact with the data
\n\/\/field and other agents<\/p>\n
class Agent {<\/p>\n
Vec3D loc,vel,origLoc;
\nArrayList trail;
\nfloat gVal;
\nfloat dropNum;
\nboolean alive;
\nint id;
\nfloat aX,aY,aZ,lastBest;
\nint easy;<\/p>\n
Agent(Vec3D _loc, Vec3D _vel,int _id){
\nloc = _loc;
\nvel = _vel;
\ntrail = new ArrayList();
\ngVal = 255;
\ndropNum = 0;
\nalive = true;
\nid = _id;
\neasy = 0;
\nlastBest = 0;
\n}<\/p>\n
void run(){
\nif (alive){
\ngetGridPos();
\npushPullToData(true,0.5);<\/p>\n
\/\/branch();<\/p>\n
volBrushLocation();<\/p>\n
update();
\n}
\nrender();
\n}<\/p>\n
\/\/————————————————————————————-<\/p>\n
\/\/Behaviour for avoiding the worst (highest greyscale value) data within a proximity,
\n\/\/and interpolating the agent velocity away that point.<\/p>\n
\/\/————————————————————————————-<\/p>\n
void getGridPos(){
\naX = constrain(loc.x,0,SCALE.x);
\naX = map(aX,0,SCALE.x,0,DX);
\naY = constrain(loc.y,0,SCALE.y);
\naY = map(aY,0,SCALE.y,0,DY);
\naZ = constrain(loc.z,0,SCALE.z);
\naZ = map(aZ,0,SCALE.z,0,DZ);
\nif(loc.x>SCALE.x || loc.y >SCALE.y || loc.z>SCALE.z)vel.invert();
\nif(loc.x }<\/p>\n
void pushPullToData(Boolean push, float t){<\/p>\n
float bestThreshDiff = 100;
\nVec3D bestThreshPt = new Vec3D(0,0,0);<\/p>\n
\/\/iterate over neighbouring grid pos
\nfor (int i =-searchRad;i for (int j =-searchRad;j for (int k =-searchRad; k=0&&(aX+i)=0&&(aY+j)=0&&(aZ+k)0.0001 &&val float threshVal = abs(val-t);<\/p>\n
if(threshVal cutoffToDie){
\nalive=false;
\nprint(“dead”);
\n}<\/p>\n
\/\/if in a good posy, see if we should branch
\nif(bestThreshDiff easy++;
\n}else{
\neasy = 0;
\n}
\n}<\/p>\n
Boolean moveToThreshold(int searchRad){<\/p>\n
getGridPos();<\/p>\n
for (int i =-searchRad;ifor (int j =-searchRad;jfor (int k =-searchRad;k=0&&(aX+i)=0&&(aY+j)=0&&(aZ+k)=1&&valnumEasyMovesBeforeBranch){
\neasy =0;
\nfor (int i = 1;i0) brushV.drawAtGridPos(aX+5,aY,aZ,-density*2);<\/p>\n
\/\/ if(aZ-5>0) brushV.drawAtGridPos(aX,aY,aZ-5,-density*2);
\n\/\/ brushV.drawAtGridPos(aX,aY,aZ,density);
\n}<\/p>\n
\/\/————————————————————————————-<\/p>\n
\/\/Functions for managing trails and rendering the agent.<\/p>\n
\/\/————————————————————————————-<\/p>\n
void update(){
\n\/\/move the agent randomly
\nvel.normalizeTo(agentSpeed);
\n\/\/vel.interpolateToSelf(new Vec3D (0,0,1),vert);
\nloc.addSelf(vel);<\/p>\n
\/\/update the trail geom every 5 moves
\nif(dropNum%dropSpacing==0){
\nif(trail.size() trail.add(new Vec3D(loc.x, loc.y, loc.z));
\n}else{
\ntrail.remove(0);
\n}
\ndropNum=0;
\n}
\ndropNum++;<\/p>\n
}<\/p>\n
void render(){
\n\/\/render agent trail as lines
\nVec3D ls = null;
\nfor (Vec3D l : trail){
\nif(ls !=null){
\nstroke(0);
\nline (l.x,l.y,l.z,ls.x,ls.y,ls.z);
\n}
\nls = l;
\n}
\n}<\/p>\n
}
\n\/\/\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013Controllers \u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013\/\/
\nvoid setupP5 (){
\ncontrolP5 = new ControlP5(this);
\ncontrolP5.setAutoDraw(false);
\ncontrolWindow = controlP5.addControlWindow(“controlP5window”,100,100,400,300);
\ncontrolWindow.hideCoordinates();<\/p>\n
controlWindow.setTitle(“Spring GUI”);
\n\/\/\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013Oscillating Volumes \u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013\/\/
\nclass OscillatingVolume extends VolumetricSpaceArray {<\/p>\n
int z;<\/p>\n
public OscillatingVolume (Vec3D _S, int _DX, int _DY,int _DZ) {
\nsuper(_S,_DX,_DY,_DZ);
\nz = 0;
\n}<\/p>\n
\/\/—————————<\/p>\n
void setLandscape(){<\/p>\n
for (int i=0;ifor (int j =0;jfor (int k=0;k<\/p>\n
\/\/create the noise field – dividing k steps the field.
\n\/\/used to emphasise the difference generated by the agents
\n\/\/setVoxelAt(i,j,k,nI.getValType(i,j,k\/4));<\/p>\n
\/\/use this function to make the entire field a particular type
\nsetVoxelAt(i,j,k,nI.getValType(i,j,k\/4,0));
\n}
\n}
\n}
\n}<\/p>\n
void drawPoints() {
\nint r=4;
\nfor(int i=0;i for(int j=0;j for(int k=0;k=90&&v stroke(100);
\nline(i,j,k,i,j,k+1);
\n}
\n}
\n}
\n}<\/p>\n
}<\/p>\n
}
\n\/\/\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013Voxel Functions \u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013\/\/
\npublic void GsaveData(OscillatingVolume v,String fn) {
\nprint(“saving volume data…”);
\nint c=0;
\nint tot = v.getData().length;
\ntry {
\nBufferedOutputStream ds = new BufferedOutputStream(new FileOutputStream(fn));
\n\/\/ ds.writeInt(volumeData.length);
\nfor (float element : v.getData()) {
\nif(element ds.write((int)(element*100));<\/p>\n
}
\nds.flush();
\nds.close();
\n} catch (IOException e) {
\ne.printStackTrace();
\n}
\n}
\n\/\/\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013World\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013\/\/
\n\/\/Class to manage addition and subtraction of agents from
\n\/\/the environment<\/p>\n
class World {<\/p>\n
ArrayListagentPop;
\nArrayList newPop;<\/p>\n
World(){<\/p>\n
agentPop = new ArrayList();
\nnewPop = new ArrayList();<\/p>\n
}<\/p>\n
\/\/————————–main function<\/p>\n
void update(){
\nnewPop = new ArrayList();<\/p>\n
runAgents();<\/p>\n
if(newPop.size()>0) addAgents();
\n}<\/p>\n
void runAgents(){
\nint c = 0;
\nfor (Agent a :agentPop){
\na.run();
\n}
\n}<\/p>\n
void addAgents(){
\nagentPop.addAll(newPop);
\n}<\/p>\n
void spawn(){
\nnewPop = new ArrayList();
\nagentPop = new ArrayList();
\nfor (int i=0;i \/\/ Vec3D loc = new Vec3D(random(0,gW*gRes),random(0,gH*gRes),0);
\nVec3D loc = new Vec3D(random(0,DX),random(0,DY),SCALE.z);
\nVec3D vel = new Vec3D(0,0,-3);
\nAgent a = new Agent(loc,vel,i);
\na.origLoc = loc.copy();
\nagentPop.add(a);
\n}
\n}<\/p>\n
void spawnOnThreshold(){
\nnewPop = new ArrayList();
\nagentPop = new ArrayList();
\nfor (int i=0;i50) searching = false;
\n}<\/p>\n
}
\n}
\n}
\n\/\/\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013Noise Interpolator \u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013\/\/
\nclass NoiseInterpolator {<\/p>\n
int aLimit,bLimit,cLimit,dLimit,eLimit;
\nfloat aScale,bScale,cScale,dScale,eScale;
\nfloat aContrast, bContrast, cContrast, dContrast,eContrast;
\nfloat aCompression,bCompression,cCompression,dCompression,eCompression;<\/p>\n
NoiseInterpolator(){<\/p>\n
aLimit = 0;
\naScale = 1.75;
\naContrast = 0.2;
\naCompression = 1;<\/p>\n
bLimit = (int)SCALE.z\/4;
\nbScale = 4;
\nbContrast = 0;
\nbCompression = 3;<\/p>\n
cLimit = (int)SCALE.z\/2;
\ncScale = 7;
\ncContrast = 0.35;
\ncCompression = 8;<\/p>\n
dLimit = (int)(SCALE.z\/4)*3;
\ndScale = 4;
\ndContrast = 0;
\ndCompression = 1;<\/p>\n
eLimit = (int)SCALE.z;
\neScale = 10;
\neContrast = 0;
\neCompression = 1;
\n}<\/p>\n
float getVal(int x, int y, int z){<\/p>\n
int type = floor(((float)z\/DZ)*4);
\nfloat f = ((float)z%(DZ\/4));
\nf = (f\/(DZ\/4));<\/p>\n
float ns = noiseScale;
\nfloat nContrast = 0;
\nfloat nCompression =1;<\/p>\n
if(type == 0){
\n\/\/between a and b<\/p>\n
ns = aScale + ( bScale – aScale ) * f;
\nnContrast = aContrast + ( bContrast – aContrast ) * f;
\nnCompression = aCompression + ( bCompression – aCompression ) * f;<\/p>\n
}else if (type ==1){
\n\/\/between b and c
\nns = bScale + ( cScale – bScale ) * f;
\nnContrast = bContrast + ( cContrast – bContrast ) * f;
\nnCompression = bCompression + ( cCompression – bCompression ) * f;
\n}
\nelse if (type ==2){
\n\/\/between c and d
\nns = cScale + ( dScale – cScale ) * f;
\nnContrast = cContrast + ( dContrast – cContrast ) * f;
\nnCompression = cCompression + ( dCompression – cCompression ) * f;
\n}
\nelse if (type ==3){
\n\/\/between d and …
\nns = dScale + ( eScale – dScale ) * f;
\nnContrast = dContrast + ( eContrast – dContrast ) * f;
\nnCompression = dCompression + ( eCompression – dCompression ) * f;
\n}
\nns = ns\/100;
\nfloat n = noise(x*ns,y*ns,z*ns*nCompression);
\nn = (map(n,nContrast,1-nContrast,0,1));<\/p>\n
return n;
\n}<\/p>\n
float getValType(int x, int y, int z,int type){<\/p>\n
float ns = noiseScale;
\nfloat nContrast = 0;
\nfloat nCompression =1;<\/p>\n
if(type == 0){
\n\/\/between a and b<\/p>\n
ns = aScale;
\nnContrast = aContrast;
\nnCompression = aCompression;<\/p>\n
}else if (type ==1){
\n\/\/between b and c
\nns = bScale;
\nnContrast = bContrast;
\nnCompression = bCompression;
\n}
\nelse if (type ==2){
\n\/\/between c and d
\nns = cScale;
\nnContrast = cContrast;
\nnCompression = cCompression;
\n}
\nelse if (type ==3){
\n\/\/between d and …
\nns = dScale;
\nnContrast = dContrast ;
\nnCompression = dCompression;
\n}
\nns = ns\/100;
\nfloat n = noise(x*ns,y*ns,z*ns*nCompression);
\nn = (map(n,nContrast,1-nContrast,0,1));<\/p>\n
return n;
\n}<\/p>\n
}
\n\/\/\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013saving \u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2014\u2013\/\/
\nvoid savePts(){
\nString[] lines = new String[world.agentPop.size()];
\nint c=0;
\nfor (int i =0;i Agent a = (Agent) world.agentPop.get(i);
\nlines[c] = “”;<\/p>\n
for(Vec3D t:a.trail){
\nlines[c] =lines[c]+ (t.x +”,” + t.y + “,” + t.z + “\/”);
\n}
\nc++;
\n}<\/p>\n
saveStrings(“G:\/ncertainties\/export\/”+day()+”_”+hour()+”_”+minute()+”_”+second()+”lines.txt”, lines);
\n}<\/p>\n
void saveParams(){
\nString[] lines = new String[23];
\nlines[0] = “agentSpeed=”+agentSpeed;
\nlines[1] = “trailLength=”+trailLength;
\nlines[2] = “confineToThreshold=”+confineToThreshold;
\nlines[3] = “dropSpacing=”+dropSpacing;
\nlines[4] = “numAgents=”+numAgents;
\nlines[5] = “cutoffToDie=”+cutoffToDie;
\nlines[6] = “cutoffToThinkAboutBranching=”+cutoffToThinkAboutBranching;
\nlines[7] = “numEasyMovesBeforeBranch=”+numEasyMovesBeforeBranch;
\nlines[8] = “splitIntoAgents=”+splitIntoAgents;
\nlines[10] = “rootRadius=”+rootRadius;
\nlines[11] = “density=”+density;
\nlines[20] = “DX=”+DX;
\nlines[21] = “DY=”+DY;
\nlines[22] = “DZ=”+DZ;
\nsaveStrings(“G:\/ncertainties\/export\/”+day()+”_”+hour()+”_”+minute()+”_”+second()+”paramaters.txt”, lines);
\n}<\/p>\n","protected":false},"excerpt":{"rendered":"
\/* * Copyright (c) 2012 Gwyllim Jahn * http:\/\/scripts.crida.net\/gh * Developed By Gwyllim Jahn, Mark Kowalyov and ZHU Teng Fei * during the 2012 RMIT nCertainties studio. * * This code makes extensive use of several freely distributed * processing libraries \u2013 thankyou especially to toxi. * This code is free software; you can redistribute […]<\/p>\n","protected":false},"author":1,"featured_media":0,"parent":20,"menu_order":0,"comment_status":"closed","ping_status":"open","template":"","meta":{"footnotes":""},"class_list":["post-585","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/www.new-territories.com\/blog\/rmit.uts\/wp-json\/wp\/v2\/pages\/585","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.new-territories.com\/blog\/rmit.uts\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/www.new-territories.com\/blog\/rmit.uts\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/www.new-territories.com\/blog\/rmit.uts\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.new-territories.com\/blog\/rmit.uts\/wp-json\/wp\/v2\/comments?post=585"}],"version-history":[{"count":3,"href":"https:\/\/www.new-territories.com\/blog\/rmit.uts\/wp-json\/wp\/v2\/pages\/585\/revisions"}],"predecessor-version":[{"id":1012,"href":"https:\/\/www.new-territories.com\/blog\/rmit.uts\/wp-json\/wp\/v2\/pages\/585\/revisions\/1012"}],"up":[{"embeddable":true,"href":"https:\/\/www.new-territories.com\/blog\/rmit.uts\/wp-json\/wp\/v2\/pages\/20"}],"wp:attachment":[{"href":"https:\/\/www.new-territories.com\/blog\/rmit.uts\/wp-json\/wp\/v2\/media?parent=585"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}