The autograd branch, which will be merged soon, supports repeated application of .backward (or, more conveniently, a new autograd.differentiate operator) and can compute the exact Hessian-vector product.