My Sketch - OpenProcessing

mySketch
Particle
tab3
xxxxxxxxxx
 
let particles = [];
​
const intendedDT = 1;
​
const substeps = 8;
​
const heating = 0.5;
​
const dt = intendedDT / (substeps);
​
const colors = [
  "#DDD",
  "#D88",
  "#88D",
  "#444",
  "#8D8",
  "#FD8",
  "#D8F",
  "#8DF"
]
​
​
function mouseClicked() {
  particles.push(new Particle(mouseX, mouseY, 6, particles.length))
}
​
function setup() {
  createCanvas(windowWidth, windowHeight);
  background(20);
  
  noStroke();
​
  for(let i = 0; i < 150; i++) {
    
    let x = random(0,width);
    let y = random(0,height);
    // let bonds = (random(0,1)**8)*8+1;
    let bonds = x > width/2 ? 1:1;
    
    particles.push(new Particle(x, y, bonds,i));
    particles[i].addForce(random(50,-50),random(50,-50)); 
  }
  
}
​
function keyPressed() { console.log(particles)}
​
function draw() {
  
  
  for(let i = 0; i < substeps; i++) {
    step();
  }
  
  background(20);
  
  for(let i = 0; i < particles.length; i++) {
    fill(255)
    
    fill(colors[particles[i].bonds-1])
    
    circle(particles[i].x, particles[i].y, 5+particles[i].bonds);
  }
}
​
function step () {
  for(let i = 0 ; i < particles.length; i++) {
    particles[i].addForce(random(heating,-heating),random(heating,-heating)); 
    
    for(let j = 0; j < particles.length; j++) {
      if(j === i) continue;
      
      let dist = sqrt((particles[i].x - particles[j].x)**2 + (particles[i].y - particles[j].y)**2);
      
      if(dist > 120) continue;
      
      let normalX = (particles[i].x - particles[j].x) / dist;
      let normalY = (particles[i].y - particles[j].y) / dist;
      
      if(dist < 3) dist = 3;
      
      
      
      if(dist < 20 && particles[i].availableBonds >= 1 ) {
        if(particles[j].availableBonds >= 1) {
          particles[i].bond(particles[j], 1);
        } else {
          particles[j].collapse(particles[i]);
        }
      } 
      
​
​
      
      if(particles[i].bondedTo(particles[j].id)) {
        
        let avgX = 0.5*(particles[i].x + particles[j].x);
        let avgY = 0.5*(particles[i].y + particles[j].y);
        
        
        particles[i].nx = avgX + normalX * (4+(particles[i].bonds+particles[j].bonds)*0.5);
        particles[i].ny = avgY + normalY * (4+(particles[i].bonds+particles[j].bonds)*0.5);
        
        particles[j].nx = avgX - normalX * (4+(particles[i].bonds+particles[j].bonds)*0.5);
        particles[j].ny = avgY - normalY * (4+(particles[i].bonds+particles[j].bonds)*0.5);
        
        
        particles[i].addForce (0.5*(particles[i].nx - particles[i].x), 0.5*(particles[i].ny - particles[i].y)) 
        particles[j].addForce (0.5*(particles[j].nx - particles[j].x), 0.5*(particles[j].ny - particles[j].y))
        
        if(particles[i].ax**2 + particles[i].ay**2 > 7) {
          particles[i].debond(particles[j])
        }
        
      } 
      
      particles[i].addForce(
        70*normalX / (dist**3), 
        70*normalY / (dist**3)
      );
​
      particles[j].addForce(
        -70*normalX / (dist**3), 
        -70*normalY / (dist**3)
      );
      
      
      
      
    }
    
    
    if(isNaN(particles[i].ax+particles[i].ay)) {
      console.log(particles[i])
      
      step = null;
      draw = null;
    }
  }
  
  for(let i = 0 ; i < particles.length; i++) {
    particles[i].bound();
    
    particles[i].updateVel();
    particles[i].updatePos();
  }
}