gtfold/gtfold — pull request #129 by againerdewar (2014-12-31 21:04:40+00:00)
Title: Fix bogus h bounds in WMPrime loop

In the computation of the WMPrime auxillary array, the index variable h is allowed to range from i+TURN+1 to j-TURN-2. However, referencing Andronescu's thesis, it appears that this range should go all the way up to j-TURN. I have test cases where the current code gives the wrong energy which are fixed by this revision.