An optimal non-linear method for simulating relic neutrinos