def test_alpha(): assert_equal(1.0 / k, alpha(0.0, 1.0, 0.0, m_n, 1.0)) assert_almost_equal(1.0, (1.0 / (12.0*k)) / alpha(0.0, 1.0, 0.0, 12*m_n, 1.0)) assert_almost_equal(1.0, (1.5 / (12.0*k)) / alpha(0.5, 1.0, np.pi/2, 12*m_n, 1.0)) assert_almost_equal(1.0, (1.5 / (12.0*k*2.0)) / alpha(0.5, 1.0, np.pi/2, 12*m_n, 2.0)) assert_almost_equal(1.0, ((1.5 - 2*np.sqrt(0.5)*np.cos(np.pi/4)) / (12.0*k*2.0)) / alpha(0.5, 1.0, np.pi/4, 12*m_n, 2.0))
def test_alpha_at_theta_pi(): E_prime = np.linspace(0.5, 1.5, 101) E = np.linspace(0.75, 1.25, 101) M_A = np.linspace(1, 300, 101) T = np.linspace(1, 1800, 101) assert_array_almost_equal(alpha_at_theta_pi(E_prime, 1.0, 12*m_n, 2.0), alpha(E_prime, 1.0, np.pi, 12*m_n, 2.0)) assert_array_almost_equal(alpha_at_theta_pi(E_prime, E, 12*m_n, 2.0), alpha(E_prime, E, np.pi, 12*m_n, 2.0)) assert_array_almost_equal(alpha_at_theta_pi(E_prime, E, M_A, 2.0), alpha(E_prime, E, np.pi, M_A, 2.0)) assert_array_almost_equal(alpha_at_theta_pi(E_prime, E, M_A, T), alpha(E_prime, E, np.pi, M_A, T))
def test_alpha(): assert_equal(1.0 / k, alpha(0.0, 1.0, 0.0, m_n, 1.0)) assert_almost_equal(1.0, (1.0 / (12.0 * k)) / alpha(0.0, 1.0, 0.0, 12 * m_n, 1.0)) assert_almost_equal(1.0, (1.5 / (12.0 * k)) / alpha(0.5, 1.0, np.pi / 2, 12 * m_n, 1.0)) assert_almost_equal(1.0, (1.5 / (12.0 * k * 2.0)) / alpha(0.5, 1.0, np.pi / 2, 12 * m_n, 2.0)) assert_almost_equal(1.0, ((1.5 - 2 * np.sqrt(0.5) * np.cos(np.pi / 4)) / (12.0 * k * 2.0)) / alpha(0.5, 1.0, np.pi / 4, 12 * m_n, 2.0))
def test_alpha_at_theta_pi(): E_prime = np.linspace(0.5, 1.5, 101) E = np.linspace(0.75, 1.25, 101) M_A = np.linspace(1, 300, 101) T = np.linspace(1, 1800, 101) assert_array_almost_equal(alpha_at_theta_pi(E_prime, 1.0, 12 * m_n, 2.0), alpha(E_prime, 1.0, np.pi, 12 * m_n, 2.0)) assert_array_almost_equal(alpha_at_theta_pi(E_prime, E, 12 * m_n, 2.0), alpha(E_prime, E, np.pi, 12 * m_n, 2.0)) assert_array_almost_equal(alpha_at_theta_pi(E_prime, E, M_A, 2.0), alpha(E_prime, E, np.pi, M_A, 2.0)) assert_array_almost_equal(alpha_at_theta_pi(E_prime, E, M_A, T), alpha(E_prime, E, np.pi, M_A, T))