Function sp2md 

pub fn sp2md<T: ComplexScalar>(
    op_a: Trans,
    descr_a: &MatDescr,
    a: &ComplexSparseMatrix<T>,
    op_b: Trans,
    descr_b: &MatDescr,
    b: &ComplexSparseMatrix<T>,
    alpha: T,
    beta: T,
    c: &mut [T],
    layout: Order,
    ldc: usize,
) -> Result<()>

Compute C := α · op_A(A) · op_B(B) + β · C for complex sparse A and B, dense complex C.