def test_interrelations_i2_bessel(self):
n = 9
x = 6.083765036260925
left = np.complex(sph_i2n_bessel(n, x))
right = np.complex(mpc(0, 1) ** (-n - 1) * sph_yn_bessel(n, 1j * x))
assert_allclose(left, right)
def test_interrelations_h1_bessel(self):
n = 6
x = 8.160875740148962
left = np.complex(sph_h1n_bessel(n, x))
right = np.complex(sph_jn_bessel(n, x) + 1j * sph_yn_bessel(n, x))
assert_allclose(left, right)
def test_interrelations_h2_bessel(self):
n = 5
x = 2.826012588052626
left = np.complex(sph_h2n_bessel(n, x))
right = np.complex(sph_jn_bessel(n, x) - 1j * sph_yn_bessel(n, x))
assert_allclose(left, right)
def test_sph_y1_bessel(self):
x = 29.2
desired = -np.cos(x) / x ** 2 - np.sin(x) / x
actual = np.complex(sph_yn_bessel(1, x))
assert_allclose(actual, desired)
def test_sph_y2_bessel(self):
x = 0.321
desired = (1.0 / x - 3 / x ** 3) * np.cos(x) - 3 / x ** 2 * np.sin(x)
actual = np.complex(sph_yn_bessel(2, x))
assert_allclose(actual, desired)
def test_interrelations_i2_bessel(self):
n = 9
x = 6.083765036260925
left = np.complex(sph_i2n_bessel(n, x))
right = np.complex(mpc(0, 1)**(-n - 1) * sph_yn_bessel(n, 1j * x))
assert_allclose(left, right)
def test_sph_y0_bessel(self):
x = 4.3
desired = -np.cos(x) / x
actual = np.complex(sph_yn_bessel(0, x))
assert_allclose(actual, desired)
def test_interrelations_h2_bessel(self):
n = 5
x = 2.826012588052626
left = np.complex(sph_h2n_bessel(n, x))
right = np.complex(sph_jn_bessel(n, x) - 1j * sph_yn_bessel(n, x))
assert_allclose(left, right)
def test_interrelations_h1_bessel(self):
n = 6
x = 8.160875740148962
left = np.complex(sph_h1n_bessel(n, x))
right = np.complex(sph_jn_bessel(n, x) + 1j * sph_yn_bessel(n, x))
assert_allclose(left, right)
def test_sph_y2_bessel(self):
x = 0.321
desired = (1.0 / x - 3 / x**3) * np.cos(x) - 3 / x**2 * np.sin(x)
actual = np.complex(sph_yn_bessel(2, x))
assert_allclose(actual, desired)
def test_sph_y1_bessel(self):
x = 29.2
desired = -np.cos(x) / x**2 - np.sin(x) / x
actual = np.complex(sph_yn_bessel(1, x))
assert_allclose(actual, desired)
def test_sph_y0_bessel(self):
x = 4.3
desired = -np.cos(x) / x
actual = np.complex(sph_yn_bessel(0, x))
assert_allclose(actual, desired)
