Advertisement
Guest User

Untitled

a guest
Jun 16th, 2019
87
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.57 KB | None | 0 0
  1. n = 1000000;
  2. A = SparseArray[{
  3. {1, 1} -> 1. n,
  4. {-1, -1} -> 1. n,
  5. Band[{2, 1}] -> -1. n,
  6. Band[{1, 1}] -> 2. n,
  7. Band[{1, 2}] -> -1. n
  8. },
  9. {n + 1, n + 1}, 0.
  10. ]; // AbsoluteTiming // First
  11. M = SparseArray[{
  12. {1, 1} -> 1./(3 n),
  13. {-1, -1} -> 1./(3 n),
  14. Band[{2, 1}] -> 1./(6 n),
  15. Band[{1, 1}] -> 2./(3 n),
  16. Band[{1, 2}] -> 1./(6 n)
  17. },
  18. {n + 1, n + 1}, 0.
  19. ]; // AbsoluteTiming // First
  20. {[Lambda], U} = Eigensystem[{A, M}, -10, Method -> {"Arnoldi", "Shift" -> 1.}]; //
  21. AbsoluteTiming // First
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement