def test_fd_derivative(): x = np.linspace(-1, 1, 25) fx = np.exp(x) for n in range(1, 7): df = fd_derivative(fx, x, n=n) m = n // 2 + 2 np.testing.assert_allclose(df[m:-m], fx[m:-m], atol=1e-5) np.testing.assert_allclose(df[-m:], fx[-m:], atol=1e-4) np.testing.assert_allclose(df[:m], fx[:m], atol=1e-4)
def test_fd_derivative(): x = np.linspace(-1, 1, 25) h = np.diff(x).mean() fx = np.exp(x) for n in range(1, 7): df = fd_derivative(fx, x, n=n) m = n // 2 + 2 np.testing.assert_allclose(df[m:-m], fx[m:-m], atol=1e-5) np.testing.assert_allclose(df[-m:], fx[-m:], atol=1e-4) np.testing.assert_allclose(df[:m], fx[:m], atol=1e-4)