A symmetric primal-dual method with double extrapolation for composite convex optimization involving three functions