RF: continuing torturing the beast while also trying to figure out its functioning
ENH: first step of big RF -- moved two methods into 2 separate functions

produces the same (divergant from Matlab for lsall?) results are before
added an option to extract only a single EV -- dramastically changes the results bringing them closer to matlab's
final fixes to reach 100% match with matlab implementation for lsall

previous commit on extracing only some good evs did not add other good evs into
nuisances.  also centering of designs should happen after filtering (although
not critical at all) and we should include matching starting time points
(were excluded before)
fixed my own bug for marking derivatives + critical for lsone was to add  EVs own derivative into "others" to match matlab performance

plus added same extract_evs handling into lsone
plus fixed up selection of timepoints to include starting points
plus no collapsing_other should have been done
no krons please for demeaning -- rely on Python powerful broadcasting