Hessian-vector product or jacobian-vector product with autodiff

It is often convenient for such large-scale problems to evaluate matrix-vector products instead of evaluating the whole matrix every time. Have you considered using automatic differentiation for such operations?