Add tests

Author: janbruedigam

Project.toml

@@ -8,10 +8,14 @@ LightGraphs = "093fc24a-ae57-5d10-9952-331d41423f4d"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 
 [compat]
-julia = "1"
+LightGraphs = "1.3"
+StaticArrays = "1.0"
+julia = "1.5, 1.6"
 
 [extras]
+LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
+SafeTestsets = "1bc83da4-3b8d-516f-aca4-4fe02f6d838f"
 
 [targets]
-test = ["Test"]
+test = ["LinearAlgebra", "Test", "SafeTestsets"]

src/GraphBasedSystems.jl

@@ -8,8 +8,6 @@ export System,
     full_matrix,
     full_vector,
 
-    ldu_factorization!,
-    ldu_backsubstitution!,
     ldu_solve!
 
 
@@ -19,6 +17,9 @@ include("entry.jl")
 include("system.jl")
 include("graph_functions.jl")
 
+# include("lu.jl")
+# include("llt.jl")
+# include("ldlt.jl")
 include("ldu.jl")
 
 end

src/entry.jl

@@ -11,5 +11,5 @@ end
 
 function randomize!(entry::Entry)
     value = entry.value
-    entry.value = rand(eltype(value), size(value))
+    entry.value = randn(eltype(value), size(value))
 end

src/ldlt.jl

@@ -0,0 +1 @@
+# LDLᵀ factorization for block-symmetric systems

src/ldu.jl

@@ -1,3 +1,5 @@
+# LDU factorization for unsymmetric systems
+
 function ldu_factorization_l!(offdiagonal, diagonal)
     offdiagonal.value = offdiagonal.value/diagonal.value
     return

src/llt.jl

@@ -0,0 +1 @@
+# LLᵀ (Choleski) factorization for symmetric systems

src/lu.jl

@@ -0,0 +1 @@
+# LU factorization for unsymmetric systems

test/ldu_test.jl

@@ -0,0 +1,46 @@
+using GraphBasedSystems
+using LinearAlgebra
+GBS = GraphBasedSystems
+
+Z = zeros(Int64,10,10)
+A = [
+    0 1 0 1 1 0 1 0 1 0
+    1 0 1 0 0 0 0 0 0 0
+    0 1 0 0 0 0 0 0 0 0
+    1 0 0 0 1 0 0 0 0 0
+    1 0 0 1 0 1 0 0 0 0
+    0 0 0 0 1 0 0 0 0 0
+    1 0 0 0 0 0 0 1 0 0
+    0 0 0 0 0 0 1 0 1 1
+    1 0 0 0 0 0 0 1 0 0
+    0 0 0 0 0 0 0 1 0 0
+]
+
+A = [
+    A Z Z Z
+    Z A Z Z
+    Z Z A Z
+    Z Z Z A
+]
+
+A[3,11] = A[11,3] = 1
+A[6,21] = A[21,6] = 1
+A[10,31] = A[31,10] = 1
+
+
+system = System{Float64}(A, [rand(1:3) for i=1:size(A)[1]])
+
+for i=1:10
+    for entry in collect(values(system.matrix_entries))
+        GBS.randomize!(entry)
+    end
+    for entry in collect(values(system.vector_entries))
+        GBS.randomize!(entry)
+    end
+
+    B = full_matrix(system)
+    b = full_vector(system)
+    ldu_solve!(system)
+    @test maximum(abs.(full_vector(system)-B\b)) < 1e-5
+end
+

test/runtests.jl

@@ -1,6 +1,6 @@
-using GraphBasedSystems
 using Test
+using SafeTestsets
 
-@testset "GraphBasedSystems.jl" begin
-    # Write your tests here.
+@safetestset "LDU Tests" begin
+    include("ldu_test.jl")
 end