An interior-point implementation developed and tuned for radiation therapy treatment planning

Sebastiaan Breedveld (), Bas Berg and Ben Heijmen ()
Additional contact information
Sebastiaan Breedveld: Erasmus University Medical Center - Cancer Institute
Bas Berg: Erasmus University Medical Center - Cancer Institute
Ben Heijmen: Erasmus University Medical Center - Cancer Institute

Computational Optimization and Applications, 2017, vol. 68, issue 2, No 2, 209-242

Abstract: Abstract While interior-point methods share the same fundamentals, the implementation determines the actual performance. In order to attain the highest efficiency, different applications may require differently tuned implementations. In this paper we describe an implementation specifically designed for optimisation in radiation therapy. These problems are large-scale nonlinear (and sometimes nonconvex) constrained optimisation problems, consisting of both sparse and dense data. Several application-specific properties are exploited to enhance efficiency. Permuting, tiling and mixed precision arithmetic allow the algorithm to optimally process the mixed dense and sparse data matrices (making this step 2.2 times faster, and overall runtime reduction of $$55\%$$ 55 % ) and scalability (16 threads resulted in a speed-up factor of 9.8 compared to singlethreaded performance, against a speed-up factor of 7.7 for the less optimised implementation). Predefined cost-functions are hard-coded and the computationally expensive second derivatives are written in canonical form, and combined if multiple cost-functions are defined for the same clinical structure. The derivatives are then computed using a scaled matrix–matrix product. A cheap initialisation strategy based on the background knowledge reduces the number of iterations by $$11\%$$ 11 % . We also propose a novel combined Mehrotra–Gondzio approach. The algorithm is extensively tested on a dataset consisting of 120 patients, distributed over 6 tumour sites/approaches. This test dataset is made publicly available.

Keywords: Nonlinear optimisation; Radiation therapy; Large-scale; Performance optimised; Interior-point; Tiled matrix algebra; Multiple precision arithmetic; Initialisation; Higher-order methods; 90C06; 90C25; 90C26; 90C30; 90C51; 90C90; 97M40; 97M60 (search for similar items in EconPapers)
Date: 2017
References: View references in EconPapers View complete reference list from CitEc
Citations:

Downloads: (external link)
http://link.springer.com/10.1007/s10589-017-9919-4 Abstract (text/html)
Access to the full text of the articles in this series is restricted.

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:spr:coopap:v:68:y:2017:i:2:d:10.1007_s10589-017-9919-4

Ordering information: This journal article can be ordered from
http://www.springer.com/math/journal/10589

DOI: 10.1007/s10589-017-9919-4

Access Statistics for this article

Computational Optimization and Applications is currently edited by William W. Hager

More articles in Computational Optimization and Applications from Springer
Bibliographic data for series maintained by Sonal Shukla () and Springer Nature Abstracting and Indexing ().