在约束条件下,该方阵简化为Jacobian矩阵,每一行的元素沿着方阵的对角线,等于fn',所有其他元素都等于零。
因此,有必要为每个向量行fn构造对角。为此,我从较大的A(m )矩阵中提取了一个定义为A(hA X 1)的目标向量。然后,我准备了一个定义为C(hA X hA)的零矩阵,它将用来保存对角线。可能有更有效的方法来使用一些预先构建的例程来实现这一
为了使Gauss和Jacobi方法收敛,有必要检查系数矩阵是否为对角占优的,即对角元素在其列的所有元素中应具有最大的值。如果还不是对角线占优势,则采用旋转。如果矩阵是对角占优的,则应满足以下条件:(这也称为收敛)abs(A[i][i]) > summation(abs(A[i][j]),j=1 to n) where j !i...n
//swapping