Function add 

pub fn add<T: ComplexScalar>(
    op: Trans,
    a: &ComplexSparseMatrix<T>,
    alpha: T,
    b: &ComplexSparseMatrix<T>,
) -> Result<ComplexSparseMatrix<T>>

Compute C := α · op(A) + B returning a fresh complex CSR matrix.