The extraction of hadron spectra from COMPASS data is a challenging task not only from the physics but also from the computing point of view. Large data sets for a multitude of final states and production mechanisms require fast and flexible tools that make efficient use of the available computing infrastructure. We present first results from a new framework for decay amplitude calculation and the status of the acceleration of PWA maximum-likelihood fits using GPUs.