1 #ifndef VIENNACL_COORDINATE_MATRIX_OPERATIONS_HPP_
2 #define VIENNACL_COORDINATE_MATRIX_OPERATIONS_HPP_
47 template<
class SCALARTYPE,
unsigned int ALIGNMENT,
unsigned int VECTOR_ALIGNMENT>
48 vector_expression<const coordinate_matrix<SCALARTYPE, ALIGNMENT>,
49 const vector<SCALARTYPE, VECTOR_ALIGNMENT>,
67 template<
class SCALARTYPE,
unsigned int ALIGNMENT,
unsigned int VECTOR_ALIGNMENT>
76 viennacl::op_prod >(mat, vec);
88 template<
class TYPE,
unsigned int ALIGNMENT,
unsigned int VECTOR_ALIGNMENT>
100 unsigned int thread_num = 256;
123 template <
typename SCALARTYPE,
unsigned int ALIGNMENT>
124 template <
unsigned int MAT_ALIGNMENT>
128 viennacl::op_prod> & proxy)
131 if (proxy.rhs().handle() == this->
handle())
151 template <
typename SCALARTYPE,
unsigned int ALIGNMENT>
152 template <
unsigned int MAT_ALIGNMENT>
168 template <
typename SCALARTYPE,
unsigned int ALIGNMENT>
169 template <
unsigned int MAT_ALIGNMENT>
187 template <
typename SCALARTYPE,
unsigned int ALIGNMENT>
188 template <
unsigned int MAT_ALIGNMENT>
194 assert(proxy.get_lhs().size1() ==
size());
205 template <
typename SCALARTYPE,
unsigned int ALIGNMENT>
206 template <
unsigned int MAT_ALIGNMENT>
212 assert(proxy.get_lhs().size1() ==
size());
215 result = *
this - result;