14 bool areNotApprox(const MatrixBase<Derived1>& m1, const MatrixBase<Derived2>& m2, typename Derived1…
16   return !((m1-m2).cwiseAbs2().maxCoeff() < epsilon * epsilon
17                           * (std::max)(m1.cwiseAbs2().maxCoeff(), m2.cwiseAbs2().maxCoeff()));
40   MatrixType m1 = MatrixType::Random(rows, cols),  in product()  local
52   OtherMajorMatrixType tm1 = m1;  in product()
62   VERIFY_IS_APPROX((m1*m1.transpose())*m2,  m1*(m1.transpose()*m2));  in product()
63   m3 = m1;  in product()
64   m3 *= m1.transpose() * m2;  in product()
65   VERIFY_IS_APPROX(m3,                      m1 * (m1.transpose()*m2));  in product()
66   VERIFY_IS_APPROX(m3,                      m1 * (m1.transpose()*m2));  in product()
69   VERIFY_IS_APPROX(square*(m1 + m2),        square*m1+square*m2);  in product()
70   VERIFY_IS_APPROX(square*(m1 - m2),        square*m1-square*m2);  in product()
73   VERIFY_IS_APPROX(s1*(square*m1),          (s1*square)*m1);  in product()
74   VERIFY_IS_APPROX(s1*(square*m1),          square*(m1*s1));  in product()
83      VERIFY_RAISES_ASSERT(m3 = m1*m1);  in product()
89     VERIFY(areNotApprox(m1.transpose()*m2,m2.transpose()*m1));  in product()
94   res.noalias() += m1 * m2.transpose();  in product()
95   VERIFY_IS_APPROX(res, square + m1 * m2.transpose());  in product()
98     VERIFY(areNotApprox(res,square + m2 * m1.transpose()));  in product()
101   vcres.noalias() += m1.transpose() * v1;  in product()
102   VERIFY_IS_APPROX(vcres, vc2 + m1.transpose() * v1);  in product()
106   res.noalias() -= m1 * m2.transpose();  in product()
107   VERIFY_IS_APPROX(res, square - (m1 * m2.transpose()));  in product()
110     VERIFY(areNotApprox(res,square - m2 * m1.transpose()));  in product()
113   vcres.noalias() -= m1.transpose() * v1;  in product()
114   VERIFY_IS_APPROX(vcres, vc2 - m1.transpose() * v1);  in product()
116   tm1 = m1;  in product()
117   VERIFY_IS_APPROX(tm1.transpose() * v1, m1.transpose() * v1);  in product()
118   VERIFY_IS_APPROX(v1.transpose() * tm1, v1.transpose() * m1);  in product()
122     res.row(i) = m1.row(i) * m2.transpose();  in product()
123   VERIFY_IS_APPROX(res, m1 * m2.transpose());  in product()
126     res.col(i) = m1 * m2.transpose().col(i);  in product()
127   VERIFY_IS_APPROX(res, m1 * m2.transpose());  in product()
130   res2.noalias() += m1.transpose() * m2;  in product()
131   VERIFY_IS_APPROX(res2, square2 + m1.transpose() * m2);  in product()
134     VERIFY(areNotApprox(res2,square2 + m2.transpose() * m1));  in product()