Method update

src/simulation.rs:54–75

(&mut self)

    }

    pub fn update(&mut self) {
        for i in 0..self.bodies.len() {
            let p1 = self.bodies[i].pos;
            let m1 = self.bodies[i].mass;
            for j in (i+1)..self.bodies.len() {
                let p2 = self.bodies[j].pos;
                let m2 = self.bodies[j].mass;

                let r = p2 - p1;
                let mag_sq = r.x * r.x + r.y * r.y;
                let mag = mag_sq.sqrt();
                let tmp = r / (mag_sq.max(MIN) * mag);

                self.bodies[i].acc += m2 * tmp;
                self.bodies[j].acc -= m1 * tmp;
            }
        }

        for i in 0..self.bodies.len() {
            self.bodies[i].update(DT);
        }
    }
}