An O(log N) Parallel Algorithm for Newton Step Computations with Applications to Moving Horizon Estimation
2016 (English)In: Proceedings of the 2016 European Control Conference, 2016, 1630-1636 p.Conference paper (Refereed)
In Moving Horizon Estimation (MHE) the computed estimate is found by solving a constrained finite-time optimal estimation problem in real-time at each sample in a receding horizon fashion. The constrained estimation problem can be solved by, e.g., interior-point (IP) or active-set (AS) methods, where the main computational effort in both methods is known to be the computation of the search direction, i.e., the Newton step. This is often done using generic sparsity exploiting algorithms or serial Riccati recursions, but as parallel hardware is becoming more commonly available the need for parallel algorithms for computing the Newton step is increasing. In this paper a newly developed tailored, non-iterative parallel algorithm for computing the Newton step using the Riccati recursion for Model Predictive Control (MPC) is extended to MHE problems. The algorithm exploits the special structure of the Karush-Kuhn-Tucker system for the optimal estimation problem. As a result it is possible to obtain logarithmic complexity growth in the estimation horizon length, which can be used to reduce the computation time for IP and AS methods when applied to what is today considered as challenging estimation problems. Furthermore, promising numerical results have been obtained using an ANSI-C implementation of the proposed algorithm, which uses Message Passing Interface (MPI) together with InfiniBand and is executed on true parallel hardware. Beyond MHE, due to similarities in the problem structure, the algorithm can be applied to various forms of on-line and off-line smoothing problems.
Place, publisher, year, edition, pages
2016. 1630-1636 p.
Parallel MHE, Moving Horizon Estimation, Riccati Recursion
IdentifiersURN: urn:nbn:se:liu:diva-129995OAI: oai:DiVA.org:liu-129995DiVA: diva2:945951
2016 European Control Conference, Aalborg, Denmark, June 29 - July 1, 2016.