Add function `bonds` to extract bonds of sublattice

[LinjianMa]

Currently, there's a function bonds that returns the bonds of each column/row of a square. We can generalize the version to extract bonds on a general hypercube lattice, with the sublattice being a general SubHyperCubic type.

Reference: #30