  We investigate the algorithms for dynamical overlap fermions aiming at
improving the performance for large-scale simulations. We look for the best
combination of Hybrid Monte Carlo options and iterative quark solvers with
respect to the numerical costs. Our main target is a $N_f=2$ simulation with
overlap fermion on a $16^3\times 32$ lattice at lattice spacing around 0.12 fm.
