Skip to content

Commit

Permalink
Update krylov_utils.jl
Browse files Browse the repository at this point in the history
  • Loading branch information
amontoison authored Nov 5, 2024
1 parent 21d584d commit 906537b
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion src/krylov_utils.jl
Original file line number Diff line number Diff line change
Expand Up @@ -390,7 +390,7 @@ function to_boundary(n :: Int, x :: AbstractVector{FC}, d :: AbstractVector{FC},

# q₂ = ‖d‖², q₁ = xᴴd + dᴴx, q₀ = ‖x‖² - Δ²
# ‖x‖² ≤ Δ² ⟹ (q₁)² - 4 * q₂ * q₀ ≥ 0
roots = roots_quadratic(dNorm2, 2 * rxd, xNorm2 - radius2)
roots = roots_quadratic(dNorm2, 2 * real(rxd), xNorm2 - radius2)

return roots # `σ1` and `σ2`
end
Expand Down

0 comments on commit 906537b

Please sign in to comment.