Computing reachable sets via barrier methods on SIMD architectures

Title data

Grüne, Lars ; Jahn, Thomas U.:
Computing reachable sets via barrier methods on SIMD architectures.
In: Eberhardsteiner, Josef ; Böhm, Helmut J. ; Rammerstorfer, Franz G. (eds.): Proceedings of the 6th European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 2012) Held at the University of Vienna, Austria, September 10-14, 2012. - Vienna : Vienna University of Technology , 2012 . - pp. 2076-2095
ISBN 9783950248197

Related URLs

Abstract in another language

We consider the problem of computing reachable sets of ODE-based control systems parallely on CUDA hardware. To this end, we modify an existing algorithm based on solving optimal control problems.
The idea is to simplify the optimal control problems to pure feasibility problems instead of minimizing an objective function. We show that an interior point algorithm is well suited for solving the resulting feasibility problems and leads to a sequence of linear systems of equations with identical matrix layout. If the problem is defined properly, these matrices are sparse and can be transformed into a hierarchical lower arrow form which can be solved on CUDA with sparse linear algebra and Cholesky’s method.
We demonstrate the performance of our new algorithm by computing the reachable sets of two test problems on a CPU implementation using several explicit and implicit Runge-Kutta methods of different order. The experiments reveal a significant speedup compared to the original optimal control algorithm.