Product integration methods based on discrete spline quasi-interpolants and application to weakly singular integral equations