1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
# This file is a part of Julia. License is MIT: https://julialang.org/license
module TestLDLT
using Test, LinearAlgebra, Random
Random.seed!(123)
@testset "Factorization conversions of LDLT" begin
S = SymTridiagonal(randn(5), randn(4))
F = ldlt(S)
@test Factorization{eltype(S)}(F) === F
@test Array(Factorization{complex(eltype(S))}(F)) ≈ Array(ldlt(complex(S)))
@test eltype(Factorization{complex(eltype(S))}) == complex(eltype(S))
end
@testset "REPL printing of LDLT" begin
S = SymTridiagonal(randn(5), randn(4))
F = ldlt(S)
ldltstring = sprint((t, s) -> show(t, "text/plain", s), F)
lstring = sprint((t, s) -> show(t, "text/plain", s), F.L)
dstring = sprint((t, s) -> show(t, "text/plain", s), F.D)
@test ldltstring == "$(summary(F))\nL factor:\n$lstring\nD factor:\n$dstring"
end
end # module TestLDLT
|
