[Code] Cuda csr matrix: accum logic fixes

Author: EgorOrachyov

cubool/include/cubool/cubool.h

@@ -91,7 +91,9 @@ typedef enum cuBool_Hint {
     /** Logging hint: log includes all types of messages */
     CUBOOL_HINT_LOG_ALL = 0x128,
     /** No duplicates in the build data */
-    CUBOOL_HINT_NO_DUPLICATES = 0x256
+    CUBOOL_HINT_NO_DUPLICATES = 0x256,
+    /** Performs time measurement and logs elapsed operation time */
+    CUBOOL_HINT_TIME_CHECK = 0x512
 } cuBool_Hint;
 
 /** Hit mask */

cubool/sources/cuda/matrix_csr_multiply.cu

@@ -40,26 +40,27 @@ namespace cubool {
         assert(this->getNrows() == M);
         assert(this->getNcols() == N);
 
+        if (!accumulate) {
+            // Clear all values
+            this->mMatrixImpl.zero_dim();
+        }
+
         if (a->isMatrixEmpty() || b->isMatrixEmpty()) {
-            // A or B has no values
+            // Return empty matrix
             return;
         }
 
-        CHECK_RAISE_ERROR(accumulate, NotImplemented, "Supported only accumulated multiplication");
+        // Ensure csr proper csr format even if empty
+        a->resizeStorageToDim();
+        b->resizeStorageToDim();
+        this->resizeStorageToDim();
 
-        if (accumulate) {
-            // Ensure csr proper csr format even if empty
-            a->resizeStorageToDim();
-            b->resizeStorageToDim();
-            this->resizeStorageToDim();
+        // Call backend r = c + a * b implementation, as C this is passed
+        nsparse::spgemm_functor_t<bool, index, DeviceAlloc<index>> spgemmFunctor;
+        auto result = spgemmFunctor(mMatrixImpl, a->mMatrixImpl, b->mMatrixImpl);
 
-            // Call backend r = c + a * b implementation, as C this is passed
-            nsparse::spgemm_functor_t<bool, index, DeviceAlloc<index>> spgemmFunctor;
-            auto result = spgemmFunctor(mMatrixImpl, a->mMatrixImpl, b->mMatrixImpl);
-
-            // Assign result to this
-            this->mMatrixImpl = std::move(result);
-        }
+        // Assign result to this
+        this->mMatrixImpl = std::move(result);
     }
 
 }