Function numpy_view

include/bind.h:93–106

#define Float typename Solver::Float
#define Vector typename Solver::Vector
std::function<py::array_t<Float>(Solver &)> numpy_view(std::shared_ptr<std::vector<Vector>> Solver::*field) {
#undef Vector
#undef Float
    const size_t dim = Solver::dim;
    typedef typename Solver::Float Float;

    return [field](Solver &solver) {
        Float *data = reinterpret_cast<Float *>((solver.*field)->data());
        py::capsule dummy(data, [](void *ptr) {});
        return py::array_t<Float>({(solver.*field)->size(), dim},
                                  {sizeof(Float) * dim, sizeof(Float)},
                                  data, dummy);
    };
}

// Class interface
template<class Index = size_t>