Fast numerical solvers for parameter identification problems in mathematical biology

Karolína Benková, John W. Pearson, Mariya Ptashnyk

Research output: Working paperPreprint

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.
Original languageEnglish
PublisherarXiv
DOIs
Publication statusPublished - 27 Aug 2024

Keywords

  • math.NA
  • cs.NA
  • math.OC
  • 49M41, 92C15, 65M22, 65M60, 65F08, 65F10
  • PDE-constrained optimization
  • parameter idenitification
  • pattern formation
  • time-stepping
  • Krylov subspace methods
  • preconditioning
  • 49M41
  • 92C15
  • 65M22
  • 65M60
  • 65F08
  • 65F10

Fingerprint

Dive into the research topics of 'Fast numerical solvers for parameter identification problems in mathematical biology'. Together they form a unique fingerprint.

Cite this