Exemplo n.º 1
0
def test_mist_basic(bands=['G','B','V','J','H','K','W1','W2','W3']):
    ic = MIST_Isochrone(bands)
    _basic_ic_checks(ic)
    print('{} ({})'.format(ic.radius(1.0, 9.5, 0.0), 0.9764494078461442))
    assert np.isclose(ic.radius(1.0, 9.5, 0.0), 0.9764494078461442)
    assert np.isclose(ic.radius(1.01, 9.72, 0.02), 1.0671791635014685)
    assert np.isclose(ic.radius(1.21, 9.38, 0.11), 1.2836469028034225)
    assert np.isclose(ic.radius(0.61, 9.89, -0.22), 0.59475269177846402)
Exemplo n.º 2
0
def test_mist_basic(bands=['J']):
    ic = MIST_Isochrone(bands, version='1.0')
    ic2 = MIST_Isochrone(bands + ['TESS', 'BP', 'RP'], version='1.1')

    _basic_ic_checks(ic)
    _basic_ic_checks(ic2)

    assert np.isclose(ic.radius(1.0, 9.5, 0.0), 0.9764494078461442)
    assert np.isclose(ic.radius(1.01, 9.72, 0.02), 1.0671791635014685)
    assert np.isclose(ic.radius(1.21, 9.38, 0.11), 1.2963342261673843)
    assert np.isclose(ic.radius(0.61, 9.89, -0.22), 0.5873830516268735)

    assert np.isclose(ic2.radius(1.0, 9.5, 0.0), 0.9765234978729515)
    assert np.isclose(ic2.radius(1.01, 9.72, 0.02), 1.0671845393364638)
    assert np.isclose(ic2.radius(1.21, 9.38, 0.11), 1.2963536270911573)
    assert np.isclose(ic2.radius(0.61, 9.89, -0.22), 0.5873849015685695)