Changed comparison to account for floating point errors.

Author: loiseaujc

src/stdlib_sparse_spmv.fypp

@@ -620,9 +620,9 @@ contains
         #:endif
 
         ! Deal with optional arguments.
-        alpha_ = 1.0_${k1}$ ; if (present(alpha)) alpha_ = alpha
-        beta_  = 0.0_${k1}$ ; if (present(beta))  beta_  = beta
-        op_    = "N"        ; if (present(op))    op_    = op
+        alpha_ = 1.0_${k1}$     ; if (present(alpha)) alpha_ = alpha
+        beta_  = 0.0_${k1}$     ; if (present(beta))  beta_  = beta
+        op_    = sparse_op_none ; if (present(op))    op_    = op
 
         ! Prepare Lapack arguments.
         n = A%nrows ; ldx = n ; ldy = n ; y = 0.0_${k1}$

test/linalg/test_linalg_sparse.fypp

@@ -4,6 +4,8 @@
 module test_sparse_spmv
     use testdrive, only : new_unittest, unittest_type, error_type, check, skip_test
     use stdlib_kinds, only: sp, dp, xdp, qp, int8, int16, int32, int64
+    use stdlib_linalg, only: hermitian
+    use stdlib_math, only: all_close
     use stdlib_sparse
 
     implicit none
@@ -408,14 +410,22 @@ contains
 
             ! Test y = A @ x
             y1 = matmul(Amat, x) ; call spmv(A, x, y2)
-            call check(error, all(y1 == y2))
+            call check(error, all_close(y1, y2))
             if (allocated(error)) return
 
             ! Test y = A.T @ x
             y1 = 0.0_wp ; y2 = 0.0_wp
             y1 = matmul(transpose(Amat), x) ; call spmv(A, x, y2, op="T")
-            call check(error, all(y1 == y2))
+            call check(error, all_close(y1, y2))
             if (allocated(error)) return
+
+            #:if t1.startswith('complex') 
+            ! Test y = A.H @ x
+            y1 = 0.0_wp ; y2 = 0.0_wp
+            y1 = matmul(hermitian(Amat), x) ; call spmv(A, x, y2, op="H")
+            call check(error, all_close(y1, y2))
+            if (allocated(error)) return
+            #:endif
         end block
         #:endif
         #:endfor