Simulation of linear stationary Hawkes processes: the subexponential case