This clustering method computes the clusters via performing replicator dynamics. It has been shown that the solutions of a replicator dynamics correspond to the solutions of the following quadratic program [1]}, [2]}.
\(\max _\mathbf {v} \; f(\mathbf {v}) = \mathbf {v}^{\texttt {T}} \mathbf {X} \mathbf {v}, \;\; \texttt {s.t.} \; \mathbf {v} \ge \mathbf {0}\, , \sum _{i=1}^{n}\mathbf {v}_i = 1\, ,\) 
