[rol97]

Hi everyone,

I am writing a routine in OpenFOAM that requires me to perform the following operations:

- calculate a scalarField (volScalarField works fine as well, I do not need boundaries thought), let us call it myField.

- perform the following operation on the field

for (label i = 1; i < myField.size(); i++)
{
    myField[i] += myField[i-1]; 
}

- then from myField I have to find the closest element to a scalar myScalar, i was thinking of the following loop

forAllReverse (myField, i)
{
    if (myScalar >= myField[i])
    {
    myElement = i;
    break; 
    }
}

I have already implemented and works perfectly in single core, however, when moving to multicore this process obviously does not work, as the same process is repeated separately by each core, thus I am not finding an "absolute" value of myElement, but several indeces.

Any suggestion on how to make it multicore-proof?

Thanks!

Edit: I forgot to say that I am using OpenFOAM v9, I am not sure it that makes any difference!