[Shibi]

Hello to all,

In a functionObject I am coding, I need to access the values of pressure of cellID {20, 70, 30, 50, 90}, in this order.


I using a 2D square domain with 10x10 cells for testing decomposed in 2 processors.


I have the following code.

   List<scalar> cellIds{20, 70, 30, 50, 90};
   List<scalar> pValues(cellIds.size() );

   const volScalarField& sField = lookupObject<volScalarField>("p");

    forAll(cellIds, cellI)
    {
          pValues[cellI] = sField [ cellIds [cellI] ];

    }

In serial, this will run without a problem. In parallel the domain is decomposed into 50 cells for each processor. So, when I try to access
sField[70], it will give me the following error:

[1] 
[1] 
[1] --> FOAM FATAL ERROR: (openfoam-2012)
[1] index 70 out of range [0,50]
[1] 
[1]     From void Foam::UList<T>::checkIndex(Foam::label) const [with T = double; Foam::label = int]
[1]     in file /home/pc/OpenFOAM/OpenFOAM-v2012/src/OpenFOAM/lnInclude/UListI.H at line 138.
[1] 
FOAM parallel run aborting
[1] 
[0] 
[0] 
[0] --> FOAM FATAL ERROR: (openfoam-2012)
[0] index 70 out of range [0,50]
[0] 
[0]     From void Foam::UList<T>::checkIndex(Foam::label) const [with T = double; Foam::label = int]
[0]     in file /home/pc/OpenFOAM/OpenFOAM-v2012/src/OpenFOAM/lnInclude/UListI.H at line 138.
[0] 
FOAM parallel run aborting

Is it possible of access the local cellID on each processor based on the global cellID?


Any builtin function to do this?


Best Regards