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))