pep8
pythonize is_adjacent
bugfix: transitions matrix and all of discretize use reversed directionality

than the standard convention
bugfix: plot si, sj from the start of the iteration
improve discretize logging, print progress ratio