I introduce the HESSIAN (highly efficient simulation smoothing in a nutshell) method for numerically efficient simulation smoothing in state space models with univariate states. Given a vector θ of parameters, the vector of states α=(α1,…,αn) is Gaussian and the observed vector...