import matplotlib.pyplot as plt
from tenpy.models import lattice
fig, axes = plt.subplots(1, 2, sharex=True, sharey=True, figsize=(5, 6))
orders = ['default', 'snake']
lat = lattice.Honeycomb(4, 3, None, bc='periodic')
for order, ax in zip(orders, axes.flatten()):
    lat.order = lat.ordering(order)
    lat.plot_order(ax, linestyle=':')
    lat.plot_sites(ax)
    lat.plot_basis(ax, origin=-0.25*(lat.basis[0] + lat.basis[1]))
    ax.set_title(repr(order))
    ax.set_aspect('equal')
    ax.set_xlim(-1)
    ax.set_ylim(-1)
plt.show()