Function compute_dos_parallel 

pub fn compute_dos_parallel(
    orbital_energies: &[f64],
    sigma: f64,
    e_min: f64,
    e_max: f64,
    n_points: usize,
) -> DosResult

Compute total DOS using rayon over the energy grid.