Function first_brillouin_1d 

pub fn first_brillouin_1d(k: f64, b: f64) -> Result<f64>

Reduce a scalar wavevector to (-b/2, b/2], for positive reciprocal pitch b.

Upstream: treams.misc.firstbrillouin1d.