LduMatrix;;negSumDiag

Вычитает недиагональные записи из диагонали.

for (register label face=0; face<l.size(); face++)
{
    Diag[l[face]] -= Lower[face];
    Diag[u[face]] -= Upper[face];
}

Где:

• Diag это диагональный коэффициент;
• l это lduAddressing для нижнего треугольника;
• u это lduAddressing для верхнего треугольника;
• Lower это нижний треугольник; и
• Upper это верхний треугольник.

С точки зрения discretization coefficients (коэффициентов дискретизации), это эквивалентно:

$ A_p = \sum\limits_r A_r $

а в терминах matrix coefficients (матричных коэффициентов), это:

$ A_{ii} = -\sum\limits_o A_{o,i} -\sum\limits_n A_{n,i} $

Где:

• A это matrix coefficients; и
• o и n это owner and neighbour indices (индексы владельца и соседа).

Ссылка на оригинал - http://openfoamwiki.net/index.php/OpenFOAM_guide/LduMatrix;;negSumDiag