Beispiel #1
0
    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))
Beispiel #2
0
    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))