pub fn add<F>(a: &CsrMatrix<F>, b: &CsrMatrix<F>) -> SparseResult<CsrMatrix<F>>where F: Float + NumAssign + Sum + SparseElement + 'static + Debug,
Matrix addition