code: def find_hessian(point, vars=None, model=None):
      """
      Returns Hessian of logp at the point passed.
  
      Parameters
      ----------
      model: Model (optional if in `with` context)
      point: dict
      vars: list
          Variables for which Hessian is to be calculated.
      """
      model = modelcontext(model)
      H = model.compile_d2logp(vars)
      return H(Point(point, filter_model_vars=True, model=model))
  
file_path: pymc/pymc/tuning/scaling.py
imports: import numpy
  from numpy import exp, log, sqrt
  from pymc.blocking import DictToArrayBijection
  from pymc.model import Point, modelcontext
  from pymc.pytensorf import hessian_diag
  from pymc.util import get_var_name
label_desc: Use image node_id 2 for calling a global function with example usage: find_hessian(point, vars, model) and returns: H