|
1 | 1 | using StaticArrays, Test, LinearAlgebra, Random |
2 | 2 |
|
3 | | -Base.randn(::Type{BigFloat}) = BigFloat(randn(Float64)) |
4 | | -Base.randn(::Type{BigFloat}, I::Integer) = [randn(BigFloat) for i=1:I] |
5 | | -Base.randn(::Type{Int}) = rand(-9:9) |
6 | | -Base.randn(::Type{Int}, I::Integer) = [randn(Int) for i=1:I] |
7 | | -Base.randn(::Type{Complex{T}}) where T = Complex{T}(randn(T,2)...) |
8 | | -Base.randn(::Type{Complex}) = randn(Complex{Float64}) |
| 3 | +broadenrandn(::Type{BigFloat}) = BigFloat(randn(Float64)) |
| 4 | +broadenrandn(::Type{Int}) = rand(-9:9) |
| 5 | +broadenrandn(::Type{Complex{T}}) where T = Complex{T}(broadenrandn(T), broadenrandn(T)) |
| 6 | +broadenrandn(::Type{T}) where T = randn(T) |
9 | 7 |
|
10 | 8 | Random.seed!(42) |
11 | 9 | @testset "QR decomposition" begin |
@@ -46,7 +44,7 @@ Random.seed!(42) |
46 | 44 | for eltya in (Float32, Float64, BigFloat, Int), |
47 | 45 | rel in (real, complex), |
48 | 46 | sz in [(3,3), (3,4), (4,3)] |
49 | | - arr = SMatrix{sz[1], sz[2], rel(eltya), sz[1]*sz[2]}( [randn(rel(eltya)) for i = 1:sz[1], j = 1:sz[2]] ) |
| 47 | + arr = SMatrix{sz[1], sz[2], rel(eltya), sz[1]*sz[2]}( [broadenrandn(rel(eltya)) for i = 1:sz[1], j = 1:sz[2]] ) |
50 | 48 | test_qr(arr) |
51 | 49 | end |
52 | 50 | # some special cases |
|
0 commit comments