quickdudley/eocalon - Change DMOMPCJ2SNZ7SEAVPROANSOHQDGVL2OKD223EC24ZF5GJJ53SYTQC

Orbit picker: calcuate accelleration on a per-body basis since I'm now planning to use the DSKD algorithm

Created by quickdudley on June 16, 2021

DMOMPCJ2SNZ7SEAVPROANSOHQDGVL2OKD223EC24ZF5GJJ53SYTQC

Dependencies

In discussions

#1 Calendar

Change contents

Insertion in orbitpick/src/vec3.rs at line 101 [3.10]

[3.1984]

impl std::iter::Sum<V3> for V3 {
    fn sum<I: Iterator<Item = V3>>(iter: I) -> V3 {
        iter.fold(
            V3 {
                x: 0.,
                y: 0.,
                z: 0.,
            },
            |a, b| a + b,
        )
    }
}

Replacement in orbitpick/src/vec3.rs at line 116 [3.10]

B:BD[3.2030] → [2.59:128]

        (self.x * self.x + self.y * self.y + self.z * self.z).sqrt()

[3.2030]

[3.2093]

        self.squared_magnitude().sqrt()
    }
    pub fn squared_magnitude(self) -> f64 {
        self.x * self.x + self.y * self.y + self.z * self.z

Insertion in orbitpick/src/nbody.rs at line 21 [2.140]

[2.557]

    }
    pub fn accelleration<N: Iterator<Item = Self>>(&self, others: N) -> V3 {
        others
            .filter(|other| other.position != self.position)
            .map(|other| {
                let d = other.position - self.position;
                d.normalize() * other.mass / d.squared_magnitude()
            })
            .sum::<V3>()
            * G

Deletion in orbitpick/src/nbody.rs at line 47 [2.140]

B:BD[2.923] → [2.923:1534]

}
pub fn forces(state: Vec<Body>) -> Vec<V3> {
    let mut result = vec![
        V3 {
            x: 0.,
            y: 0.,
            z: 0.
        };
        state.len()
    ];
    for (ix, xb) in state.iter().enumerate() {
        for (iy, yb) in state.iter().enumerate() {
            if ix != iy {
                let d = yb.position - xb.position;
                let r2 = d.x * d.x + d.y * d.y + d.z * d.z;
                let f = xb.mass * yb.mass / r2;
                result[ix] += d.normalize() * f;
            }
        }
    }
    for b in result.iter_mut() {
        *b *= G;
    }
    result