Closed talassio closed 5 years ago
Consider adding a fallback call for type_of_method = 2 and nlls_method = 3|4 if call regularization_solver fails, then don't terminate with error but try call solve_galahad.
type_of_method = 2
nlls_method = 3|4
call regularization_solver
call solve_galahad
https://github.com/ralna/RALFit/blob/b9f438c25874379a71922f7c18f352722d7cf55c/libRALFit/src/ral_nlls_internal.f90#L1192-L1207
The issue is raised when for some reason Delta is such that Ashift = A + 1/Delta does not guarantee Ashift is SPD and the call to SPD_solve fails.
Ashift = A + 1/Delta
Ashift
SPD_solve
Consider adding a fallback call for
type_of_method = 2andnlls_method = 3|4ifcall regularization_solverfails, then don't terminate with error but trycall solve_galahad.https://github.com/ralna/RALFit/blob/b9f438c25874379a71922f7c18f352722d7cf55c/libRALFit/src/ral_nlls_internal.f90#L1192-L1207
The issue is raised when for some reason Delta is such that
Ashift = A + 1/Deltadoes not guaranteeAshiftis SPD and the call toSPD_solvefails.