Parallel multigrid waveform relaxation for parabolic problems