Hello everyone, i would appreciate any suggestion you could give me to solve this

Im trying to plot the conductance for a two dimensional array of 3 basis atoms, i have defined the atoms positions, the onsite energies and hoppings matrices and the leads but it is showing me this error:
(i search in the mailing list but theres no info about it)

```ValueError                                Traceback (most recent call last)
<ipython-input-1-ee484ce5ea50> in <module>()
111
112 if __name__ == '__main__':
--> 113     main()

<ipython-input-1-ee484ce5ea50> in main()
107                site_color=family_colors)
108     syst = syst.finalized()
--> 109     plot_conductance(syst, energies=[0.01 * i - 0.3 for i in range(100)])
110
111

<ipython-input-1-ee484ce5ea50> in plot_conductance(syst, energies)
85     data = []
86     for energy in energies:
---> 87         smatrix = kwant.smatrix(syst, energy)
88         data.append(smatrix.transmission(1, 0))
89     pyplot.figure()

371                                                   check_hermiticity, False,
--> 372                                                   params=params)
373
374         kept_vars = np.concatenate([coords for i, coords in

~/anaconda3/lib/python3.6/site-packages/kwant/solvers/common.py in _make_linear_sys(self, sys, in_leads, energy, args, check_hermiticity, realspace, params)
162         lhs, norb = syst.hamiltonian_submatrix(args, sparse=True,
163                                                return_norb=True,
--> 164                                                params=params)[:2]
165         lhs = getattr(lhs, 'to' + self.lhsformat)()
166         lhs = lhs - energy * sp.identity(lhs.shape[0], format=self.lhsformat)

kwant/_system.pyx in kwant._system.hamiltonian_submatrix()

kwant/_system.pyx in kwant._system.make_sparse_full()

ValueError: Hopping from site 407 to site 814 does not match the dimensions of onsite Hamiltonians of these sites.```

Here is the code so far:

import kwant
import tinyarray
import numpy
from matplotlib import pyplot
%matplotlib notebook

VecPrim = [(3.176, 0, 0), (1.588, 2.7504, 0), (0, 0, 17.49)]
base =  [(0, 0, 0), (1.588, 0.9168, 0.8745), (0, 0, 1.749)]

lat = kwant.lattice.general(prim_vecs=VecPrim, basis=base)
c1, c2, c3 = lat.sublattices

c1_ons = numpy.array([[-1.888842, -0.014064-0.000409j, 0.026908+0.003422j],
[-0.014064+0.000409j, -1.784936, -0.031384+0.021565j],
[0.026908-0.003422j, -0.031384-0.021565j, 0.710443]])

c2_ons = numpy.array([[7.168134, -0.005189-0.047376j, 0.001352+0.000103j, -0.001189-0.00095j, -0.048166-0.007663j, 0.429882+0.00052j],
[-0.005189+0.047376j, 0.245358, 0.113162-0.072937j, 0.660849-0.06348j, -0.017571+0.035587j, 0.020961-0.017679j],
[0.001352-0.000103j, 0.113162+0.072937j, 1.338059, -0.003126-0.016285j, -0.624702-0.065652j, 0.002353+0.001639j],
[-0.001189+0.00095j, 0.660849+0.06348j, -0.003126+0.016285j, 1.325337, 0.110694+0.01548j, 0.009041+0.004483j],
[-0.048166+0.007663j, -0.017571-0.035587j, -0.624702+0.065652j, 0.110694-0.01548j, 0.362102, -0.030934+0.002401j],
[0.429882-0.00052j, 0.020961+0.017679j, 0.002353-0.001639j, 0.009041-0.004483j, -0.030934-0.002401j, -0.650839]])

c3_ons = numpy.array([[-3.398186, 0.001839+0.001631j, -0.117366+0.00498j],
[0.001839-0.001631j, -3.401815, 0.023339+0.009067j],
[-0.117366-0.00498j, 0.023339-0.009067j, -0.637207]])

c1_c3_hop = numpy.array([[-1.852931-0.147524j, -0.274338+0.015064j, -0.103127-0.001438j],
[0.274586+0.034451j, -1.796449+0.055331j, 0.004816+0.00237j],
[-0.015051+0.001677j, -0.034396-0.019407j, 3.410475-0.001248j]])

