def test_computes_the_rates_of_exceedance(self): expected_rates = numpy.array([0.124, 0.0344, 0.0189, 0.0133, 0.0078, 0.0078, 0.0056, 0.0044, 0.0044, 0.0044, 0.0011, 0.0011, 0.0011, 0.0011, 0.011, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]) self.assertTrue(numpy.allclose(expected_rates, prob.compute_rates_of_exceedance( self.cum_histogram, self.gmfs["TSES"]), atol=0.01))
def test_computes_probs_of_exceedance(self): expected_probs = [0.9980, 0.8213, 0.6111, 0.4866, 0.3222, 0.3222, 0.2425, 0.1993, 0.1993, 0.1993, 0.0540, 0.0540, 0.0540, 0.0540, 0.0540, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0] self.assertTrue(numpy.allclose(expected_probs, prob.compute_probs_of_exceedance( prob.compute_rates_of_exceedance( self.cum_histogram, self.gmfs["TSES"]), self.gmfs["TimeSpan"]), atol=0.0001))