def test_Chelyabinsk(): final_state = run.run_asteroid(planet="Earth", asteroid="Chelyabinsk", show=False, anal_assumption=False, tol=1e-8) KE = final_state[6] z = final_state[4] assert np.all(abs(np.diff(KE) / np.diff(z)) * 1e3 / 4.184E12 <= 10**4)
def test_dist2(): # assert x is increasing final_state = run.run_asteroid(planet="Earth", asteroid="Chelyabinsk", show=False, anal_assumption=False, tol=1e-8) x = final_state[5] assert all(np.diff(x) >= 0)
def test_KE(): # assert KE is bounded final_state = run.run_asteroid(planet="Earth", asteroid="Chelyabinsk", show=False, anal_assumption=False, tol=1e-8) KE = final_state[6] assert np.all(KE / (4.184E12) <= 10**4)
def test_speed1(): # assert that v is decreasing final_state = run.run_asteroid(planet="Earth", asteroid="Chelyabinsk", show=False, anal_assumption=False, tol=1e-8) v = final_state[0] assert all(np.around(np.diff(v))) <= 0
def test_alt1(): # assert that altitude is bounded final_state = run.run_asteroid(planet="Earth", asteroid="Chelyabinsk", show=False, anal_assumption=False, tol=1e-8) z = final_state[4] assert 0 <= all(z) <= 1E5
def test_mass(): # assert that m is bounded final_state = run.run_asteroid(planet="Earth", asteroid="Chelyabinsk", show=False, anal_assumption=False, tol=1e-8) m = final_state[1] assert all(0 < m) and all(m <= 1E9)
def test_speed2(): # assert that v is bounded final_state = run.run_asteroid(planet="Earth", asteroid="Chelyabinsk", show=False, anal_assumption=False, tol=1e-8) v = final_state[0] assert all(v >= 0) and all(v <= 100E3)
def test_alt2(): # assert that altitude is decreasing final_state = eroscode.main() final_state = run.run_asteroid(planet="Earth", asteroid="Chelyabinsk", show=False, anal_assumption=False, tol=1e-8) z = final_state[4] assert all(np.diff(z) <= 0)
def test_Tunguska(): #Call Tungaska Parameters final_state = run.run_asteroid(planet="Earth", asteroid="Chelyabinsk", show=False, anal_assumption=False, tol=1e-8) KE = final_state[6] z = final_state[4] assert all(abs(np.diff(KE) / np.diff(z)) * 1e3 / 4.184E12 <= 10**4)
def test_alt3(): # assert crater event for large asteroids final_state = eroscode.main() final_state = run.run_asteroid(planet="Earth", asteroid="Chelyabinsk", show=False, anal_assumption=False, tol=1e-8) r = final_state[7] burst_index = final_state[8] if r[0] > 100: assert burst_index == 0