Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Added documentation for Trust region radius update schemes #199

Merged
merged 3 commits into from
Aug 17, 2023
Merged
Show file tree
Hide file tree
Changes from 2 commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
10 changes: 10 additions & 0 deletions docs/src/api/nonlinearsolve.md
Original file line number Diff line number Diff line change
Expand Up @@ -8,3 +8,13 @@ These are the native solvers of NonlinearSolve.jl.
NewtonRaphson
TrustRegion
```

### Radius Update Schemes for Trust Region (RadiusUpdateSchemes)

```@docs
RadiusUpdateSchemes.Simple
RadiusUpdateSchemes.Hei
RadiusUpdateSchemes.Yuan
RadiusUpdateSchemes.Bastin
RadiusUpdateSchemes.Fan
```
44 changes: 44 additions & 0 deletions src/trustRegion.jl
Original file line number Diff line number Diff line change
@@ -1,8 +1,52 @@
EnumX.@enumx RadiusUpdateSchemes begin
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Add a docstring to the enum type as well describing what these do.

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

okay got it

"""
`RadiusUpdateSchemes.Simple`

The simple or conventional radius update scheme. This scheme is chosen by default
and follows the conventional approach to update the trust region radius, i.e. if the
trial step is accepted it increases the radius by a fixed factor (bounded by a maximum radius)
and if the trial step is rejected, it shrinks the radius by a fixed factor.
"""
Simple

"""
`RadiusUpdateSchemes.Hei`

This scheme is proposed by [Hei, L.] (https://www.jstor.org/stable/43693061). The trust region radius
depends on the size (norm) of the current step size. The hypothesis is to let the radius converge to zero
as the iterations progress, which is more reliable and robust for ill-conditioned as well as degenerate
problems.
"""
Hei

"""
`RadiusUpdateSchemes.Yuan`

This scheme is proposed by [Yuan, Y.] (https://www.researchgate.net/publication/249011466_A_new_trust_region_algorithm_with_trust_region_radius_converging_to_zero).
Similar to Hei's scheme, the trust region is updated in a way so that it converges to zero, however here,
the radius depends on the size (norm) of the current gradient of the objective (merit) function. The hypothesis
is that the step size is bounded by the gradient size, so it makes sense to let the radius depend on the gradient.
"""
Yuan

"""
`RadiusUpdateSchemes.Bastin`

This scheme is proposed by [Bastin, et al.] (https://www.researchgate.net/publication/225100660_A_retrospective_trust-region_method_for_unconstrained_optimization).
The scheme is called a retrospective update scheme as it uses the model function at the current
iteration to compute the ratio of the actual reduction and the predicted reduction in the previous
trial step, and use this ratio to update the trust region radius. The hypothesis is to exploit the information
made available during the optimization process in order to vary the accuracy of the objective function computation.
"""
Bastin

"""
`RadiusUpdateSchemes.Fan`

This scheme is proposed by [Fan, J.] (https://link.springer.com/article/10.1007/s10589-005-3078-8). It is very much similar to
Hei's and Yuan's schemes as it lets the trust region radius depend on the current size (norm) of the objective (merit)
function itself. These new update schemes are known to improve local convergence.
"""
Fan
end

Expand Down