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)