test_form_hlrf_correlation.py 文件源码

def test_form_hlrf_correlation2():
    def limit_state(x1, x2, x3):
        """From choi 2007 p. 224
        x1 = Ma
        x2 = P1
        x3 = P2
        """
        return x1 - x2 - 2*x3

    X = StochasticModel(['norm', 50, 5],
                        ['norm', 10, 2],
                        ['norm', 15, 3])

    x, beta, i = form_hlrf_correlation(limit_state, X, tol=1e-5)
    pf = stats.norm.cdf(-beta)
    assert pytest.approx(pf, rel=1e-2) == 0.1073