The complex Langevin method holds great potential for addressing theories with complex actions. However, real-time quantum field theories remain largely unsolved due to the challenging sign problem. In this work, we present our recent advancements in applying the complex Langevin method to SU(2) Yang-Mills theory in 3+1 dimensions. To enhance stability, we introduce a novel anisotropic kernel that allows for real-time simulations beyond the inverse temperature. We explicitly show the method's ability to reproduce thermal relations among various propagators on the Schwinger-Keldysh contour. This approach may open the door to calculate transport coefficients and other real-time observables from first principles.