GLISSE: A GPU-optimized planetary system integrator with application to orbital stability calculations.

•A graphical processing unit (GPU)-optimized integrator, about 100 times faster than a single central processing unit (CPU) core, is presented•The hardware architecture and algorithm allow orbital stability calculations for large numbers of test particles•We illustrate the integrator’s usage on the example problems of the Uranus-Neptune belt andtwo transneptunian mean motion resonances