Accelerating total variation regularization for matrix-valued images on GPUs