fn read_matrix<const R: usize, const C: usize>( reader: &mut Reader<'_>, ) -> Result<SMatrix<Complex, R, C>>