| 69 | } |
| 70 | |
| 71 | template<typename MatrixType> void diagonal_assert(const MatrixType& m) { |
| 72 | Index rows = m.rows(); |
| 73 | Index cols = m.cols(); |
| 74 | |
| 75 | MatrixType m1 = MatrixType::Random(rows, cols); |
| 76 | |
| 77 | if (rows>=2 && cols>=2) |
| 78 | { |
| 79 | VERIFY_RAISES_ASSERT( m1 += m1.diagonal() ); |
| 80 | VERIFY_RAISES_ASSERT( m1 -= m1.diagonal() ); |
| 81 | VERIFY_RAISES_ASSERT( m1.array() *= m1.diagonal().array() ); |
| 82 | VERIFY_RAISES_ASSERT( m1.array() /= m1.diagonal().array() ); |
| 83 | } |
| 84 | } |
| 85 | |
| 86 | void test_diagonal() |
| 87 | { |