Timeline for Efficient way to do ( or Matrix free) Arnoldi
Current License: CC BY-SA 4.0
11 events
| when toggle format | what | by | license | comment | |
|---|---|---|---|---|---|
| Mar 7 at 12:37 | comment | added | Henrik Schumacher | Ah, great. Good that it works for you! Citing Tensors would be misleading as an Arnoldi eigensolver is not really part of it, yet. Just cite this post if you used this code. (On the bottom left of the post click on "Share" and then "Copy link" in the freshly opened window. | |
| Mar 7 at 10:54 | comment | added | Erosannin | I would like to ask for permission if that is okay for you. I have done some modifications and detailed testing of your code, and now I can get a few (hundreds for system size of 2^{18}, the trick is to choose proper Krylov dimension) eigenvalues within machine precision (I used some basic restarting techniques). By benchmarking it, I am going to use it for further calculations for my paper. Is it okay if I cite your "Tensors" package or if there is some other way to cite this? | |
| Mar 4 at 2:25 | comment | added | Henrik Schumacher | Uh yeah, I renamed that function a couple of times. I corrected that and changed ArnoldiBasis to Arnoldi. And yes, eArnoldi[[-1]] should have been some random vector (or some initial guess for an eigenvector). I corrected that as well. | |
| Mar 4 at 2:22 | history | edited | Henrik Schumacher | CC BY-SA 4.0 | deleted 10 characters in body |
| Mar 4 at 0:54 | comment | added | Erosannin | I just started checking it. It seems the Arnoldi function should instead be ArnoldiBasis? Also eArnoldi[[-1]] seems to be undefined as well, which I believe can be some random normalized vector?) | |
| Mar 3 at 22:30 | comment | added | Henrik Schumacher | Glad to be of help. The lack of a matrix-free Arnoldi API has come up sevaral times on this site. And I was curious myself how difficult this would be. Furtunately I already had some code for the Arnoldi function; I had written it as a prototype for the GMRES solver in my package Tensors. | |
| Mar 3 at 22:21 | comment | added | Erosannin | Thank you so much for your time and great answer. It exactly is what I needed. | |
| Mar 3 at 22:13 | vote | accept | Erosannin | ||
| Mar 3 at 20:21 | history | edited | Henrik Schumacher | CC BY-SA 4.0 | added 730 characters in body |
| Mar 3 at 20:09 | history | edited | Henrik Schumacher | CC BY-SA 4.0 | added 684 characters in body |
| Mar 3 at 20:01 | history | answered | Henrik Schumacher | CC BY-SA 4.0 |