Estimating upcrossing FPT densities via simulation of Gaussian processes