[liuhuafei]

Files:
fvMatrix.C
fvMatrixSolve.C
lduMatrixATmul.C

Some question about the coupled boundary.
In fvMatrix.C
In subroutine Foam::lduMatrix::solverPerformance Foam::fvMatrix<Type)::solve。

Firstly we call addBoundarySource(source)，then the 'implicit' part of coupled boundary has been added to the source.

void Foam::fvMatrix<type>::addBoundarySource
(
....
if (!ptf.coupled()){
addToInternalField(lduAddr().patchAddr (patchI), pbc, source);
}
else if (couples)
{
...
forAll(addr, facei)
{
source[addr[facei]] += cmptMultiply(pbc[facei], pnf[facei]);
}

But in the following code, the updateMatrixInterface will add the 'implicit' part to the source once again. If the updateMatrixInterface has the same function, the addBoundarySource doesn't need including the coupled boundary part. Why need add the 'implicit' part twice?

Another question,what the function of the initMatrixinterfaces

In addition,In lduMatrix's AMul and TMul also call the initMatrixInterfaces and UpdateMatrixInterfaces，it seems that we needn't call the initMatrixInterfaces and Updatematrix in the solve subroutine.