Create LICENSE

speed up: use cpp pixel estimator

update readme

very significant speed up: use cpp root finder

added polynomial root finder

significant speed up: 'continue' ASAP

slightly speed up: use kx,ky rather than Point

significant speed up: no need to compute psi

speed up: kill sign multiplication

update readme

update readme

update readme

clean-ups

added test

organize modules