Parallel and accurate Poisson disk sampling on arbitrary surfaces