Changed a incorrect J to a Jz in hams/Heisenberg.h
Z3 Model constructor now assures that QN::Nmax() is set properly.
Added IQTensor::diag() method.