A partially inexact generalized primal-dual hybrid gradient method for saddle point problems with bilinear couplings