Hybrid Jacobian and Gauss-Seidel proximal block coordinate update methods for linearly constrained convex programming