1 #ifndef VIENNACL_LINALG_HANKEL_MATRIX_OPERATIONS_HPP_
2 #define VIENNACL_LINALG_HANKEL_MATRIX_OPERATIONS_HPP_
49 template<
class SCALARTYPE,
unsigned int ALIGNMENT,
unsigned int VECTOR_ALIGNMENT>
50 vector_expression<const hankel_matrix<SCALARTYPE, ALIGNMENT>,
51 const vector<SCALARTYPE, VECTOR_ALIGNMENT>,
69 template<
class SCALARTYPE,
unsigned int ALIGNMENT,
unsigned int VECTOR_ALIGNMENT>
78 viennacl::op_prod >(mat, vec);
89 template<
class SCALARTYPE,
unsigned int ALIGNMENT,
unsigned int VECTOR_ALIGNMENT>
98 viennacl::detail::fft::reverse(result);
109 template <
typename SCALARTYPE,
unsigned int ALIGNMENT>
110 template <
unsigned int MAT_ALIGNMENT>
114 viennacl::op_prod> & proxy)
117 if (proxy.rhs().handle() == this->
handle())
137 template <
typename SCALARTYPE,
unsigned int ALIGNMENT>
138 template <
unsigned int MAT_ALIGNMENT>
154 template <
typename SCALARTYPE,
unsigned int ALIGNMENT>
155 template <
unsigned int MAT_ALIGNMENT>
173 template <
typename SCALARTYPE,
unsigned int ALIGNMENT>
174 template <
unsigned int MAT_ALIGNMENT>
180 assert(proxy.get_lhs().size1() ==
size());
191 template <
typename SCALARTYPE,
unsigned int ALIGNMENT>
192 template <
unsigned int MAT_ALIGNMENT>
198 assert(proxy.get_lhs().size1() ==
size());
201 result = *
this - result;