| 15 | template<typename Derived> |
| 16 | template<typename OtherDerived> |
| 17 | typename internal::traits<Derived>::Scalar |
| 18 | SparseMatrixBase<Derived>::dot(const MatrixBase<OtherDerived>& other) const |
| 19 | { |
| 20 | EIGEN_STATIC_ASSERT_VECTOR_ONLY(Derived) |
| 21 | EIGEN_STATIC_ASSERT_VECTOR_ONLY(OtherDerived) |
| 22 | EIGEN_STATIC_ASSERT_SAME_VECTOR_SIZE(Derived,OtherDerived) |
| 23 | EIGEN_STATIC_ASSERT((internal::is_same<Scalar, typename OtherDerived::Scalar>::value), |
| 24 | YOU_MIXED_DIFFERENT_NUMERIC_TYPES__YOU_NEED_TO_USE_THE_CAST_METHOD_OF_MATRIXBASE_TO_CAST_NUMERIC_TYPES_EXPLICITLY) |
| 25 | |
| 26 | eigen_assert(size() == other.size()); |
| 27 | eigen_assert(other.size()>0 && "you are using a non initialized vector"); |
| 28 | |
| 29 | internal::evaluator<Derived> thisEval(derived()); |
| 30 | typename internal::evaluator<Derived>::InnerIterator i(thisEval, 0); |
| 31 | Scalar res(0); |
| 32 | while (i) |
| 33 | { |
| 34 | res += numext::conj(i.value()) * other.coeff(i.index()); |
| 35 | ++i; |
| 36 | } |
| 37 | return res; |
| 38 | } |
| 39 | |
| 40 | template<typename Derived> |
| 41 | template<typename OtherDerived> |