Interacting electrons and non-adiabatic Holstein phonons: The numerical treatment of the neutral-ionic phase transition