diff --git a/test/test_pmtools.jl b/test/test_pmtools.jl index ee08a42..6b2200f 100644 --- a/test/test_pmtools.jl +++ b/test/test_pmtools.jl @@ -156,6 +156,11 @@ W = zeros(3,3); A2,E2,B2,C2,D2 = spm2ls(D,D,D,W,atol = 1.e-7, minimal=true) @test D2 ≈ -D +D = rand(3,3); +W = zeros(3,3,1); +A2,E2,B2,C2,D2 = spm2ls(D,D,D,W,atol = 1.e-7, minimal=true) +@test D2 ≈ -D + A2,E2,B2,C2,D2 = spm2ls(D,D,D,W,atol = 1.e-7) @test -C2*inv(A2)*B2+D2 ≈ -D && iszero(E2) @@ -165,6 +170,12 @@ W = zeros(3,3); sys = spm2lps(D,D,D,W,atol = 1.e-7, minimal=true) @test sys[7] ≈ -D && sys[8] ≈ 0*D +# simple test +D = rand(3,3); +W = zeros(3,3,1); +sys = spm2lps(D,D,D,W,atol = 1.e-7, minimal=true) +@test sys[7] ≈ -D && sys[8] ≈ 0*D + # D = rand(3,3,2); V = zeros(3,3,1); V[:,:,1] = Matrix{eltype(V)}(I,3,3); W = zeros(3,3); A2,E2,B2,C2,D2 = spm2ls(D,D,V,W,atol = 1.e-7, minimal=true)