Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- n = 10000;
- M = N[SparseArray[
- Join[Table[{i, i} -> 1.0, {i, 1, n}],
- Table[{i, i - 1} -> 2.0, {i, 2, n}],
- Table[{i, i + 1} -> 2.0, {i, 1, n - 1}]]]];
- val = Eigenvalues[M, -1][[1]]
- A=coo_matrix((data, (row, col)), shape=(nn, nn)).toarray()
- eigs(A, 1, sigma=0)[0][0]
- (-0.00013788630102868301+0j)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement