Repository navigation
copy a view with non-trivial strides along with its underlying container #2803
Description
Activity
You code does not compile. I had to update it to the following:
template <class T, class... S, size_t... I> auto copy_view_to_new_container(const xt::xview<T, S...>& view, std::index_sequence<I...>) { const auto& container = view.expression(); xt::xarray<float> new_container = xt::eval(container); auto slices = view.slices(); auto new_view = xt::view(new_container, std::get<I>(slices)...); return std::make_pair(std::move(new_container), std::move(new_view)); } TEST(gpu_proxy, dtype_float_view1) { xt::xarray<double> array = xt::arange<double>(100); array.reshape({10, 10}); auto view = xt::view(array, xt::range(2, 8), xt::range(3, 7)); auto [new_container, new_view] = copy_view_to_new_container(view, std::make_index_sequence<2>()); std::cout << new_container << "\n"; std::cout << new_view << "\n"; // core dump here. }The issue you're encountering here is
new_containeris captured by reference. It's life time is bound to the code ofcopy_view_to_new_containeryou then movenew_containerthis means thatnew_viewnow has a dangling reference which when called results in a segmentation fault.I'm not really sure what you're trying to accomplish to be honest.
I have a need to copy a view along with its underlying container. Could someone show me the right way to do it?
My current code doesn't propagate the strides info to the view expression. The code crashed when I printed the new view out.
The above code crashed at xview.hpp at the function of compute_strides_impl(). When I printed out xview's m_e, I noticed that m_strides and m_backstrides are not properly initialized. It seems to be the case that the view does not have trivial_strides, so the initialization of m_strides are skipped in xview constructor.