def METHOD_NAME():
    """ """
    problem = {
        "bounds": [[0.0, 2.5], [0.0, 1.0], [0.0, 1.0], [0.0, 1.0]],
        "num_vars": 4,
        "names": ["x1", "x2", "x3", "x4"],
    }
    X = np.array(
        [
            [2.5, 1.0, 1.0, 1.0],
            [2.5, 0, 1.0, 0],
            [2.5, 1.0, 0, 0],
            [2.5, 0, 0, 1.0],
            [0, 0, 0, 0],
            [0, 1.0, 0, 1.0],
            [0, 0, 1.0, 1.0],
            [0, 1.0, 1.0, 0],
        ],
        dtype=float,
    )
    Y = X[:, 0] + (0.1 * X[:, 1]) + ((1.2 * X[:, 2]) * (1.3 + X[:, 3]))
    ie_names, ie = interactions(problem, Y)
    actual = ie
    assert_allclose(actual, [0.3, 0, 0, 0, 0, 0.3], rtol=1e-4, atol=1e-4)