Skip to content
Snippets Groups Projects
Commit 5f53bd64 authored by MichelJuillard's avatar MichelJuillard
Browse files

pivoted Qr

parent b329510b
No related branches found
No related tags found
No related merge requests found
...@@ -160,6 +160,7 @@ for (geqrf, ormqr, elty) in ...@@ -160,6 +160,7 @@ for (geqrf, ormqr, elty) in
end end
end end
end end
=#
struct QrpWs{T <: Number} <: QR struct QrpWs{T <: Number} <: QR
tau::Vector{T} tau::Vector{T}
...@@ -210,5 +211,5 @@ for (geqp3, elty) in ...@@ -210,5 +211,5 @@ for (geqp3, elty) in
t3 = Adjoint{$elty, <: StridedMatrix} t3 = Adjoint{$elty, <: StridedMatrix}
end end
end end
=#
end end
...@@ -31,7 +31,7 @@ function mycompare(alphar_::Ptr{T}, alphai_::Ptr{T}, beta_::Ptr{T})::Cint where ...@@ -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) return convert(Cint, ((alphar*alphar + alphai*alphai) < criterium*beta*beta) ? 1 : 0)
end end
mutable struct DgeesWs struct DgeesWs
jobvs::Ref{UInt8} jobvs::Ref{UInt8}
sdim::Ref{BlasInt} sdim::Ref{BlasInt}
wr::Vector{Float64} wr::Vector{Float64}
......
...@@ -11,7 +11,7 @@ n = 10 ...@@ -11,7 +11,7 @@ n = 10
QrAlgo.geqp3!(A, ws) QrAlgo.geqp3!(A, ws)
target = qr(A0, Val(true)) target = qr(A0, ColumnNorm())
#display(triu(A)) #display(triu(A))
#display(triu(target.R)) #display(triu(target.R))
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment