Sparse module: * improved performance of mat*=scalar * bug fix in cwise*
diff --git a/Eigen/src/Sparse/SparseCwise.h b/Eigen/src/Sparse/SparseCwise.h index 0697439..2206883 100644 --- a/Eigen/src/Sparse/SparseCwise.h +++ b/Eigen/src/Sparse/SparseCwise.h
@@ -73,8 +73,16 @@ operator*(const SparseMatrixBase<OtherDerived> &other) const; template<typename OtherDerived> - const EIGEN_SPARSE_CWISE_BINOP_RETURN_TYPE(ei_scalar_quotient_op) - operator/(const SparseMatrixBase<OtherDerived> &other) const; + const EIGEN_SPARSE_CWISE_PRODUCT_RETURN_TYPE + operator*(const MatrixBase<OtherDerived> &other) const; + +// template<typename OtherDerived> +// const EIGEN_SPARSE_CWISE_BINOP_RETURN_TYPE(ei_scalar_quotient_op) +// operator/(const SparseMatrixBase<OtherDerived> &other) const; +// +// template<typename OtherDerived> +// const EIGEN_SPARSE_CWISE_BINOP_RETURN_TYPE(ei_scalar_quotient_op) +// operator/(const MatrixBase<OtherDerived> &other) const; template<typename OtherDerived> const EIGEN_SPARSE_CWISE_BINOP_RETURN_TYPE(ei_scalar_min_op) @@ -96,27 +104,13 @@ // const EIGEN_SPARSE_CWISE_UNOP_RETURN_TYPE(ei_scalar_sin_op) sin() const; // const EIGEN_SPARSE_CWISE_UNOP_RETURN_TYPE(ei_scalar_pow_op) pow(const Scalar& exponent) const; + template<typename OtherDerived> + inline ExpressionType& operator*=(const SparseMatrixBase<OtherDerived> &other); + +// template<typename OtherDerived> +// inline ExpressionType& operator/=(const SparseMatrixBase<OtherDerived> &other); + /* - const ScalarAddReturnType - operator+(const Scalar& scalar) const; - - friend const ScalarAddReturnType - operator+(const Scalar& scalar, const Cwise& mat) - { return mat + scalar; } - - ExpressionType& operator+=(const Scalar& scalar); - - const ScalarAddReturnType - operator-(const Scalar& scalar) const; - - ExpressionType& operator-=(const Scalar& scalar); - - template<typename OtherDerived> - inline ExpressionType& operator*=(const MatrixBase<OtherDerived> &other); - - template<typename OtherDerived> - inline ExpressionType& operator/=(const MatrixBase<OtherDerived> &other); - template<typename OtherDerived> const EIGEN_CWISE_BINOP_RETURN_TYPE(std::less) operator<(const MatrixBase<OtherDerived>& other) const;
diff --git a/Eigen/src/Sparse/SparseCwiseBinaryOp.h b/Eigen/src/Sparse/SparseCwiseBinaryOp.h index 99bbae0..bdd1901 100644 --- a/Eigen/src/Sparse/SparseCwiseBinaryOp.h +++ b/Eigen/src/Sparse/SparseCwiseBinaryOp.h
@@ -42,40 +42,6 @@ // 4 - dense op dense product dense // generic dense -// template<typename BinaryOp, typename Lhs, typename Rhs> class ei_sparse_cwise_binary_op -// : ei_no_assignment_operator -// { -// typedef CwiseBinaryOp<BinaryOp,Lhs,Rhs> CwiseBinaryXpr; -// public: -// -// typedef typename ei_traits<CwiseBinaryXpr>::LhsNested LhsNested; -// typedef typename ei_traits<CwiseBinaryXpr>::RhsNested RhsNested; -// -// EIGEN_STRONG_INLINE ei_sparse_cwise_binary_op(const Lhs& lhs, const Rhs& rhs, const BinaryOp& func = BinaryOp()) -// : m_lhs(lhs), m_rhs(rhs), m_functor(func) -// { -// EIGEN_STATIC_ASSERT((ei_functor_allows_mixing_real_and_complex<BinaryOp>::ret -// ? int(ei_is_same_type<typename Lhs::RealScalar, typename Rhs::RealScalar>::ret) -// : int(ei_is_same_type<typename Lhs::Scalar, typename Rhs::Scalar>::ret)), -// YOU_MIXED_DIFFERENT_NUMERIC_TYPES__YOU_NEED_TO_USE_THE_CAST_METHOD_OF_MATRIXBASE_TO_CAST_NUMERIC_TYPES_EXPLICITLY) -// // require the sizes to match -// EIGEN_STATIC_ASSERT_SAME_MATRIX_SIZE(Lhs, Rhs) -// ei_assert(lhs.rows() == rhs.rows() && lhs.cols() == rhs.cols()); -// } -// -// EIGEN_STRONG_INLINE int rows() const { return m_lhs.rows(); } -// EIGEN_STRONG_INLINE int cols() const { return m_lhs.cols(); } -// -// EIGEN_STRONG_INLINE const LhsNested& _lhs() const { return m_lhs; } -// EIGEN_STRONG_INLINE const RhsNested& _rhs() const { return m_rhs; } -// EIGEN_STRONG_INLINE const UnaryOp& _functor() const { return m_functor; } -// -// protected: -// const LhsNested m_lhs; -// const RhsNested m_rhs; -// const BinaryOp m_functor; -// } - template<typename BinaryOp, typename Lhs, typename Rhs> struct ei_traits<SparseCwiseBinaryOp<BinaryOp, Lhs, Rhs> > { @@ -132,103 +98,15 @@ EIGEN_STRONG_INLINE int rows() const { return m_lhs.rows(); } EIGEN_STRONG_INLINE int cols() const { return m_lhs.cols(); } - EIGEN_STRONG_INLINE const _LhsNested& _lhs() const { return m_lhs; } - EIGEN_STRONG_INLINE const _RhsNested& _rhs() const { return m_rhs; } - EIGEN_STRONG_INLINE const BinaryOp& _functor() const { return m_functor; } + EIGEN_STRONG_INLINE const _LhsNested& lhs() const { return m_lhs; } + EIGEN_STRONG_INLINE const _RhsNested& rhs() const { return m_rhs; } + EIGEN_STRONG_INLINE const BinaryOp& functor() const { return m_functor; } protected: const LhsNested m_lhs; const RhsNested m_rhs; const BinaryOp m_functor; }; -/* -template<typename BinaryOp, typename Lhs, typename Rhs> -class CwiseBinaryOp<BinaryOp, SparseMatrixBase<Lhs>, Rhs> - : public ei_sparse_cwise_binary_op<BinaryOp, Lhs, Rhs>, - public SparseMatrixBase<CwiseBinaryOp<BinaryOp, SparseMatrixBase<Lhs>, Rhs> > -{ - public: - EIGEN_GENERIC_PUBLIC_INTERFACE(CwiseBinaryOp) - EIGEN_STRONG_INLINE CwiseBinaryOp(const Lhs& lhs, const Rhs& rhs, const BinaryOp& func = BinaryOp()) - : ei_sparse_cwise_binary_op<BinaryOp, Lhs, Rhs>(lhs, rhs, func) - {} -}; -template<typename BinaryOp, typename Lhs, typename Rhs> -class CwiseBinaryOp<BinaryOp, SparseMatrixBase<Lhs>, SparseMatrixBase<Rhs> > - : public ei_sparse_cwise_binary_op<BinaryOp, Lhs, Rhs>, - public SparseMatrixBase<CwiseBinaryOp<BinaryOp, Lhs, Rhs> > -{ - public: - EIGEN_GENERIC_PUBLIC_INTERFACE(CwiseBinaryOp) - EIGEN_STRONG_INLINE CwiseBinaryOp(const Lhs& lhs, const Rhs& rhs, const BinaryOp& func = BinaryOp()) - : ei_sparse_cwise_binary_op<BinaryOp, Lhs, Rhs>(lhs, rhs, func) - {} -};*/ - -// template<typename BinaryOp, typename Lhs, typename RhsDerived> -// class CwiseBinaryOp<BinaryOp, Lhs, SparseMatrixBase<Rhs> > : ei_no_assignment_operator, -// public SparseMatrixBase<CwiseBinaryOp<BinaryOp, SparseMatrixBase<LhsDerived>, SparseMatrixBase<Rhs> > > -// { -// public: -// -// EIGEN_GENERIC_PUBLIC_INTERFACE(CwiseBinaryOp) -// typedef typename ei_traits<CwiseBinaryOp>::LhsNested LhsNested; -// typedef typename ei_traits<CwiseBinaryOp>::RhsNested RhsNested; -// -// EIGEN_STRONG_INLINE CwiseBinaryOp(const Lhs& lhs, const Rhs& rhs, const BinaryOp& func = BinaryOp()) -// : m_lhs(lhs), m_rhs(rhs), m_functor(func) -// { -// EIGEN_STATIC_ASSERT((ei_functor_allows_mixing_real_and_complex<BinaryOp>::ret -// ? int(ei_is_same_type<typename Lhs::RealScalar, typename Rhs::RealScalar>::ret) -// : int(ei_is_same_type<typename Lhs::Scalar, typename Rhs::Scalar>::ret)), -// YOU_MIXED_DIFFERENT_NUMERIC_TYPES__YOU_NEED_TO_USE_THE_CAST_METHOD_OF_MATRIXBASE_TO_CAST_NUMERIC_TYPES_EXPLICITLY) -// // require the sizes to match -// EIGEN_STATIC_ASSERT_SAME_MATRIX_SIZE(Lhs, Rhs) -// ei_assert(lhs.rows() == rhs.rows() && lhs.cols() == rhs.cols()); -// } -// -// EIGEN_STRONG_INLINE int rows() const { return m_lhs.rows(); } -// EIGEN_STRONG_INLINE int cols() const { return m_lhs.cols(); } -// -// protected: -// const LhsNested m_lhs; -// const RhsNested m_rhs; -// const BinaryOp m_functor; -// }; -// -// template<typename BinaryOp, typename LhsDerived, typename RhsDerived> -// class CwiseBinaryOp<BinaryOp, SparseMatrixBase<LhsDerived>, SparseMatrixBase<Rhs> > : ei_no_assignment_operator, -// public SparseMatrixBase<CwiseBinaryOp<BinaryOp, SparseMatrixBase<LhsDerived>, SparseMatrixBase<Rhs> > > -// { -// public: -// -// EIGEN_GENERIC_PUBLIC_INTERFACE(CwiseBinaryOp) -// typedef typename ei_traits<CwiseBinaryOp>::LhsNested LhsNested; -// typedef typename ei_traits<CwiseBinaryOp>::RhsNested RhsNested; -// -// EIGEN_STRONG_INLINE CwiseBinaryOp(const Lhs& lhs, const Rhs& rhs, const BinaryOp& func = BinaryOp()) -// : m_lhs(lhs), m_rhs(rhs), m_functor(func) -// { -// EIGEN_STATIC_ASSERT((ei_functor_allows_mixing_real_and_complex<BinaryOp>::ret -// ? int(ei_is_same_type<typename Lhs::RealScalar, typename Rhs::RealScalar>::ret) -// : int(ei_is_same_type<typename Lhs::Scalar, typename Rhs::Scalar>::ret)), -// YOU_MIXED_DIFFERENT_NUMERIC_TYPES__YOU_NEED_TO_USE_THE_CAST_METHOD_OF_MATRIXBASE_TO_CAST_NUMERIC_TYPES_EXPLICITLY) -// // require the sizes to match -// EIGEN_STATIC_ASSERT_SAME_MATRIX_SIZE(Lhs, Rhs) -// ei_assert(lhs.rows() == rhs.rows() && lhs.cols() == rhs.cols()); -// } -// -// EIGEN_STRONG_INLINE int rows() const { return m_lhs.rows(); } -// EIGEN_STRONG_INLINE int cols() const { return m_lhs.cols(); } -// -// protected: -// const LhsNested m_lhs; -// const RhsNested m_rhs; -// const BinaryOp m_functor; -// }; - -// template<typename T> struct ei_is_scalar_product { enum { ret = false }; }; -// template<typename T> struct ei_is_scalar_product<ei_scalar_product_op<T> > { enum { ret = true }; }; template<typename BinaryOp, typename Lhs, typename Rhs, typename Derived, int _LhsStorageMode = int(Lhs::Flags) & SparseBit, @@ -276,7 +154,7 @@ public: EIGEN_STRONG_INLINE ei_sparse_cwise_binary_op_inner_iterator_selector(const CwiseBinaryXpr& xpr, int outer) - : m_lhsIter(xpr._lhs(),outer), m_rhsIter(xpr._rhs()), m_functor(xpr._functor()) + : m_lhsIter(xpr.lhs(),outer), m_rhsIter(xpr.rhs(),outer), m_functor(xpr.functor()) { this->operator++(); } @@ -323,24 +201,6 @@ int m_id; }; -// template<typename BinaryOp, typename Lhs, typename Rhs> -// class SparseCwiseBinaryOp<BinaryOp, Lhs, Rhs>::InnerIterator -// : public ei_sparse_cwise_binary_op_inner_iterator_selector<BinaryOp,Lhs,Rhs, typename SparseCwiseBinaryOp<BinaryOp,Lhs,Rhs>::InnerIterator> -// { -// typedef CwiseBinaryOpInnerIterator< -// BinaryOp,Lhs,Rhs, typename CwiseBinaryOp<BinaryOp,Lhs,Rhs>::InnerIterator> Base; -// public: -// typedef typename CwiseBinaryOp::Scalar Scalar; -// typedef typename ei_traits<CwiseBinaryOp>::_LhsNested _LhsNested; -// typedef typename _LhsNested::InnerIterator LhsIterator; -// typedef typename ei_traits<CwiseBinaryOp>::_RhsNested _RhsNested; -// typedef typename _RhsNested::InnerIterator RhsIterator; -// // public: -// EIGEN_STRONG_INLINE InnerIterator(const CwiseBinaryOp& binOp, int outer) -// : Base(binOp.m_lhs,binOp.m_rhs,binOp.m_functor,outer) -// {} -// }; - // sparse - sparse (product) template<typename T, typename Lhs, typename Rhs, typename Derived> class ei_sparse_cwise_binary_op_inner_iterator_selector<ei_scalar_product_op<T>, Lhs, Rhs, Derived, IsSparse, IsSparse> @@ -355,9 +215,9 @@ public: EIGEN_STRONG_INLINE ei_sparse_cwise_binary_op_inner_iterator_selector(const CwiseBinaryXpr& xpr, int outer) - : m_lhsIter(xpr._lhs(),outer), m_rhsIter(xpr._rhs()), m_functor(xpr._functor()) + : m_lhsIter(xpr.lhs(),outer), m_rhsIter(xpr.rhs(),outer), m_functor(xpr.functor()) { - while (m_lhsIter && m_rhsIter && m_lhsIter.index() != m_rhsIter.index()) + while (m_lhsIter && m_rhsIter && (m_lhsIter.index() != m_rhsIter.index())) { if (m_lhsIter.index() < m_rhsIter.index()) ++m_lhsIter; @@ -368,15 +228,11 @@ EIGEN_STRONG_INLINE Derived& operator++() { - while (m_lhsIter && m_rhsIter) + ++m_lhsIter; + ++m_rhsIter; + while (m_lhsIter && m_rhsIter && (m_lhsIter.index() != m_rhsIter.index())) { - if (m_lhsIter.index() == m_rhsIter.index()) - { - ++m_lhsIter; - ++m_rhsIter; - break; - } - else if (m_lhsIter.index() < m_rhsIter.index()) + if (m_lhsIter.index() < m_rhsIter.index()) ++m_lhsIter; else ++m_rhsIter; @@ -388,7 +244,7 @@ EIGEN_STRONG_INLINE int index() const { return m_lhsIter.index(); } - EIGEN_STRONG_INLINE operator bool() const { return m_lhsIter && m_rhsIter; } + EIGEN_STRONG_INLINE operator bool() const { return (m_lhsIter && m_rhsIter); } protected: LhsIterator m_lhsIter; @@ -419,7 +275,8 @@ } EIGEN_STRONG_INLINE Scalar value() const - { return m_functor(m_lhsIter.value(), m_xpr.rhs().coeff(IsRowMajor?m_outer:m_lhsIter.index(),IsRowMajor?m_lhsIter.index():m_outer)); } + { return m_functor(m_lhsIter.value(), + m_xpr.rhs().coeff(IsRowMajor?m_outer:m_lhsIter.index(),IsRowMajor?m_lhsIter.index():m_outer)); } EIGEN_STRONG_INLINE int index() const { return m_lhsIter.index(); } @@ -513,22 +370,38 @@ template<typename ExpressionType> template<typename OtherDerived> -EIGEN_STRONG_INLINE const EIGEN_SPARSE_CWISE_BINOP_RETURN_TYPE(ei_scalar_quotient_op) -SparseCwise<ExpressionType>::operator/(const SparseMatrixBase<OtherDerived> &other) const +EIGEN_STRONG_INLINE const EIGEN_SPARSE_CWISE_PRODUCT_RETURN_TYPE +SparseCwise<ExpressionType>::operator*(const MatrixBase<OtherDerived> &other) const { - return EIGEN_SPARSE_CWISE_BINOP_RETURN_TYPE(ei_scalar_quotient_op)(_expression(), other.derived()); + return EIGEN_SPARSE_CWISE_PRODUCT_RETURN_TYPE(_expression(), other.derived()); } // template<typename ExpressionType> // template<typename OtherDerived> -// inline ExpressionType& Cwise<ExpressionType>::operator*=(const SparseMatrixBase<OtherDerived> &other) +// EIGEN_STRONG_INLINE const EIGEN_SPARSE_CWISE_BINOP_RETURN_TYPE(ei_scalar_quotient_op) +// SparseCwise<ExpressionType>::operator/(const SparseMatrixBase<OtherDerived> &other) const // { -// return m_matrix.const_cast_derived() = *this * other; +// return EIGEN_SPARSE_CWISE_BINOP_RETURN_TYPE(ei_scalar_quotient_op)(_expression(), other.derived()); // } +// +// template<typename ExpressionType> +// template<typename OtherDerived> +// EIGEN_STRONG_INLINE const EIGEN_SPARSE_CWISE_BINOP_RETURN_TYPE(ei_scalar_quotient_op) +// SparseCwise<ExpressionType>::operator/(const MatrixBase<OtherDerived> &other) const +// { +// return EIGEN_SPARSE_CWISE_BINOP_RETURN_TYPE(ei_scalar_quotient_op)(_expression(), other.derived()); +// } + +template<typename ExpressionType> +template<typename OtherDerived> +inline ExpressionType& SparseCwise<ExpressionType>::operator*=(const SparseMatrixBase<OtherDerived> &other) +{ + return m_matrix.const_cast_derived() = *this * other; +} // template<typename ExpressionType> // template<typename OtherDerived> -// inline ExpressionType& Cwise<ExpressionType>::operator/=(const SparseMatrixBase<OtherDerived> &other) +// inline ExpressionType& SparseCwise<ExpressionType>::operator/=(const SparseMatrixBase<OtherDerived> &other) // { // return m_matrix.const_cast_derived() = *this / other; // }
diff --git a/Eigen/src/Sparse/SparseCwiseUnaryOp.h b/Eigen/src/Sparse/SparseCwiseUnaryOp.h index 2b2050e..e66baa0 100644 --- a/Eigen/src/Sparse/SparseCwiseUnaryOp.h +++ b/Eigen/src/Sparse/SparseCwiseUnaryOp.h
@@ -162,14 +162,20 @@ EIGEN_STRONG_INLINE Derived& SparseMatrixBase<Derived>::operator*=(const Scalar& other) { - return *this = *this * other; + for (int j=0; j<outerSize(); ++j) + for (typename Derived::InnerIterator i(derived(),j); i; ++i) + i.valueRef() *= other; + return derived(); } template<typename Derived> EIGEN_STRONG_INLINE Derived& SparseMatrixBase<Derived>::operator/=(const Scalar& other) { - return *this = *this / other; + for (int j=0; j<outerSize(); ++j) + for (typename Derived::InnerIterator i(derived(),j); i; ++i) + i.valueRef() /= other; + return derived(); } #endif // EIGEN_SPARSE_CWISE_UNARY_OP_H
diff --git a/Eigen/src/Sparse/SparseMatrix.h b/Eigen/src/Sparse/SparseMatrix.h index 07fc0be..a5d91f1 100644 --- a/Eigen/src/Sparse/SparseMatrix.h +++ b/Eigen/src/Sparse/SparseMatrix.h
@@ -419,6 +419,7 @@ inline InnerIterator& operator++() { m_id++; return *this; } inline Scalar value() const { return m_matrix.m_data.value(m_id); } + inline Scalar& valueRef() { return const_cast<Scalar&>(m_matrix.m_data.value(m_id)); } inline int index() const { return m_matrix.m_data.index(m_id); }
diff --git a/Eigen/src/Sparse/SparseVector.h b/Eigen/src/Sparse/SparseVector.h index b189cb5..49f5bde 100644 --- a/Eigen/src/Sparse/SparseVector.h +++ b/Eigen/src/Sparse/SparseVector.h
@@ -316,6 +316,7 @@ inline InnerIterator& operator++() { m_id++; return *this; } inline Scalar value() const { return m_vector.m_data.value(m_id); } + inline Scalar& valueRef() { return const_cast<Scalar&>(m_vector.m_data.value(m_id)); } inline int index() const { return m_vector.m_data.index(m_id); }
diff --git a/test/sparse_basic.cpp b/test/sparse_basic.cpp index 07a38dd..84cf4cd 100644 --- a/test/sparse_basic.cpp +++ b/test/sparse_basic.cpp
@@ -52,6 +52,7 @@ SparseMatrix<Scalar> m(rows, cols); DenseMatrix refMat = DenseMatrix::Zero(rows, cols); DenseVector vec1 = DenseVector::Random(rows); + Scalar s1 = ei_random<Scalar>(); std::vector<Vector2i> zeroCoords; std::vector<Vector2i> nonzeroCoords; @@ -191,6 +192,35 @@ // std::cerr << m.transpose() << "\n\n" << refMat.transpose() << "\n\n"; // VERIFY_IS_APPROX(m, refMat); + // test basic computations + { + DenseMatrix refM1 = DenseMatrix::Zero(rows, rows); + DenseMatrix refM2 = DenseMatrix::Zero(rows, rows); + DenseMatrix refM3 = DenseMatrix::Zero(rows, rows); + DenseMatrix refM4 = DenseMatrix::Zero(rows, rows); + SparseMatrix<Scalar> m1(rows, rows); + SparseMatrix<Scalar> m2(rows, rows); + SparseMatrix<Scalar> m3(rows, rows); + SparseMatrix<Scalar> m4(rows, rows); + initSparse<Scalar>(density, refM1, m1); + initSparse<Scalar>(density, refM2, m2); + initSparse<Scalar>(density, refM3, m3); + initSparse<Scalar>(density, refM4, m4); + + VERIFY_IS_APPROX(m1+m2, refM1+refM2); + VERIFY_IS_APPROX(m1+m2+m3, refM1+refM2+refM3); + VERIFY_IS_APPROX(m3.cwise()*(m1+m2), refM3.cwise()*(refM1+refM2)); + VERIFY_IS_APPROX(m1*s1-m2, refM1*s1-refM2); + + VERIFY_IS_APPROX(m1*=s1, refM1*=s1); + VERIFY_IS_APPROX(m1/=s1, refM1/=s1); + + refM4.setRandom(); + // sparse cwise* dense + VERIFY_IS_APPROX(m3.cwise()*refM4, refM3.cwise()*refM4); +// VERIFY_IS_APPROX(m3.cwise()/refM4, refM3.cwise()/refM4); + } + // test innerVector() { DenseMatrix refMat2 = DenseMatrix::Zero(rows, rows); @@ -198,8 +228,8 @@ initSparse<Scalar>(density, refMat2, m2); int j0 = ei_random(0,rows-1); int j1 = ei_random(0,rows-1); -// VERIFY_IS_APPROX(m2.innerVector(j0), refMat2.col(j0)); -// VERIFY_IS_APPROX(m2.innerVector(j0)+m2.innerVector(j1), refMat2.col(j0)+refMat2.col(j1)); + VERIFY_IS_APPROX(m2.innerVector(j0), refMat2.col(j0)); + VERIFY_IS_APPROX(m2.innerVector(j0)+m2.innerVector(j1), refMat2.col(j0)+refMat2.col(j1)); } // test transpose @@ -227,13 +257,13 @@ VERIFY_IS_APPROX(m4=m2.transpose()*m3, refMat4=refMat2.transpose()*refMat3); VERIFY_IS_APPROX(m4=m2.transpose()*m3.transpose(), refMat4=refMat2.transpose()*refMat3.transpose()); VERIFY_IS_APPROX(m4=m2*m3.transpose(), refMat4=refMat2*refMat3.transpose()); - + // sparse * dense VERIFY_IS_APPROX(dm4=m2*refMat3, refMat4=refMat2*refMat3); VERIFY_IS_APPROX(dm4=m2*refMat3.transpose(), refMat4=refMat2*refMat3.transpose()); VERIFY_IS_APPROX(dm4=m2.transpose()*refMat3, refMat4=refMat2.transpose()*refMat3); VERIFY_IS_APPROX(dm4=m2.transpose()*refMat3.transpose(), refMat4=refMat2.transpose()*refMat3.transpose()); - + // dense * sparse VERIFY_IS_APPROX(dm4=refMat2*m3, refMat4=refMat2*refMat3); VERIFY_IS_APPROX(dm4=refMat2*m3.transpose(), refMat4=refMat2*refMat3.transpose());
diff --git a/test/sparse_vector.cpp b/test/sparse_vector.cpp index 18ef748..0a66af6 100644 --- a/test/sparse_vector.cpp +++ b/test/sparse_vector.cpp
@@ -77,7 +77,9 @@ VERIFY_IS_APPROX(v1*s1-v2, refV1*s1-refV2); -// std::cerr << v1.dot(v2) << " == " << refV1.dot(refV2) << "\n"; + VERIFY_IS_APPROX(v1*=s1, refV1*=s1); + VERIFY_IS_APPROX(v1/=s1, refV1/=s1); + VERIFY_IS_APPROX(v1.dot(v2), refV1.dot(refV2)); }