这里的基本策略是利用两个矩阵（分别为下三角和上三角）的已知结构。这意味着可以通过前向/后向替换有效地找到反转。例如，对的第一行的运算给出，它（对于您被告知考虑不需要旋转的情况）可以通过简单的标量除法来解决。鉴于您对的了解，您可以为 。在这个过程的最后，你有 N-1 个未知数，一个 N-1 维的右手边和一个$L$$v$$L_{00}x_0=f_0$$x_0$$\tilde{f}_i=f_i-L_{i0}x_0$$i\in [1,N-1]$$N-1\times N-1$下三角矩阵工作。这个递归形成了求解的外循环。关键是存储和以便替换过程访问连续的静态内存。$L$$U$

如果您还没有，我强烈建议您阅读 Press 等人的数字食谱的相关小节。最新版本甚至是 C++。我还建议添加作业标签，就像您在上一个问题中所做的那样。