Numerical Methods for Total Variation-based Deblurring: a Comparison in a Parallel Environmen