We study aspects concerning numerical simulations of lattice QCD with two f
lavors of dynamical Ginsparg-Wilson quarks with degenerate masses. A hybrid
Monte Carlo algorithm is described and a formula for the fermionic force i
s derived for two specific implementations. The implementation with the opt
imal rational approximation method is favored in both CPU time and memory c
onsumption. (C) 1999 Published by Elsevier Science B.V. All rights reserved
.