Whitening function now returns the filters used
Made helper function getBMXprod static instead of sub function
tweaked parameters of initial sorting

convergence criterion and drift rate were too high
Cleanup of BPSorter.fit()
Implemented clockwise model estimation of full dataset