mirror of
https://github.com/fjosw/pyerrors.git
synced 2025-03-15 06:40:24 +01:00
* [fix] Fixed intendation for flake8 * [feat] special module with Bessel functions added. * [feat] init adapted. * [tests] First test for special module added. * [tests] Check kn special function derivative explicitly vs numerical derivative of scipy function. * [feat] Added remaining autograd scipy special functions to the special module. * [feat] Imports corrected, docstring added.
12 lines
395 B
Python
12 lines
395 B
Python
import numpy as np
|
|
import scipy
|
|
import pyerrors as pe
|
|
import pytest
|
|
|
|
from autograd import jacobian
|
|
from numdifftools import Jacobian as num_jacobian
|
|
|
|
def test_kn():
|
|
for n in np.arange(0, 10):
|
|
for val in np.linspace(0.1, 7.3, 10):
|
|
assert np.isclose(num_jacobian(lambda x: scipy.special.kn(n, x))(val), jacobian(lambda x: pe.special.kn(n, x))(val), rtol=1e-10, atol=1e-10)
|