[vkoppejan]

I trying to follow the topic below:

transforming a function defined on mesh().points() into a pointScalarField

I'd like to do something similar.

1 loop over the the points of a cell and add values to the value at the points
2 Sum up the values off the points on the located on the proc boundaries. In other words, if a point resides on two procs, addup the values on each domain and scatter the result back to the procs again.

How can I perform action 2 using a collective operation on the point field?

Thanks in advance,

Victor

[vkoppejan]

To (partly) answer my own question, this is done usin the syncPointData, as is shown below for (taken from AveragingMethods/Dual/Dual.C ).

Note how "plusEqOp<Type>()" allows summing up of the values.

template<class Type>
void Foam::AveragingMethods::Dual<Type>::syncDualData()
{
    this->mesh_.globalData().syncPointData
    (
        dataDual_,
        plusEqOp<Type>(),
        mapDistribute::transform()
    );
}

[vkoppejan]

To (partly) answer my own question, this is done usin the syncPointData, as is shown below for (taken from AveragingMethods/Dual/Dual.C ).

Note how "plusEqOp<Type>()" allows summing up of the values.

template<class Type>
void Foam::AveragingMethods::Dual<Type>::syncDualData()
{
    this->mesh_.globalData().syncPointData
    (
        dataDual_,
        plusEqOp<Type>(),
        mapDistribute::transform()
    );
}