本次参考:Yulan Fang。这纯属我个人的自学笔记,因为palabos的参考资料数时有限。Yulan Fang大佬真的很强,参考他(她)的资料,一步一步来。没有任何商业用途,如果介意,联系我,立删!!
第四个(tutorial_4),我也不清楚为啥没有3。总体上和第二个是一个功能,但是好像增加了速度范数的记录。在这之前,我尝试着改一下总的迭代步数和间隔多少步数保存一次图像(切记更改完一定先make再“./”)。同时我把这个程序中输出范数文件和图像放到一个文件中。
唯一增加了:
“Array<T,2> velocity”定义了一个两个元素的velocity;“lattice.get(nx/2,ny/2).computeVelocity(velocity)”则是计算中心点的速度。palabos中pcout表示输出:输出的 是中心点处的速度。“Box2D line“来定义一条线;"*computeVelocityNorm(*extractSubDomain(lattice, line)) "用来计算线上的范数;”setprecision(3)“表示三位有效数字;
"在Palabos,输出文件使用plb_ofstream而非ofstream,以确保并行运行正常。computeAverageDensity(lattice)与computeAverageEnergy(lattice)是Palabos预先就有的功能,计算平均密度和平均能量。" -------Yulan Fang.
范数显示: