Я пытаюсь реализовать трехмерный тензор, используя матрицу uBlas из boost в качестве бэкэнда. Одна из функций - получить ссылку на срез и позволяет легко назначить матрицу.
Ниже приведен фрагмент тензорного класса:
template<class L, class M>
class tensor {
public:
typedef L layout_type;
typedef std::size_t size_type;
typedef M array_type;
private:
size_type size1_;
size_type size2_;
size_type size3_;
array_type data_;
public:
/**
* Return a constant reference to the internal storage of a dense tensor, i.e. the raw data
* It's type depends on the type used by the tensor to store its data
*/
BOOST_UBLAS_INLINE
const array_type &data() const {
return data_;
}
/**
* Return a reference to the internal storage of a dense tensor, i.e. the raw data
* It's type depends on the type used by the tensor to store its data
*/
BOOST_UBLAS_INLINE
array_type &data() {
return data_;
}
/**@}*/
/**
* @name Slices' access
* Accessors for slices across each dimension
*/
/**@{*/
BOOST_UBLAS_INLINE
typename layout_type::template slice<1>::matrix_slice_type
at_dim1_slice(uint32_t i) {
return ublas::trans(ublas::project(data(), layout_type::template slice<1>::coord1(i, size1_, size2_, size3_),
layout_type::template slice<1>::coord2(i, size1_, size2_, size3_)));
}
}
layout_type выглядит так:
template<class M>
struct basic_dim2_major {
typedef M matrix_type;
template<int DIM, class DUMMY = void>
struct slice {
};
template<class DUMMY> struct slice<1, DUMMY> {
struct trans {
template<class T>
auto operator()(T &x) ->decltype(ublas::trans(x))
{
return ublas::trans(x);
}
};
typedef ublas::matrix_slice<matrix_type> matrix_slice_type;
typedef typename std::result_of<trans(matrix_slice_type&)>::type Type;
static BOOST_UBLAS_INLINE
Type
orient(matrix_slice_type &data){
return ublas::trans(data);
}
static BOOST_UBLAS_INLINE
ublas::slice coord1(size_type i, size_type size_i, size_type size_j, size_type size_k) {
return ublas::slice(i, size_i, size_k);
}
static BOOST_UBLAS_INLINE
ublas::slice coord2(size_type i, size_type size_i, size_type size_j, size_type size_k) {
return ublas::slice(0, 1, size_j);
}
}
}
И сценарий использования, как показано ниже:
ublas::matrix slice1(3,4);
tensor<> t(2,3,4);
t.at_dim1_slice(0) = slice1;
Проблема существует в этой строке:
return ublas::trans(ublas::project(data(), layout_type::template slice<1>::coord1(i, size1_, size2_, size3_),
layout_type::template slice<1>::coord2(i, size1_, size2_, size3_)));
Когда функции trans и project используются вместе, компилятор выбирает const overload проекта и trans, и я не могу выполнить присваивание, как описано выше. Однако, если я оставляю только проект, используется неконстантный метод, и все работает. К сожалению, транспонирование среза необходимо из-за спроектированной схемы хранения (сопоставления с двумерной матрицей).
const matrix_slice<const M> project (const M &data, const typename matrix_slice<M>::slice_type &s1, const typename matrix_slice<M>::slice_type &s2);
matrix_slice<M> project (matrix_slice<M> &data, const typename matrix_slice<M>::slice_type &s1, const typename matrix_slice<M>::slice_type &s2);
Есть ли какое-либо решение для указания правильной перегрузки функции? Или я где-то допустил ошибку?