前言
在并行开发时我们经常会用到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也是相当有趣的一件事,非常欢迎私信讨论
指正的价值要比打赏更重要,下面是个人联系方式,希望能结交到志同道合的朋友