We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 649ef09 commit e679b31Copy full SHA for e679b31
src/hessenberg.jl
@@ -29,7 +29,7 @@ is the first element of the output vector of H*x.
29
end
30
xnorm = norm(x[2:end])
31
α = x[1]
32
- if !iszero(xnorm) || (n == 1 && !iszero(α))
+ if !iszero(xnorm) || (!iszero(α))
33
β = (real(α) > 0 ? -1 : +1) * hypot(α,xnorm)
34
τ = 1 - α / β
35
α = 1 / (α - β)
@@ -38,7 +38,7 @@ is the first element of the output vector of H*x.
38
else
39
τ = T(0)
40
x .= 0
41
- β = real(T)(abs(α))
+ β = real(T)(0)
42
43
return τ, β
44
0 commit comments