importsparsefromsparse._compressedimportCSR,CSCgcxs_a=sparse.random((10000,20000),.01,format="gcxs")gcxs_b=sparse.random((10000,20000),.01,format="gcxs")csr_a=CSR(gcxs_a)csr_b=CSR(gcxs_b)scipy_a=csr_a.to_scipy_sparse()scipy_b=csr_b.to_scipy_sparse()

%timeitgcxs_a+gcxs_b# 999 ms ± 10.2 ms per loop (mean ± std. dev. of 7 runs, 1 loop each)%timeitcsr_a+csr_b# 32.2 ms ± 1.9 ms per loop (mean ± std. dev. of 7 runs, 10 loops each)%timeitscipy_a+scipy_b# 27.5 ms ± 894 µs per loop (mean ± std. dev. of 7 runs, 10 loops each)