This function uses a Greedy Randomized Adaptive Search Procedure
(GRASP) to obtain a partition of m_bin.
Given a phytosociological table (m_bin, with rows corresponding to taxa
and columns corresponding to relevés) this function searches for a
k-partition (k, defined by the user) aiming at high values of the TDV.
See tdv() for an explanation on the TDV of a phytosociological table.
With thr = 1, the algorithm corresponds to the Greedy algorithm.