Lossy Compression in Optimal Control of Cardiac Defibrillation
- This paper presents efficient computational techniques for solving an optimization problem in cardiac defibrillation governed by the monodomain equations. Time-dependent electrical currents injected at different spatial positions act as the control. Inexact Newton-CG methods are used, with reduced gradient computation by adjoint solves. In order to reduce the computational complexity, adaptive mesh refinement for state and adjoint equations is performed. To reduce the high storage and bandwidth demand imposed by adjoint gradient and Hessian-vector evaluations, a lossy compression technique for storing trajectory data is applied. An adaptive choice of quantization tolerance based on error estimates is developed in order to ensure convergence. The efficiency of the proposed approach is demonstrated on numerical examples.
Author: | Sebastian GötschelORCiD, Chamakuri Nagaiah, Karl Kunisch, Martin WeiserORCiD |
---|---|
Document Type: | Article |
Parent Title (English): | J. Sci. Comput. |
Volume: | 60 |
Issue: | 1 |
First Page: | 35 |
Last Page: | 59 |
Year of first publication: | 2014 |
Page Number: | 25 |
Preprint: | urn:nbn:de:0297-zib-18566 |
DOI: | https://doi.org/10.1007/s10915-013-9785-x |