The problem is that this array has size 2 and in hsmt_nonsph the maximum
value of this array is used for some array allocation. If the system under
investigation only has a single spin the 2nd value is undefined.

(found with valgrind, only theoretical error source so far)