Fix bugs. Now (when not reoptimizing the GS), the NI
limit is beautifully reproduced. Also agreement with GKLC
for interacting case.
Something weird (numerics?) going on when reoptimize GS.