Use of generic constructor

[quosta]

I am trying to use this code to achieve different tasks: 1) simulating dwave quantum annealing process. Some questions on the parameters:

• for Gamma and Jchain can I use the values that they provide for A(s) and B(s) ?
• what about invTempJ0 and invTempH0? As far as I understand they compensate for the SQUIDs error, but maybe I don't need this level of detail, how can I get rid of them?
• can I obtain the energies? it would be nice to track the eigenvalue spectrum over an s range to see how the gap changes

2) solve Ising models: s = 1, no transverse field, it should be enough to set Gamma = 0. Am I missing anything? Is it sensible to use this library as such?

[jackraymond]

This code was published in support of this paper (and similar lattice equilibrium studies): https://www.nature.com/articles/s41467-021-20901-5. Although it allows generalization, it is not currently very user friendly - apologies.

1. This code simulates a single model (i.e some fixed values of A(s) and B(s)). If you need variation of A(s) and B(s) as per annealing make a feature request (it should be added in the near term). I think the fact you cannot do annealing with this code probably makes your other questions mute for now? In any case invTempJ0 and invTempH0 encode the problem you want to solve (e.g. in the application this code was optimized for (but not exclusive to), basically a triangular lattice with all couplers of 1, but you may want to do a more complicated problem type).
2. If you use this code with Gamma=0 it works fine, and does do classical simulation (Gibbs Updates of classical Ising model). Some of the updates (those relating to non-trivial path integrals) would do nothing, so that adds some inefficiencies, but the algorithm is otherwise not ill-behaved. Again, it doesn't do annealing, so you would be stuck at a fixed temperature.