  We study aspects concerning numerical simulations of Lattice QCD with two
flavors of dynamical Ginsparg-Wilson quarks with degenerate masses. A Hybrid
Monte Carlo algorithm is described and the formula for the fermionic force is
derived for two specific implementations. The implementation with optimal
rational approximation method is favored both in CPU time and memory
consumption.
