better comments for BD diagonalization.
Making the length PCA code less silly.

Lost a number of allocations and conversions between GSL matrices and
float array arrays by moving some intermediates to arrays of GSL
vectors.

Passes tests.
doc tweak