Algorithms involving Gaussian processes or determinantal point processes typically require computing the determinant of a kernel matrix. Frequently, the latter is computed from the Cholesky decomposition, an algorithm of cubic complexity in the size of the matrix. We show that, under mild assumptions, it is possible to estimate the determinant from only a sub-matrix, with probabilistic guarantee on the relative error. We present an augmentation of the Cholesky decomposition that stops under certain conditions before processing the whole matrix. Experiments demonstrate that this can save a considerable amount of time while having an overhead of less than $5\%$ when not stopping early. More generally, we present a probabilistic stopping strategy for the approximation of a sum of known length where addends are revealed sequentially. We do not assume independence between addends, only that they are bounded from below and decrease in conditional expectation.
翻译:涉及高斯进程或决定点过程的算法通常要求计算内核矩阵的决定因素。后者通常从Cholesky分解法中计算,后者是矩阵大小的立方复杂算法。我们表明,在轻度假设下,只能从子矩阵中估计决定因素,对相对差错有概率保证。我们展示了Cholesky分解法的增强,在特定条件下止于处理整个矩阵之前。实验表明,这可以节省大量时间,而间接费用在不早停时则低于5美元。更一般地说,我们为已知增殖量按顺序显示的总和的近似值提出了一个概率停止战略。我们并不在添加值之间保持独立性,只是它们与下限有关,而有条件期望则减少。