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)
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)