Replace subset_phenotype_columns with generic alternative

 - new function is subset_columns, for any matrix, in math/matrix.d
Remove trailing whitespace