add tests

Author: lkdvos

test/Project.toml

@@ -5,6 +5,8 @@ ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
 DerivableInterfaces = "6c5e35bf-e59e-4898-b73c-732dcc4ba65f"
 Dictionaries = "85a47980-9c8c-11e8-2b9f-f7ca1fa99fb4"
 JLArrays = "27aeb0d3-9eb9-45fb-866b-73c2ecf80fcb"
+LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
+Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 SafeTestsets = "1bc83da4-3b8d-516f-aca4-4fe02f6d838f"
 SparseArraysBase = "0d5efcca-f356-4864-8770-e1ed8d78f208"
 Suppressor = "fd094767-a336-5f1f-9728-57cf17d0bbfb"

test/basics/test_linalg.jl

@@ -0,0 +1,35 @@
+using SparseArraysBase: SparseArrayDOK
+using LinearAlgebra: mul!
+using Random
+
+# TODO: add this to main package
+function sprand(::Type{T}, sz::Base.Dims; p::Real=0.5) where {T}
+  A = SparseArrayDOK{T}(undef, sz)
+  for I in eachindex(A)
+    if rand() < p
+      A[I] = rand(T)
+    end
+  end
+  return A
+end
+
+@testset "mul!" begin
+  T = Float64
+  szA = (2, 2)
+  szB = (2, 2)
+  szC = (szA[1], szB[2])
+
+  for p in 0.0:0.25:1
+    C = sprand(T, szC; p)
+    A = sprand(T, szA; p)
+    B = sprand(T, szB; p)
+
+    check1 = mul!(Array(C), Array(A), Array(B))
+    @test mul!(copy(C), A, B) ≈ check1
+
+    α = rand()
+    β = rand()
+    check2 = mul!(Array(C), Array(A), Array(B), α, β)
+    @test mul!(copy(C), A, B, α, β) ≈ check2
+  end
+end