[sguillot]

Dear Foamers


I would like to know how to loop over faces of a cell in paralele.
I currently have the following code :

    forAll(psi.internalField(), celli)         
    {
        const labelList& faces = mesh.cells()[celli];
        
        forAll(faces, fi)                   // Local numbering
        {
            const label facei = faces[fi];  // Global numbering
            
            if (mesh.isInternalFace(facei))
            {
                dgPsi[celli] = max( mag(  d[facei] & gPsi[celli] ), dgPsi[celli]);
            }
            else
            {

            }
        }
 
    volScalarField::Boundary psiBf = psi.boundaryField();
    
    forAll(psiBf,patchi)
    {
        scalarField& psip = psiBf[patchi];
        fvPatchScalarField& pDX = dXbf[patchi];
        fvPatchScalarField& pEps = epsbf[patchi];
        
        forAll(psip,i)
        {
            pDX[i] = Foam::pow(V[patchi], 1.0/3.0);
            pEps[i] = (3.0/2.0) * pDX[i];            
        }    
    }

I don't have the same results for serial and paralele when I want to plot dX for example. I think I have to do something in the "else{ }" (test if patch.coupled() for example) but everything I use is quite bad.



Thank you for your help !