def _test_first_derivative(name): x = np.linspace(0.0001, 0.98, 5) h = _default_base_step(x, scale=2) f, df = get_function(name, n=1) der = f(bicomplex(x + h * 1j, 0)).imag1 / h der_true = df(x) np.testing.assert_allclose(der, der_true, err_msg=('%s' % name))
def _test_second_derivative(name): x = np.linspace(0.01, 0.98, 5) h = _default_base_step(x, scale=2.5) # h = 1e-8 f, df = get_function(name, n=2) der = f(bicomplex(x + h * 1j, h)).imag12 / h**2 der_true = df(x) np.testing.assert_allclose(der, der_true, err_msg=('%s' % name))