[doubtsincfd]

In the code for gauss gradient (gaussGrad.c), the grad calcuation at the boundary is done as:

forAll(vsf.boundaryField(), patchi)
{
if (!vsf.boundaryField()[patchi].coupled())
{
const vectorField n
(
vsf.mesh().Sf().boundaryField()[patchi]
/ vsf.mesh().magSf().boundaryField()[patchi]
);

gGrad.boundaryField()[patchi] += n *
(
vsf.boundaryField()[patchi].snGrad()
- (n & gGrad.boundaryField()[patchi])
);
}
}

So is there no correction for coupled patches?? I am getting embarrassingly wrong values of gradient at cyclic patch.

Second question: In the same code, calculation of grad has been done for three sets of mesh entities:
1. On face centers on patch
2. On cell centers adjacent to patch
3. On all inside cell centers

Is my understanding correct?