在此基础上回答了稠密矩阵的特征指数满足条件，并利用访客得到了结果。这次我要检查这些值是否是，而不是大于1。

template<typename Func>
struct lambda_as_visitor_wrapper : Func {
    lambda_as_visitor_wrapper(const Func& f) : Func(f) {}
    template<typename S, typename I>
    void init(const S& v, I i, I j) { return Func::operator()(v, i, j); }
};

template<typename Mat, typename Func>
void visit_lambda(const Mat& m, const Func& f)
{
    lambda_as_visitor_wrapper<Func> visitor(f);
    m.visit(visitor);
}


...
    std::vector<int> indices;
    int th = 1;

    visit_lambda(open.data,
        [&indices, th](double v, int i, int j) {
            if (std::isnan(v))
                indices.push_back(i);
        });

虽然我不打算对所用的数据进行太多的详细讨论，但是这个的xtensor等价物

是慢得多，本征打卡大约快8倍（366浮点）。