def getF_hess(**kwargs):
ia = IdealAnalytic(**kwargs)
return np.array([ia._F_x_y(x, y) for x, y in ('hh', 'ha', 'hb', 'ha', 'aa', 'ab', 'hb', 'ab', 'bb')]).reshape(
(3, 3))
def getF_jac(**kwargs):
ia = IdealAnalytic(**kwargs)
return np.array([ia._F_x(x) for x in 'hab'])