diff --git a/test/QrAlgo_test.jl b/test/QrAlgo_test.jl new file mode 100644 index 0000000000000000000000000000000000000000..5729d192ce7a2d94caafc01a33061b4a732588ce --- /dev/null +++ b/test/QrAlgo_test.jl @@ -0,0 +1,25 @@ + +n = 10 + +#for elty in (Float32, Float64, ComplexF32, ComplexF64) + +# A0 = randn(elty, n, n) + A0 = randn(n, n) + + A = copy(A0) + ws = QrAlgo.QrpWs(A) + + QrAlgo.geqp3!(A, ws) + + target = qr(A0, Val(true)) + +#display(triu(A)) +#display(triu(target.R)) + + +A = randn(n, 15) +vA = zeros(n, 2) +vA .= view(A, :, [7, 9]) +ws = QrAlgo.QrWs(vA) +geqrf_core!(vA, ws) +ormqr_core!('L', vA', A, ws)