[hrvig]

Hello,


I'm making a post processing utility where I need to access information across cyclic boundaries. At some point in my code I have something like:



forAll (mesh.cellCells()[x], celli)
{
Info << mesh.cellCells()[x][celli] << endl;
}


This gives me ID's of neighbouring cells to cell x but not across cyclic boundaries. Does anyone know how to access this information?


Best regards,
Jakob

[jherb]

You can find an example of how to do this in here:


https://github.com/OpenFOAM/OpenFOAM...hScalarField.C


The relevant lines are:

     // Get the coupling information from the mappedPatchBase
     const mappedPatchBase& mpp =
         refCast<const mappedPatchBase>(patch().patch());
     const polyMesh& nbrMesh = mpp.sampleMesh();
     const label samplePatchi = mpp.samplePolyPatch().index();
     const fvPatch& nbrPatch =
         refCast<const fvMesh>(nbrMesh).boundary()[samplePatchi];
 ...
      const turbulentTemperatureRadCoupledMixedFvPatchScalarField&
         nbrField = refCast
              <const turbulentTemperatureRadCoupledMixedFvPatchScalarField>
             (                 nbrPatch.lookupPatchField<volScalarField, scalar>(TnbrName_)
             );
 ...
     // Swap to obtain full local values of neighbour internal field
     scalarField TcNbr(nbrField.patchInternalField());
     mpp.distribute(TcNbr);

[hrvig]

Thanks for getting back. However, with my limited experience I cannot apply it to my case.

Anyone else who can provide me with a simple working code? To be more specific I'm trying to access a volScalarField across the cyclic boundary