Segfault in DMRG (user-defined site)

+1 vote
asked Oct 18, 2016 by varelse (170 points)
edited Oct 20, 2016 by varelse

Hi,

I'm trying to define my own SiteSet, basing on Hubbard (with two orbitals per site, and possibly an additional quantum number related to orbitals). The code compiles, but after running DMRG I get a segfault in first sweep. Depending on my initial state, I get it on different stages of first sweep. When running in debug mode, I get:

dmrg-g: /home/varelse/programs/itensor/itensor/itdata/qutil.h:190: void itensor::loopContractedBlocks(const BlockSparseA&, const IQIndexSet&, const BlockSparseB&, const IQIndexSet&, BlockSparseC&, const IQIndexSet&, Callable&) [with BlockSparseA = itensor::QDense; BlockSparseB = itensor::QDense; BlockSparseC = itensor::QDense; Callable = itensor::doTask(itensor::Contract&, const itensor::QDense&, const itensor::QDense&, itensor::ManageStore&) [with VA = double; VB = double]::<lambda(itensor::DataRange, const Labels&, itensor::DataRange, const Labels&, itensor::DataRange, const Labels&)>; itensor::IQIndexSet = itensor::IndexSetT]: Assertion `cblock' failed.

What could this error possibly mean?

(EDIT: I found the error - I defined my physical state IQIndex so that it had two subspaces of equal quantum numbers. )

1 Answer

0 votes
answered Oct 27, 2016 by miles (19,400 points)

Hi - glad to see you're trying some of the advanced features and that you found the error in your code. Let me know if you run into any other issues where I can help.

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.
...