并行效率疑问

[Micro]

@random_ran 在 并行效率疑问 中说：

以后多核计算前，我都会renumberMesh 的

您好！请问，在进行多核计算时，我应该先decompose,再renumber,还是应该 先renumber再decompose呢？谢谢！

[浪迹天大]

@Micro 网上见到的例子都是先分块，再并行地renumberMesh

[Micro]

@浪迹天大
好的，谢谢您！并行renumber应该是 renumberMesh -parallel 吧
能把您见到的例子分享下吗，谢谢！另外，您清楚为什么重排序后能够加速求解吗？感谢！

[浪迹天大]

@Micro 最近正在用：mpirun -np 20 renumberMesh -parallel -overwrite。运行后会告诉你：

After renumbering :
Mesh size: 2687852
Before renumbering :
    band           : 32282
    profile        : 2.13313e+10

    band           : 1914
    profile        : 2.13794e+09

具体原理不清楚。我也没对比过。不过我相信它至少不会使结果变坏，也不会使计算变慢。

[Micro]

@浪迹天大 谢谢！谢谢大佬！

[Micro]

@浪迹天大 最后一个问题大佬！如果是并行renumber，是全局（整个网格）整体重新编号，还是分配到每个进程上的网格子区域独自重新编号呢？感恩！

[Micro]

@浪迹天大 如果是全局（整个网格）整体重新编号，那么为啥非要先区域分解再renumber呢？感觉也可以先renumber再decompose吧？谢谢！

[浪迹天大]

@Micro 肯定是局部啊，不然就不需要并行模式了。

[Micro]

@浪迹天大 @浪迹天大 那运行后输出的mesh size和band(也就是矩阵带宽)是全局矩阵还是局部矩阵呢？感恩！

另外我突然发现您居然是openfoam成长之路的博主，敬仰！

[浪迹天大]

@Micro 不敢。看了一眼代码，mesh是全局的，band是局部的：

Info<< "Mesh size: " << mesh.globalData().nTotalCells() << nl
	<< "Before renumbering :" << nl
	<< "    band           : " << band << nl
	<< "    profile        : " << profile << nl;

[Micro]

@浪迹天大 好的，谢谢！band是局部的，并且info语句只有主进程执行，所以band是分配在主进程上的网格对应的矩阵的带宽，对吗？可是我认为，每个进程都保留了全局A的一部分，难道renumber只是针对这个部分矩阵进行排序吗

[浪迹天大]

@Micro 抱歉，这个工具的原理和实现我都不懂。无法做出更多判断。

[Micro]

@浪迹天大 好的，我自己看下，有进展再向您汇报

[winsway_zero]

@random_ran 您好，想咨询一下，这么多的网格您是如何进行后处理的呢？paraview 会崩溃的吧？