Method new

src/simulation.rs:29–52

(seed: u64)

impl Simulation {
    pub fn new(seed: u64) -> Self {
        fastrand::seed(seed);
        let mut bodies = Vec::new();

        let n = 3;
        for _ in 0..n {
            bodies.push(rand_body());
        }

        let vel = bodies.iter().map(|b| b.vel * b.mass).sum::<Vec2>() / n as Real;
        let pos = bodies.iter().map(|b| b.pos * b.mass).sum::<Vec2>() / n as Real;
        for b in &mut bodies {
            b.vel -= vel;
            b.pos -= pos;
        }
        let r = bodies.iter().map(|b| b.pos.mag()).max_by(Real::total_cmp).unwrap();
        for b in &mut bodies {
            b.pos /= r;
        }

        Self {
            bodies,
        }
    }

    pub fn update(&mut self) {
        for i in 0..self.bodies.len() {