38   VERIFY_IS_APPROX(sq_m1, v1.asDiagonal().toDenseMatrix());  in diagonalmatrices()
40   VERIFY_IS_APPROX(sq_m1, v1.asDiagonal().toDenseMatrix());  in diagonalmatrices()
42   VERIFY_IS_APPROX(sq_m1, sq_m2);  in diagonalmatrices()
46   VERIFY_IS_APPROX(ldm1.diagonal(), ldm3.diagonal());  in diagonalmatrices()
48   VERIFY_IS_APPROX(ldm1.diagonal(), ldm4.diagonal());  in diagonalmatrices()
51   VERIFY_IS_APPROX(sq_m1, ldm1.toDenseMatrix());  in diagonalmatrices()
53   VERIFY_IS_APPROX(sq_m1, ldm1.toDenseMatrix());  in diagonalmatrices()
58   VERIFY_IS_APPROX( ((ldm1 * m1)(i,j))  , ldm1.diagonal()(i) * m1(i,j) );  in diagonalmatrices()
59   VERIFY_IS_APPROX( ((ldm1 * (m1+m2))(i,j))  , ldm1.diagonal()(i) * (m1+m2)(i,j) );  in diagonalmatrices()
60   VERIFY_IS_APPROX( ((m1 * rdm1)(i,j))  , rdm1.diagonal()(j) * m1(i,j) );  in diagonalmatrices()
61   VERIFY_IS_APPROX( ((v1.asDiagonal() * m1)(i,j))  , v1(i) * m1(i,j) );  in diagonalmatrices()
62   VERIFY_IS_APPROX( ((m1 * rv1.asDiagonal())(i,j))  , rv1(j) * m1(i,j) );  in diagonalmatrices()
63   VERIFY_IS_APPROX( (((v1+v2).asDiagonal() * m1)(i,j))  , (v1+v2)(i) * m1(i,j) );  in diagonalmatrices()
64   VERIFY_IS_APPROX( (((v1+v2).asDiagonal() * (m1+m2))(i,j))  , (v1+v2)(i) * (m1+m2)(i,j) );  in diagonalmatrices()
65   VERIFY_IS_APPROX( ((m1 * (rv1+rv2).asDiagonal())(i,j))  , (rv1+rv2)(j) * m1(i,j) );  in diagonalmatrices()
66   VERIFY_IS_APPROX( (((m1+m2) * (rv1+rv2).asDiagonal())(i,j))  , (rv1+rv2)(j) * (m1+m2)(i,j) );  in diagonalmatrices()
74   VERIFY_IS_APPROX((big.block(i,j,rows,cols)) , v1.asDiagonal() * m1 );  in diagonalmatrices()
78   VERIFY_IS_APPROX((big.block(i,j,rows,cols)) , m1 * rv1.asDiagonal() );  in diagonalmatrices()
82   VERIFY_IS_APPROX(LeftDiagonalMatrix(ldm1*s1).diagonal(), ldm1.diagonal() * s1);  in diagonalmatrices()
83   VERIFY_IS_APPROX(LeftDiagonalMatrix(s1*ldm1).diagonal(), s1 * ldm1.diagonal());  in diagonalmatrices()
85   VERIFY_IS_APPROX(m1 * (rdm1 * s1), (m1 * rdm1) * s1);  in diagonalmatrices()
86   VERIFY_IS_APPROX(m1 * (s1 * rdm1), (m1 * rdm1) * s1);  in diagonalmatrices()