Function matmat 

pub fn matmat(
    operator: &(impl LinearOperator + ?Sized),
    columns: &[Vec<Complex64>],
) -> Result<Vec<Vec<Complex64>>>

Apply an operator to a column batch without constructing a dense matrix.