Function compute_bond_orders 

pub fn compute_bond_orders(
    elements: &[u8],
    positions: &[[f64; 3]],
) -> Result<BondOrderResult, String>

Compute Wiberg-like and Mayer-like bond orders from EHT.