def find_hessian_diag(point, vars=None, model=None):
 model = modelcontext(model)
 H = model.fastfn(hessian_diag(model.logpt, vars))
 return H(Point(point, model=model))

def find_hessian_diag ( point , vars = None , model = None ) : model = modelcontext ( model ) H = model . fastfn ( hessian_diag ( model . logpt , vars ) ) return H ( Point ( point , model = model ) )
returns hessian of logp at the point passed .