[T.D.]

Dear Foamers,

I am looking for a "reduce" function but on an "equal operation".

For example, we know that reduce function is already Ok for:
sumOp, minOp and maxOp

How to do it (or extend it) for an "equal operation" and a "pow operation".

I need to make the following code to work also in parallel:

            forAll(mesh.C(),celli)
            {
            Sens[celli] = Foam::pow(SSS[celli], 0.5) ;
            }

           // Sens and SSS are both volScalarField type

Thank you and regards,
T.D.