unit_matrix(r,c) can return an evaluated matrix object.

Cherry-picked from 052e9928.

diff --git a/ginac/matrix.cpp b/ginac/matrix.cpp
index 3c46b17efa87674b0ab5fbd4a3780ab9899066c0..82c2255fb76c30440a46d7de495b516d50940289 100644 (file)
--- a/ginac/matrix.cpp
+++ b/ginac/matrix.cpp
@@ -1593,7 +1593,7 @@ ex diag_matrix(const lst & l)
 ex unit_matrix(unsigned r, unsigned c)
 {
        matrix &Id = *new matrix(r, c);
-       Id.setflag(status_flags::dynallocated);
+       Id.setflag(status_flags::dynallocated | status_flags::evaluated);
        for (unsigned i=0; i<r && i<c; i++)
                Id(i,i) = _ex1;