TEST(pyarray, broadcast_xfunc)
    {
        pyarray<int> a = {5,5,5,5,5};
        a.reshape({5, 1});

        pyarray<int> b = xt::arange(25);
        b.reshape({5, 5});

        auto func = a * b;

        for (std::size_t i = 0; i < 5; ++i)
        {
            for (std::size_t j = 0; j < 5; ++j)
            {
                std::cout << func(i, j) << std::endl;
            }
        }
    }