From 5f53bd64a5e4cf5e36b309f1914d9cb0ba2b65b8 Mon Sep 17 00:00:00 2001
From: MichelJuillard <michel.juillard@mjui.fr>
Date: Thu, 14 Oct 2021 22:20:36 +0200
Subject: [PATCH] pivoted Qr
---
src/QrAlgo.jl | 5 +++--
src/SchurAlgo.jl | 2 +-
test/QrAlgo_test.jl | 2 +-
3 files changed, 5 insertions(+), 4 deletions(-)
diff --git a/src/QrAlgo.jl b/src/QrAlgo.jl
index bcd89db..a62522c 100644
--- a/src/QrAlgo.jl
+++ b/src/QrAlgo.jl
@@ -160,7 +160,8 @@ for (geqrf, ormqr, elty) in
end
end
end
-
+=#
+
struct QrpWs{T <: Number} <: QR
tau::Vector{T}
jpvt::Vector{BlasInt}
@@ -210,5 +211,5 @@ for (geqp3, elty) in
t3 = Adjoint{$elty, <: StridedMatrix}
end
end
-=#
+
end
diff --git a/src/SchurAlgo.jl b/src/SchurAlgo.jl
index 4c7ceb2..e18354b 100644
--- a/src/SchurAlgo.jl
+++ b/src/SchurAlgo.jl
@@ -31,7 +31,7 @@ function mycompare(alphar_::Ptr{T}, alphai_::Ptr{T}, beta_::Ptr{T})::Cint where
return convert(Cint, ((alphar*alphar + alphai*alphai) < criterium*beta*beta) ? 1 : 0)
end
-mutable struct DgeesWs
+struct DgeesWs
jobvs::Ref{UInt8}
sdim::Ref{BlasInt}
wr::Vector{Float64}
diff --git a/test/QrAlgo_test.jl b/test/QrAlgo_test.jl
index 5729d19..84756b0 100644
--- a/test/QrAlgo_test.jl
+++ b/test/QrAlgo_test.jl
@@ -11,7 +11,7 @@ n = 10
QrAlgo.geqp3!(A, ws)
- target = qr(A0, Val(true))
+ target = qr(A0, ColumnNorm())
#display(triu(A))
#display(triu(target.R))
--
GitLab