1pub fn reim_zero_ref(res: &mut [f64]) {
 2    res.fill(0.);
 3}
 4
 5pub fn reim_copy_ref(res: &mut [f64], a: &[f64]) {
 6    #[cfg(debug_assertions)]
 7    {
 8        assert_eq!(res.len(), a.len())
 9    }
10    res.copy_from_slice(a);
11}