  We studied the nonequilibrium aging behavior of the Random Field Ising Model
in three dimensions for various values of the disorder strength. This allowed
us to investigate how the aging behavior changes across the
ferromagnetic-paramagnetic phase transition. We investigated a large system
size of $N=256^3$ spins and up to $10^8$ Monte Carlo sweeps. To reach these
necessary long simulation times we employed an implementation running on Intel
Xeon Phi coprocessors, reaching single spin flip times as short as 6 ps. We
measured typical correlation functions in space and time to extract a growing
length scale and corresponding exponents.
