fn calculate_bond_angle( geometry: &Geometry, a: usize, b: usize, c: usize, ) -> f64
Calculates the current bond angle between three atoms.