| 152 | } |
| 153 | |
| 154 | bool IsIdent(const matrix& A) |
| 155 | { |
| 156 | unsigned int i,j,r=A.size(); |
| 157 | if (A[0].size()!=r) { throw bad_value(); } |
| 158 | for (i=0; i<r; i++) |
| 159 | { for (j=0; j<r; j++) |
| 160 | { if (i!=j && A[i][j]!=0) { return false; } } |
| 161 | if (A[i][i]!=1) { return false; } |
| 162 | } |
| 163 | return true; |
| 164 | } |
| 165 | |
| 166 | |
| 167 | void print(const matrix& S) |