如何获取相邻节点的标量值

[Mr.H]

各位大神：
我在计算的时候需要计算某一节点的相邻节点的标量场平均值。该如何操作？
T.internalField( ) 需要知道网格单元的index。
我能想到的方法是 根据 polyMesh中的 faces owner neighbor 文件重构出网格单元的相邻关系。但感觉这个方法太笨，有没有别的方法？

[李东岳]

const labelList& neighbour = mesh.cellCells()[cellI]

scalar sum(0);

forAll(neighbour, i)
{
    sum += T[neighbour[i]];
}

sum /= neighbour.size();

[Mr.H]

谢谢大神，我看到你原来的代码里有相关的。再次感谢。
http://www.cfd-china.com/topic/1228/owner-neighbour-网格寻址/2