Approximative inference for multivariate functional data on massively parallel processors