Further inputs to the algorithm are the marginal sample distribution formula_70 which has already been determined by the dominant eigenvector of formula_43 and the matrix valued Kullback–Leibler divergence function