added polynomial root finder
very significant speed up: use cpp root finder

14s -> 2s