def test_degenerate_pass(self):
""" Test degenerate call. """
haz = Hazard('TC')
haz.read_mat(HAZ_TEST_MAT)
return_period = np.array([25, 50, 100, 250])
haz.intensity = sparse.csr.csr_matrix(np.zeros(haz.intensity.shape))
inten_stats = haz.local_exceedance_inten(return_period)
self.assertTrue(np.array_equal(inten_stats, np.zeros((4, 100))))
def test_ref_all_pass(self):
"""Compare against reference."""
haz = Hazard('TC')
haz.read_mat(HAZ_TEST_MAT)
return_period = np.array([25, 50, 100, 250])
inten_stats = haz.local_exceedance_inten(return_period)
self.assertAlmostEqual(inten_stats[0][0], 55.424015590131290)
self.assertAlmostEqual(inten_stats[1][0], 67.221687644669998)
self.assertAlmostEqual(inten_stats[2][0], 79.019359699208721)
self.assertAlmostEqual(inten_stats[3][0], 94.615033842370963)
self.assertAlmostEqual(inten_stats[1][66], 70.608592953031405)
self.assertAlmostEqual(inten_stats[3][33], 88.510983305123631)
self.assertAlmostEqual(inten_stats[2][99], 79.717518054203623)
