[T.D.]

Hi Foamer,

I have the following code that uses a Foam:ow function. The solution is working fine on one processor, but the solution is not correct (distorted) if i do run the code in parallel, any ideas ?

-----------------------------------------------------------------------------
            Epow = DF;
            forAll(mesh.C(),i)
            {
              Sens[i] = Foam::pow(Epow[i], 0.5); 
            }
            alpha = Foam::max(a,Foam::max(alpha-half,Foam::min(b,Foam::min(alpha+half,alpha*Sens*lambda))));
-----------------------------------------------------------------------------

Note : The types of the different variables above are either a volScalarField or a dimensionedScalar


regards,
T.D.

[hk318i]

I don't think the problem comes from ``pow`` function. It is probably due to using ``max`` and ``min`` which are working locally on each processor. You can try ``gMin`` and ``gMax`` which work globally.

[T.D.]

Hi Hassan,

Thanks. It tried, but the problem seems to be due to looping on mesh.C(). For example, i tried to update a "volScalarField alpha" in the following simple code:

scalar k=0.;
scalar gg=mesh.C().size();
forAll(mesh.C(),i)
{
alpha[i]=k/gg;
k++;
}
alpha.write();
return 0;

It works fine in serial, but not in parallel. In parallel the field "alpha" is distorted at the processors shared interfaces ! ?

Foamers, any ideas, help ?

Thanks
T.D.