  The complex Langevin method has been attracting much attention as a solution
to the sign problem since the method was shown to work in finite density QCD in
the deconfined phase by using the so-called gauge cooling procedure. Whether it
works also in the confined phase with light quarks is still an open question,
though. In order to shed light on this question, we apply the method to the
chiral Random Matrix Theory, which describes the epsilon regime of finite
density QCD. Earlier works reported that a naive implementation of the method
fails to reproduce the known exact results and that the problem can be solved
by choosing a suitable coordinate. In this work we stick to the naive
implementation, and show that a generalized gauge cooling procedure can be used
to avoid the problem.
