- #1
oki
- 5
- 0
Matrix diagonalization
Hi, guys!
I've graduated some time ago, and by now have completely fogotten all those matrix rules. Could you please advise me. The problem is following:
I have a large square matrix (like 1000x1000) and I have to rise this matrix to a large power, like 1000. The matrix is quite a sparse one. It looks like this:
x0x00000000000000000
0x000000000000000000
000x00x0x0x0x0x0x0x0x
0xx0xx00000000000000
00x00000000000000000
000000x0000000000000
00x00000000000000000
00000000x00000000000
00x00000000000000000
0000000000x000000000
00x00000000000000000
000000000000x0000000
00x00000000000000000
00000000000000x00000
all x's are different, all zero's are just zero's.
So, my first quess would be that I have to diagonalize the matrix somehow. But I don't remember how to do this. I want to try manually first. Then if it fails - numerically (I prefer fortran).
Any suggestions on how to diagonalized? The right bottom part of the matrix looks like it is already diagonal (although the elements are not on the main diagonal).
Any help will be greatly appreciated.
Hi, guys!
I've graduated some time ago, and by now have completely fogotten all those matrix rules. Could you please advise me. The problem is following:
I have a large square matrix (like 1000x1000) and I have to rise this matrix to a large power, like 1000. The matrix is quite a sparse one. It looks like this:
x0x00000000000000000
0x000000000000000000
000x00x0x0x0x0x0x0x0x
0xx0xx00000000000000
00x00000000000000000
000000x0000000000000
00x00000000000000000
00000000x00000000000
00x00000000000000000
0000000000x000000000
00x00000000000000000
000000000000x0000000
00x00000000000000000
00000000000000x00000
all x's are different, all zero's are just zero's.
So, my first quess would be that I have to diagonalize the matrix somehow. But I don't remember how to do this. I want to try manually first. Then if it fails - numerically (I prefer fortran).
Any suggestions on how to diagonalized? The right bottom part of the matrix looks like it is already diagonal (although the elements are not on the main diagonal).
Any help will be greatly appreciated.
Last edited: