Single Band Mesh

generate_mesh(bands, W::Function, band_index::Int, n_bands::Int, T::Real, n_levels::Int, n_angles::Int)

Generate a single-band Mesh in a multiband system; the mesh will have (n_angles - 1) x (n_levels - 1) patches per quadrant centered at the Fermi surface using marching squares to find energy contours of bands[band_index] evaluated on a regular grid of $\mathbf{k} \in [0.0, 0.5]\times[0.0, 0.5]$. The width of the Fermi tube is $\pm$ α T.

The function orbital_weights is an array-valued function whose input is a momentum 2-vector, and whose output columns are the the eigenvectors of the single-electron Hamiltonian whose eigenvalues are the functions in bands.

generate_mesh(ε::Function, T::Real, n_levels::Int, n_angles::Int[, N, α])

Generate a Mesh of (n_angles - 1) x (n_levels - 1) patches per quadrant centered at the Fermi surface using marching squares to find energy contours of ε evaluated on a regular grid of $\mathbf{k} \in [0.0, 0.5]\times[0.0, 0.5]$. The width of the Fermi tube is $\pm$ α T.