c1_c2_hop = numpy.array([[0.003644-0.00088j, -0.046131-0.04229j, 0.23188-0.001924j],
[-0.305706+0.050514j, -1.963829+0.280977j, -0.597505+0.078772j],
[0.638425+0.083597j, -1.171301+0.088146j, -0.146976+0.012943j],
[-1.105645-0.086829j, -1.112788+0.056971j, -0.210584-0.002248j],
[1.463462+0.012994j, 0.174769-0.106872j, 0.177022-0.031485j],
[-0.01338+0.00264j, -0.058159+0.007987j, 0.953363+0.001855j]])

c2_c3_hop = numpy.array([[0.184152+0.002778j, 0.321743-0.010571j, 0.183538-0.00067j],
[0.275145+0.007691j, -0.54144+0.012815j, 0.045847+0.015379j],
[0.47922+0.015775j, 1.256012+0.007326j, -0.523174+0.00585j],
[1.255561-0.002255j, 1.903063+0.010647j, -0.876392+0.010246j],
[1.282837-0.091022j, -0.54898+0.001161j, 0.057341-0.004095j],
[0.02765-0.001793j, -0.002902-0.001866j, -0.746426-0.000748j]])

def make_cuboid(t=1.0, a=15, b=10, c=5):
def cuboid_shape(pos):
x, y, z = pos
return 0 <= x < a and 0 <= y < b and 0 <= z < c

syst = kwant.Builder()

syst[c1.shape(cuboid_shape, (0, 0, 0))] = c1_ons
syst[c2.shape(cuboid_shape, (0, 1, 0))] = c2_ons
syst[c3.shape(cuboid_shape, (0, 0, 1))] = c3_ons

syst[kwant.builder.HoppingKind((0, 0, 0), c2, c1)] = c1_c2_hop
syst[kwant.builder.HoppingKind((-1, 0, 0), c2, c1)] = c1_c2_hop
syst[kwant.builder.HoppingKind((0, -1, 0), c2, c1)] = c1_c2_hop
syst[kwant.builder.HoppingKind((0, 0, 0), c3, c2)] = c2_c3_hop
syst[kwant.builder.HoppingKind((1, 0, 0), c3, c2)] = c2_c3_hop
syst[kwant.builder.HoppingKind((0, 1, 0), c3, c2)] = c2_c3_hop

lead[c1.shape(cuboid_shape, (0, 0, 0))] = c1_ons
lead[c2.shape(cuboid_shape, (0, 1, 0))] = c2_ons
lead[c3.shape(cuboid_shape, (0, 0, 1))] = c3_ons

lead[kwant.builder.HoppingKind((0, 0, 0), c2, c1)] = c1_c2_hop
lead[kwant.builder.HoppingKind((-1, 0, 0), c2, c1)] = c1_c2_hop
lead[kwant.builder.HoppingKind((0, -1, 0), c2, c1)] = c1_c2_hop
lead[kwant.builder.HoppingKind((0, 0, 0), c3, c2)] = c2_c3_hop
lead[kwant.builder.HoppingKind((1, 0, 0), c3, c2)] = c2_c3_hop
lead[kwant.builder.HoppingKind((0, 1, 0), c3, c2)] = c2_c3_hop

return syst

def plot_conductance(syst, energies):
data = []
for energy in energies:
smatrix = kwant.smatrix(syst, energy)
data.append(smatrix.transmission(1, 0))
pyplot.figure()
pyplot.plot(energies, data)
pyplot.xlabel("energy [t]")
pyplot.ylabel("conductance [e^2/h]")
pyplot.show()

def main():
syst = make_cuboid()
kwant.plot(syst)
syst = make_cuboid(a=100, b=28, c=4)
def family_colors(site):
if site.family == c1:
return 'yellow'
elif site.family == c2:
return 'gray'
else:
return 'blue'
kwant.plot(syst, site_size=0.05, site_lw=0.01, hop_lw=0.01,
site_color=family_colors)
syst = syst.finalized()
plot_conductance(syst, energies=[0.01 * i - 0.3 for i in range(100)])

if __name__ == '__main__':
main()