def fixed_hessian(point, vars=None, model=None):
 model = modelcontext(model)
 if (vars is None):
  vars = model.cont_vars
 vars = inputvars(vars)
 point = Point(point, model=model)
 bij = DictToArrayBijection(ArrayOrdering(vars), point)
 dlogp = bij.mapf(model.fastdlogp(vars))
 rval = (np.ones(bij.map(point).size) / 10)
 return rval

def fixed_hessian ( point , vars = None , model = None ) : model = modelcontext ( model ) if ( vars is None ) : vars = model . cont_vars vars = inputvars ( vars ) point = Point ( point , model = model ) bij = DictToArrayBijection ( ArrayOrdering ( vars ) , point ) dlogp = bij . mapf ( model . fastdlogp ( vars ) ) rval = ( np . ones ( bij . map ( point ) . size ) / 10 ) return rval
returns a fixed hessian for any chain location .