CFD Online Discussion Forums - full cell stencil with parallelisation

Hi,

I'm trying to get the cell centered values surrounding an arbitrary given cell in a certain volField, i.e. volVectorField U. First I have used a search algorithm that finds for a given cell all neighbour cells and returns the cell values but all in all, it's not working in parallel due to missing diagonal dependencies in the mesh.

Afterwards, I got a hint that the extendedStencil stuff could possibly be used for this. So I'm trying to make it work but ...hmm ... no success till now (especially not in parallel).

Code:

 CPCCellToCellStencil wideStencil(mesh());



 forAll(mesh().C(), lcellI)

 {

       const labelList& cell27Stencil=wideStencil[lcellI];

 

       forAll(cell27Stencil, cellI)

       {

          // now it becomes tricky ...



         // that's not working since the global addressing of the volField is wrong

         Pout << " center value: " << unFilteredField().[cell27Stencil[cellI]] << endl;

        

       }

  }

Furthermore I have no idea how to get correctly the values from the boundary fields (and processor patches) using the stencil stuff.

So I would be grateful for suggestions how to solve this messy problem! Maybe somebody has already done this or has some experience with the stencil stuff.

Best regards

Matthias