The semantics for default parameters are as follows. If a value function has a
parameter that takes a default value, then an exception is raised if the user
ever tries to assign a value to this parameter via 'params'. These semantics
are chosen to eliminate the possibility that "forgotten" default parameters are
not silently overwritten.

The semantics for default parameters are as follows. If a value
function has a parameter that takes a default value, then an
exception is raised if the user ever tries to assign a value
to this parameter via 'params'. These semantics are chosen
to eliminate the possibility that "forgotten" default parameters
are not silently overwritten.

We now use a polynomial with finite support for our smoothing
function: f(r) = (1 - r**2)**2 Θ(1 - r**2).

This method has been prototyped by Adrien Sorgniard with contributions
by Michal Nowak.  The prototype has been cleaned-up, debugged, optimized,
and integrated into Kwant by Christoph Groth.

It is now impossible to reduce the number of random vectors or moments.
It is also no longer possible to set the number of energy points used by
the method -- this can be supplied by providing a sequence of energies
when calling the KPM object. We replace the 'increase_energy_resolution'
and 'increase_accuracy' methods by 'add_moments' and 'add_vectors'.
The tests have also been updated to reflect these changes.

num_rand_vecs → num_vectors
num_sampling_points → num_energies
ham → hamiltonian
epsilon → eps

e → expr
discrete_coordinates → discrete_coords
substitutions → subs

The optimized routine now replicates the graph faithfully, exactly
preserving the hopping "directions" and values.

Filling now works directly on the graph dictionary, thus avoiding many
unnecessary intermediate steps.  Since filling is essentially graph
replication, the time gains are considerable.

The optimized fill is also more faithful to the template graph.

The new fill() sets each hoppings only once.

The filling starts now at (one or multiple) sites, not at a domain.
This is more consistent, since the flood fill works on the graph of
sites.

Other changes:

* max_sites is now the maximal number of sites that is still OK
* fill() no longer raises an exception when no sites were added, it
  warns now.

Rafal Skolasinski authored 7 years ago and Joseph Weston committed 7 years ago

Previously, only onsite matrix elements were tested, whereas the
test should have also tested the hopping elements.

Additionally, make sure to test the LLL code accordingly.

Rafal Skolasinski authored 7 years ago and Bas Nijholt committed 7 years ago

Calling "kwant.continuum.discretize('A * k_x', substitutions={'A': 'A + B'})"
ends with performing substitutions twice. Fixed by removing redundant second
substitution. Fixes issue #121.

Previously 'Builder.fill' did not respect the symmetry of the target
Builder. This change also adds a test to check that target system
symmetry is respected.

Fixes #119

Rafal Skolasinski authored 7 years ago and Joseph Weston committed 7 years ago

Current implementation uses ``kwant.physics.Bands``. This example fails
when one substitute ``k_x**2`` with for example ``k_x**4`` due to farther
then nearest neighbours hoppings. Using wraparound provides also easier way
to extend it to multidimensional band structures (except of plotting).