StackOverflow2013

Note that there are some explanatory texts on larger screens.

plurals

PO
text
Body
copied!<p>I do the following (array type is container/iterator range concept):</p> <pre><code>ublas::matrix<douple> A; foreach (double & element, A.data()) { } </code></pre> <p>However, this will not work for slices: your best solution is to write an iterator for them.</p> <p>Here is an example of using <code>multi_array</code> to provide storage of a custom class. Perhaps you could do the same:</p> <pre><code>template<size_t N, typename T> struct tensor_array : boost::multi_array_ref<T,N> { typedef boost::multi_array_ref<T,N> base_type; typedef T value_type; typedef T& reference; typedef const T& const_reference; tensor_array() : base_type(NULL, extents()) { // std::cout << "create" << std::endl; } template<class A> tensor_array(const A &dims) : base_type(NULL, extents()) { //std::cout << "create" << std::endl; resize(dims); } template<typename U> void resize(const U (&dims)[N]) { boost::array<U,N> dims_; std::copy(dims, dims + N, dims_.begin()); resize(dims_); } template<typename U> void resize(const boost::array<U,N> &dims) { size_t size = 1; boost::array<size_t,N> shape; for (size_t i = 0; i < N; ++i) { size *= dims[i]; shape[N-(i+1)] = dims[i]; } data_.clear(); data_.resize(size, 0); // update base_type parent set_base_ptr(&data_[0]); this->num_elements_ = size; reshape(shape); } size_t size() const { return data_.size(); } size_t size(size_t i) const { return this->shape()[N-(i+1)]; } tensor_array& fill(const T &value) { std::fill(data_.begin(), data_.end(), value); return *this; } private: typedef boost::detail::multi_array::extent_gen<N> extents; std::vector<T> data_; }; </code></pre>

Querying!

Guidance

An individual column

Larger individual text columns get their own page to allow for proper reading.

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload