discretizer substititions argument not working as expected
Say I have the Hamiltonian c * k_x
and I want to replace c
by c * c_tunnel
. I would think that I should do the following:
kwant.continuum.discretize('c * k_x', substitutions={'c': 'c * c_tunnel'}, verbose=True)
which results in
def hopping_1(site1, site2, *, c, c_tunnel):
return (-1j*c*c_tunnel**2/2)
this is different from:
kwant.continuum.discretize('c * c_tunnel * k_x', verbose=True)
which results in
def hopping_1(site1, site2, *, c, c_tunnel):
return (-1j*c*c_tunnel/2)
@r-j-skolasinski might know how to fix it.