We present our latest results on the application of the complex Langevin method to one-and twodimensional QCD. Although the method is stable, it unfortunately converges to an incorrect result when applied as such. After applying additional gauge cooling steps, the results agree with the known analytical results in the one-dimensional case. However, in the two-dimensional case the disagreement subsists, even with gauge cooling, when the sign problem is sufficiently large.