def min_div_est(self, LU, nframes): Lu = np.zeros((self.tv_dim, self.tv_dim)) Lu[self.itril] = LU.sum(0) / nframes Lu += np.tril(Lu, -1).T self.Tm = cholesky(Lu).dot(self.Tm)