Solvers

LSQR, RRQR, BRR & ARD solvers that solve the Ax=b problem with some regularisation are available via IPFitting module. solver_params() is the way these are define and selected in ACE1pack.

ACE1pack.solver_paramsFunction

solver_params(; kwargs...) : returns a dictionary containing the complete set of parameters required to construct one of the solvers. All parameters are passed as keyword argument and the kind of parameters required depend on "type".

QR Parameters

  • 'type = "qr"`

LSQR Parameters

  • type = "lsqr"
  • damp = 5e-3
  • atol = 1e-6
  • colim = 1e8
  • maxiter = 1e5
  • verbose = false

RRQR Parameters

  • type = "rrqr"
  • rtol = 1e-5

BBR

  • type = "bbr"
  • n_iter = 300
  • tol = 1e-3

ARD

  • type = "ard"
  • n_iter = 300
  • tol = 1e-3
  • threshold_lambda = 1e4
source

Additional regularizers (currently only "laplacian") are given via regularizer_params()

ACE1pack.regularizer_paramsFunction

regularizer_params(; type = "laplacian", kwargs...) : returns a dictionary containing the complete set of parameters required to construct one of the solvers. All parameters are passed as keyword argument and the kind of parameters required depend on "type".

LSQR Parameters (default)

  • type = "laplacian"
  • rlap_scal = 3.0
source