[elden]

Hi,
I'm trying to run a few tests before converting the following coded FO in an actual OpenFoam FO.

This FO initialises a zero volScalarField in all the grid, then it changes its values ONLY in a previously generated cellZone (by means of topoSet). Each cell value of the volScalarField is now 0 or 1. The output consists of cells number where 0 becomes 1, printed in the logfile.

I noticed differences in the output running the case with mpirun. Looking at the results in paraview it seems that there is an error in forAll cycle: the coded FO changes the values of the new volScalarField outside the defined cellzone.

I thought the problem was related to the definition of the scalar "sum" and that the answer I was looking for was associated with the use of the reduce() command, but I didn't solve it.

[olesen]

An indexing error in your code.
You want something like this instead:

forAll(cz, i)
{
    const label celli = cz[i];
    ....
}

Unless you need the loop index, it is often safer and clearer to use a range-for:

for (const label celli : cz)
{
    ....
}

/mark




@Bruno - relocate to programming?

[elden]

Thanks a lot.
And sorry for the silly question.