示例#1
0
 def test_interrelations_k_bessel(self):
     n = 7
     x = 5.746161870059113
     left = np.complex(sph_kn_bessel(n, x))
     right = np.complex((-1)**(n + 1) * pi / 2 *
                        (sph_i1n_bessel(n, x) - sph_i2n_bessel(n, x)))
     assert_allclose(left, right)
示例#2
0
 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)
示例#3
0
 def test_interrelations_k_bessel(self):
     n = 7
     x = 5.746161870059113
     left = np.complex(sph_kn_bessel(n, x))
     right = np.complex((-1) ** (n + 1) * pi / 2 * (sph_i1n_bessel(n, x) - sph_i2n_bessel(n, x)))
     assert_allclose(left, right)
示例#4
0
 def test_sph_i22_bessel(self):
     x = 6.916886515975276
     desired = (1 / x + 3 / x ** 3) * np.cosh(x) - 3 / x ** 2 * np.sinh(x)
     actual = np.complex(sph_i2n_bessel(2, x))
     assert_allclose(actual, desired)
示例#5
0
 def test_sph_i21_bessel(self):
     x = 11.070885074673603
     desired = -np.cosh(x) / x ** 2 + np.sinh(x) / x
     actual = np.complex(sph_i2n_bessel(1, x))
     assert_allclose(actual, desired)
示例#6
0
 def test_sph_i20_bessel(self):
     x = 5.807556326907072
     desired = np.cosh(x) / x
     actual = np.complex(sph_i2n_bessel(0, x))
     assert_allclose(actual, desired)
示例#7
0
 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)
示例#8
0
 def test_sph_i22_bessel(self):
     x = 6.916886515975276
     desired = (1 / x + 3 / x**3) * np.cosh(x) - 3 / x**2 * np.sinh(x)
     actual = np.complex(sph_i2n_bessel(2, x))
     assert_allclose(actual, desired)
示例#9
0
 def test_sph_i21_bessel(self):
     x = 11.070885074673603
     desired = -np.cosh(x) / x**2 + np.sinh(x) / x
     actual = np.complex(sph_i2n_bessel(1, x))
     assert_allclose(actual, desired)
示例#10
0
 def test_sph_i20_bessel(self):
     x = 5.807556326907072
     desired = np.cosh(x) / x
     actual = np.complex(sph_i2n_bessel(0, x))
     assert_allclose(actual, desired)