5 May
2016
5 May
'16
9:48 p.m.
Dear all,
I want to model a random variation of hopping integrals. Simply writing:
def hopping(sitei, sitej):
return -1 + random.gauss(0,0.2)
does not work since this does not ensure that the hopping is
Hermitian. Is there a smart way to write something like:
def hopping(sitei, sitej):
if sitei