关于同位网格Lid driven cavity SIMPLE算法的一些问题，CFD

大家好，我最近在尝试使用同位网格的SIMPLE算法来解决Lid driven cavity
的问题。我使用的软件是matlab。我参考了一些代码，其中一个很符合我的预期。但是，我不能明白在function solve里面，对于动量的离散的一些内容。根据我看的一些资料，在迭代中，对速度进行了松弛，但我不能理解的是，这里面的u(i,j)为什么没有在后面进行修正，而是一直在迭代，以及u_star是一直在修正的。以下是部分代码。

希望大家可以帮我解答！万分感谢！

这是在Jacobi迭代，最近我也在学同位网格，可以分析一下这个代码吗？