Approximating Hessians in multilevel unconstrained optimization