Method row_reduce

src/complex/matrix.rs:1772–1781

(&self, f: F)

    }

    fn row_reduce<F>(&self, f: F) -> Vec<C64>
    where
        F: Fn(Vec<C64>) -> C64,
    {
        let mut v = vec![Complex::zero(); self.row];
        for (i, slot) in v.iter_mut().enumerate().take(self.row) {
            *slot = f(self.row(i));
        }
        v
    }
}

pub fn cdiag(n: usize) -> ComplexMatrix {