ITensor/LAPACK SVD

+1 vote
asked Jun 30 by brendenroberts (220 points)

It appears that the ITensor svd operation does not ultimately call a LAPACK SVD routine. What advantage does this implementation have?

1 Answer

+1 vote
answered Jul 18 by miles (18,230 points)
selected Jul 19 by brendenroberts
 
Best answer

Hi Brendan,
Thanks for the question. This is the case mainly for historical reasons, because we were having some issues with the LAPACK SVD and needed a more reliable SVD as well as a more accurate one.

But it's been on my short list to put in a backend that calls LAPACK for cases where speed is the primary concern. Based on how the tests go, I may make LAPACK the default and allow the current implementation to be turned on with a flag for cases where high accuracy is needed (the current implementation uses a recursive algorithm that corrects errors made in the smallest singular values in earlier passes; it is described in the Appendix of the real-space parallel DMRG paper).

Miles

Welcome to ITensor Support Q&A, where you can ask questions and receive answers from other members of the community.

Formatting Tips:
  • To format code, indent by four spaces
  • To format inline LaTeX, surround it by @@ on both sides
  • To format LaTeX on its own line, surround it by $$ above and below
  • For LaTeX, it may be necessary to backslash-escape underscore characters to obtain proper formatting. So for example writing \sum\_i to represent a sum over i.
If you cannot register due to firewall issues (e.g. you cannot see the capcha box) please email Miles Stoudenmire to ask for an account.

To report ITensor bugs, please use the issue tracker.
...