Numerical techniques for lattice QCD in the ϵ-regime