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()));
41              m2 = MatrixType::Random(rows, cols),  in product()  local
62   VERIFY_IS_APPROX((m1*m1.transpose())*m2,  m1*(m1.transpose()*m2));  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()
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()
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()
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()