Bilinear discretization of quadratic vector fields: integrability and geometry