Python CUBAの例

コード例 #1

ファイル: test_cuba.py プロジェクト: nhmc/pyigm

def test_zinterp():
    # 
    cuba = CUBA()
    # 
    jnu = cuba.zinterp_jnu(2.1)
    np.testing.assert_allclose(jnu[0:5].value, 
        [1.36534694e-19, 1.33934694e-19, 1.31393878e-19, 1.28893878e-19, 1.26493878e-19])
    assert(jnu.unit == u.erg/u.s/u.cm**2)

コード例 #2

def test_phi():
    cuba = CUBA()
    phi = cuba.phi(2.1)
    #
    np.testing.assert_allclose(phi.value, 381837.59540310985)
    assert (phi.unit == 1 / u.s / u.cm**2)
    #
    phi = cuba.phi(2.1, min_energy=50. * u.eV)
    np.testing.assert_allclose(phi.value, 4 * 8171.47189724)

コード例 #3

ファイル: test_cuba.py プロジェクト: nhmc/pyigm

def test_phi():
    cuba = CUBA()
    phi = cuba.phi(2.1)
    #    
    np.testing.assert_allclose(phi.value, 95459.3988508)
    assert(phi.unit == 1/u.s/u.cm**2)
    #
    phi = cuba.phi(2.1, min_energy=50.*u.eV)
    np.testing.assert_allclose(phi.value, 8171.47189724)

コード例 #4

def test_zinterp():
    #
    cuba = CUBA()
    #
    jnu = cuba.zinterp_jnu(2.1)
    np.testing.assert_allclose(jnu[0:5].value, [
        1.36534694e-19, 1.33934694e-19, 1.31393878e-19, 1.28893878e-19,
        1.26493878e-19
    ])
    assert (jnu.unit == u.erg / u.s / u.cm**2)

コード例 #5

def test_logU():
    cuba = CUBA()
    logU = cuba.logU(0.2)
    #
    np.testing.assert_allclose(logU, -6.13089862882)

コード例 #6

def test_init():
    # Class
    cuba = CUBA()
    # z values
    np.testing.assert_allclose(cuba.z[0:5],
                               [0., 0.04912, 0.1006, 0.1547, 0.2114])