test reuse sparsity pattern

Signed-off-by: Steven Hahn <[email protected]>

Author: quantumsteve

resolve/workspace/LinAlgWorkspaceCpu.cpp

@@ -108,6 +108,7 @@ namespace ReSolve
   {
     assert(scaleaddi_buffer_ == nullptr);
     scaleaddi_buffer_ = buffer;
+    scaleAddISetupDone();
   }
 
 } // namespace ReSolve

tests/unit/matrix/MatrixHandlerTests.hpp

@@ -214,7 +214,17 @@ namespace ReSolve
         {
           A->syncData(memory::HOST);
         }
-        status *= verifyScaleAddICsrMatrix(A, val);
+        // expected sum is 2*30+1
+        status *= verifyScaleAddICsrMatrix(A, 61);
+
+        // run again to reuse sparsity pattern.
+        handler_.scaleAddI(A, val, memspace_);
+        if (memspace_ == memory::DEVICE)
+        {
+          A->syncData(memory::HOST);
+        }
+        // expected sum is 2*61+1
+        status *= verifyScaleAddICsrMatrix(A, 123);
 
         delete A;
         return status.report(testname.c_str());
@@ -705,10 +715,10 @@ namespace ReSolve
        * @pre A is a valid, allocated CSR matrix
        * @invariant A
        * @param[in] A matrix::Csr* pointer to the matrix to be verified
-       * @param[in] alpha matrix scale factor
+       * @param[in] expected sum of row elements
        * @return bool true if the matrix is valid, false otherwise
        */
-      bool verifyScaleAddICsrMatrix(matrix::Csr* A, real_type alpha)
+      bool verifyScaleAddICsrMatrix(matrix::Csr* A, real_type expected)
       {
         // Check if the matrix is valid
         if (A == nullptr)
@@ -728,8 +738,6 @@ namespace ReSolve
         index_type* scaled_row_ptr  = A->getRowData(memory::HOST);
         real_type*  scaled_value    = A->getValues(memory::HOST);
 
-        const real_type expected = 30. * alpha + 1.;
-
         // Verify values - each element scaled by scale. Diagonal elements should be  1.
         for (index_type i = 0; i < n; ++i)
         {

tests/unit/matrix/runMatrixHandlerTests.cpp

@@ -76,8 +76,10 @@ void runTests(const std::string& backend, ReSolve::tests::TestingResults& result
   result += test.rightScale(1024, 2048);
   workspace.resetLinAlgWorkspace();
   result += test.rightScale(2048, 1024);
+#if !defined(RESOLVE_USE_CUDA) && !defined(RESOLVE_USE_HIP)
   workspace.resetLinAlgWorkspace();
   result += test.scaleAddI(100);
+#endif
   std::cout << "\n";
 }