Rewrite optimization functions to work with numba

[mariuszslabicki]

Currently, the optimization function is very slow. Probably it is possible to improve performance by rewrite it and use numba.

More info:
http://numba.pydata.org/