Cluster Newton method for sampling multiple solutions of underdetermined inverse problems: Application to a parameter identification problem in pharmacokinetics

Yasunori Aoki, Ken Hayami, Hans De Sterck, Akihiko Konagaya

Research output: Contribution to journalArticleResearchpeer-review

18 Citations (Scopus)


A new algorithm is proposed for simultaneously finding multiple solutions of an underdetermined inverse problem. The algorithm was developed for an ODE parameter identification problem in pharmacokinetics for which multiple solutions are of interest. The algorithm proceeds by computing a cluster of solutions simultaneously, and is more efficient than algorithms that compute multiple solutions one-by-one because it fits the Jacobian in a collective way using a least squares approach. It is demonstrated numerically that the algorithm finds accurate solutions that are suitably distributed, guided by a priori information on which part of the solution set is of interest, and that it does so much more efficiently than a baseline Levenberg Marquardt method that computes solutions one-by-one. It is also demonstrated that the algorithm benefits from improved robustness due to an inherent smoothing provided by the least-squares fitting.

Original languageEnglish
Pages (from-to)B14–B44
Number of pages31
JournalSIAM Journal on Scientific Computing
Issue number1
Publication statusPublished - 2014
Externally publishedYes


  • Inverse problems
  • Method of least squares
  • Pharmacokinetics
  • Underdetermined problems

Cite this