Function add 

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

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