前言

在并行开发时我们经常会用到Pstream::gather()函数或是全局函数reduce()或者其他,需要输入参数Binary &Op,本篇主要讨论Binary &Op是什么

templatevoid reduce(    T& Value,    const BinaryOp& bop, // 这里要输入什么参数    const int tag,    const label comm,    label& request){    NotImplemented;}

Binary &Op

单从名字上看,猜是一个二进制的操作,类似一组操作返回一个二进制的标记
然后去openfoam官网去找,找不到Binary &Op的任何释义
去网上找,发现了一点端倪
c++标准库中有应用Binary &Op这个输入参数,使用的函数也叫reduce,详情可见该网页

这里顺便说几个学习C++不错的网址,查询起来很方便

https://en.cppreference.com/w/https://learncpp-cn.github.io/https://cplusplus.com/https://www.tutorialspoint.com/index.htmhttps://awesome-cpp.readthedocs.io/en/latest/README.html#https://stackoverflow.org.cn/

标准库对输入参数binary_op的解释是:
将以未指定顺序应用于解引用输入迭代器结果、其他 binary_op 结果及 init 上的二元函数对象 (FunctionObject) 。
我们再看下类似的std::transform,详情可见该网页
可能实现的版本有

templateOutputIt transform(InputIt1 first1, InputIt1 last1, InputIt2 first2,                    OutputIt d_first, BinaryOperation binary_op){    while (first1 != last1) {        *d_first++ = binary_op(*first1++, *first2++);    }    return d_first;}

在这个程序中我们看到输入的参数可以传的参数可以函数指针,也可以是仿函数,而且这个仿函数需要传两个参数,按要求送达至返回值中
但按照C++的风格更喜欢是仿函数了
所以我们经常能看到类似这样的函数使用模式

#include #include #include int main () {   bool foo[] = {true,true,false,false};   bool bar[] = {true,false,true,false};   bool result[4];   std::transform (foo, foo+4, bar, result, std::logical_or());   std::cout << std::boolalpha << "Logical OR example as shown below:\n";   for (int i=0; i<4; i++)      std::cout << foo[i] << " OR " << bar[i] << " = " << result[i] << "\n";   return 0;}

以下是输出结果:

Logical OR example as shown below:true OR true = truetrue OR false = truefalse OR true = truefalse OR false = false

那么我类似的把std::logical_or()放到openfoam的reduce函数里可不可以呢
答案是可以的

当然啊,std::logical_or()放到reduce函数中是没什么意义的,但是我们这就可以清楚,我们需要写怎么样的函数,或者函数指针传送到Binary &Op参数中

对于reduce而言,在ops.H中内置了很丰富的仿函数以供调用

Op(sum, x + y)Op(plus, x + y)Op(minus, x - y)Op(multiply, x * y)Op(divide, x / y)Op(cmptMultiply, cmptMultiply(x, y))Op(cmptPow, cmptPow(x, y))Op(cmptDivide, cmptDivide(x, y))Op(stabilise, stabilise(x, y))Op(max, max(x, y))Op(min, min(x, y))Op(minMagSqr, (magSqr(x)=magSqr(y) ? x : y))Op(minMod, minMod(x, y))Op(and, x && y)Op(or, x || y)Op(not, x != y)Op(eqEq, x == y)Op(less, x < y)Op(lessEq, x  y)Op(greaterEq, x >= y)

如果想写自己的Binary &Op操作,可以在ops.H的基础上去创建,要轻松很多


结语

一起探索openfoam也是相当有趣的一件事,非常欢迎私信讨论
指正的价值要比打赏更重要,下面是个人联系方式,希望能结交到志同道合的朋友