Fast numerical solvers for parameter identification problems in mathematical biology
Abstract
In this paper, we consider effective discretization strategies and iterative solvers for nonlinear PDE-constrained optimization models for pattern evolution within biological processes. Upon a Sequential Quadratic Programming linearization of the optimization problem, we devise appropriate time-stepping schemes and discrete approximations of the cost functionals such that the discretization and optimization operations are commutative, a highly desirable property of a discretization of such problems. We formulate the large-scale, coupled linear systems in such a way that efficient preconditioned iterative methods can be applied within a Krylov subspace solver. Numerical experiments demonstrate the viability and efficiency of our approach.
- Publication:
-
arXiv e-prints
- Pub Date:
- August 2024
- DOI:
- 10.48550/arXiv.2408.14926
- arXiv:
- arXiv:2408.14926
- Bibcode:
- 2024arXiv240814926B
- Keywords:
-
- Mathematics - Numerical Analysis;
- Mathematics - Optimization and Control;
- 49M41;
- 92C15;
- 65M22;
- 65M60;
- 65F08;
- 65F10
- E-Print:
- 33 pages, 3 figures, 7 tables