Solving Singular Control Problems in Mathematical Biology, Using PASA (2010.06744v1)

Abstract: In this paper, we will demonstrate how to use a nonlinear polyhedral constrained optimization solver called the Polyhedral Active Set Algorithm (PASA) for solving a general singular control problem. We present methods of discretizing a general optimal control problem that involves the use of the gradient of the Lagrangian for computing the gradient of the cost functional so that PASA can be applied. When a numerical solution contains artifacts that resemble "chattering'', a phenomenon where the control oscillates wildly along the singular region, we recommend a method of regularizing the singular control problem by adding a term to the cost functional that measures a scalar multiple of the total variation of the control, where the scalar is viewed as a tuning parameter. We then demonstrate PASA's performance on three singular control problems that give rise to different applications of mathematical biology. We also provide some exposition on the heuristics that we use in determining an appropriate size for the tuning parameter.