我正在尝试为不完整的矩阵分解实现最小丢弃填充（MDF）排序算法。算法描述在第 60 页用于可压缩 Navier-Stokes 方程的 Newton-Krylov 算法的预处理技术。

有一件事我不明白。对于初始化，我应该使用算法 7（第 62 页）

使用算法 7 计算 A0 图中所有节点 vj 的丢弃值

但是这个算法使用的是稍后计算的值。我们只有具有 k-index 的变量

使用以下迭代计算丢弃值，discard(vk) = ||Fk||F

还有 k+1-index 变量

列弗(k+1)

我不明白它是递归还是关于迭代方法的东西......如果任何人都容易理解这个 MDF 算法，请告诉我如何实现这一刻。