def METHOD_NAME():
    rng = np.random.default_rng(2)
    for _ in range(5):
        T = pt.random_transform(rng)
        cov = np.diag(rng.random((6, 6)))
        T_inv, cov_inv = pu.invert_uncertain_transform(T, cov)
        T2, cov2 = pu.invert_uncertain_transform(T_inv, cov_inv)
        assert_array_almost_equal(T, T2)
        assert_array_almost_equal(cov, cov2)