4.2 并行范式与代码谱系

4.2 并行范式与开源代码谱系

本节摘要：把算法切到上千个进程上，靠的不是“并行很重要”这句空话，而是领域分解、节点内共享内存、以及负载再划分。个别时间步与均匀切域互相讨厌；AMR 的层次补丁让空洞进程无事可做。开源代码随后把这些选择冻结成方言：Gadget 一脉、RAMSES 与 Enzo 的网格、FLASH 的爆炸传统、AREPO 的运动网格、Athena 的 MHD 测试题文化。学会读谱系，比学会背模块名更能判断一篇论文的隐性假设。

本节地图

阅读完本节，你应当能够：

1. 用领域分解说明粒子或格子如何分家，边界如何交换
2. 解释个别时间步为何让负载均衡变成动态问题
3. 区分 MPI、节点内线程、加速器三种层次，各自管哪段墙
4. 把常见开源码放到 SPH/网格/运动网格/MHD 谱系上
5. 判断“我们改了一处反馈”在某方言里实际碰到的是哪一层基础设施

切开空间，而不是切开公式

宇宙学盒子最自然的切法是切空间：每个进程管一块，粒子走进走出时打包发送。网格码切的是格子，AMR 切的是补丁。消息里走的是边界上的幽灵层或穿越粒子。切得越碎，表面积越大，消息税越高。强扩展因此不是线性的道德，它有几何上限。弱扩展——问题与核数一起长大——才是宇宙学真正在用的尺子。

个别时间步把这幅图画脏。密晕里的粒子步长可以比空洞里短出许多数量级，均匀切域会让管晕的进程成为整场的节拍器。动态再划分把密集区切细，把空洞并给更少进程。再划分本身要通信，不能每步都做。星团直接 N 体更极端：共享列表或复制力计算，空间切分不总是对的。并行范式因此跟着对象走，不存在统一的“天体物理 MPI 模板”。

对象           自然切法                 最疼的负载
  均匀早期宇宙 规则域                   后来的晕集中
  AMR 爆炸     层次补丁                 细区集中在冲击波
  星团         不一定切空间             硬双星的时钟
  辐射蒙特卡罗 粒子包或域               光学厚区的随机游走
​

💡 关键直觉：并行图上的进程不是平等公民。物质去哪，公民权就跟到哪，否则你在为空洞付租金。

三层墙，以及代码如何变成方言

第一层是节点间：MPI 消息。第二层是节点内：线程吃共享缓存。第三层是加速器：把齐步核丢过去。三层都写进同一时间步，调试会变成考古。成功的方言会把三层藏在基础设施里，让写反馈的人少碰幽灵层。Gadget 把 TreePM 与 SPH 做成宇宙学联邦的普通话；其后续把运动网格接进来，方言发生过一次著名的语法改革。RAMSES 用八叉树网格把 AMR 宇宙学写成另一种普通话。Enzo、FLASH 把块结构 AMR 带到星形成与核爆炸。Athena 与 PLUTO 用测试题文化训练 MHD 共同体。ChaNGa 与 PKDGRAV 把树码并行推到极大 N。读论文方法节，看见这些名字，应读成“继承了哪一种切域与哪种流体刀”，而不是“用了著名软件所以可靠”。

软件工程在这里不是附属。版本、输入参数的可重复、测试题库、对单元物理的回归，决定方言能不能活过导师更替。没有测试题的私有码可以发高水平可视化，却很难被审计。开源也不是自动的美德：接口不稳定、文档只存在于口头，一样会让谱系断裂。真正活下来的是那些把入场考试写进仓库的。

改一处物理，会碰到哪一层

学生常以为改星形成阈值只是改一个数。在 Gadget 方言里，这个数挂在粒子密度估计上，密度又挂在核半径与邻居数上。在 RAMSES 方言里，它挂在最细层格子与杰恩斯长度上。同一物理句子，基础设施不同，收敛行为不同。黑洞吸积更明显：有的方言用邻居气体的均值，有的用键接的吸积半径。并行的域边界若把邻居切开，还要额外的消息。所以“我们实现了同样的反馈”常常是假的字面平等。

编程语言与脚本传统也是谱系的一部分。核心热循环长期是编译型语言，分析与参数扫描越来越多用脚本。把分析写进仓库、与演化版本钉在一起，才谈得上可重复。只保存最终可视化，等于只保存铸币的图案，不保存模具。

⚠️ 常见坑：跨方言对比星系质量函数，却不统一星形成阈值的操作定义。你在对比基础设施，不是在对比宇宙。

概念账本：

对象 → 切空间或切列表
聚集则动态再划分
三层：消息、线程、加速器
改物理时写下：密度如何估、邻居如何跨域
对照必须声明方言与测试题版本
​

代码把演化跑完之后，产物可能大到无法全部运回家。下一节是数据洪流：就地分析、光锥和可视化如何成为方法本身。

测试题库、回归与方言的寿命

活得过二十年的方言，几乎都有一组夜间回归：激波管、静水、宇宙学功率谱、磁转子。回归失败则禁止合并。这套纪律比许可证类型更能预测代码是否可审计。没有回归的分支，可视化可以继续，谱系却在暗中断裂——你不知道哪一次“优化”把角动量偷走了。

参数文件是方言的口语。同一物理量在不同码里名称不同、单位不同、默认值不同。社区对照项目必须先做词典。学生把默认值当物理常数，是最常见的入门事故。把关键默认值写进论文表格，比写“我们使用了某码”有用得多。

版本钉死还涉及依赖库：消息实现、数学库、编译器优化等级。宇宙学功率谱对优化等级通常不敏感，近距正则化可能敏感。可重复声明若只给一个版本号，没有编译选项，在星团问题上可能不够。这些听起来像运维，却决定第 1 章那本守恒账能不能在别人机器上再记一次。

强扩展的谎言与弱扩展的诚实，以及文件层并行

在固定问题规模上把核数加到极限，效率必然掉，因为表面积税与串行段会现身。用强扩展曲线证明“我们的码更优越”，而不报问题规模，是常见的展示术。宇宙学真正依赖弱扩展：盒子与核数一起长。弱扩展曲线上的折点，往往对应再划分或全局傅里叶的瓶颈。把折点写进架构图，比写峰值占用率有用。

输入输出是被忘记的并行层。全家福若由单进程串行写出，演化再快也会在快照处排队。并行文件、每个进程写分块再后合并，是数据管道的入口。第 4.3 节的洪流，有一半在这里已经形成。方言若把 IO 当附属，用户会为了少写快照而损害并合树。基础设施决定科学提问，再次成立。

问题：两个项目都写“实现了同样的超新星反馈”，能否直接比质量函数？

先对齐操作定义。SPH 邻居与 AMR 最细层格子的体积不同，同一焦耳的命运不同。密度如何估计、邻居如何跨域、随机种子如何同步，都属于定义。方言不同，字面平等是假的。再看测试题库与默认参数表。默认值被学生当成物理常数，是入门事故。跨码对比必须先做词典，再做宇宙。强扩展曲线若固定小问题，还可能是展示术；宇宙学真正依赖弱扩展。

弱扩展才是宇宙学的尺子：盒子与核数一起长，折点对应再划分或全局变换瓶颈。固定小问题的强扩展曲线容易变成展示术。输入输出是被忘记的并行层，全家福若串行写出，演化会在快照处排队。用户为少写快照而损害并合树，是基础设施在改科学。回归测试、参数词典与编译选项，决定方言能不能在别人机器上把守恒账再记一次。

谱系树怎么长：私有模具变成联邦普通话

1990 年代的宇宙学码多半是实验室私有模具：能跑、能出图，换一组人就断代。Gadget 把树与 SPH 做成可携带的联邦件之后，后续版本在同一棵树上接个别时间步、TreePM、再接运动网格，语法改革发生在继承里，而不是另起炉灶。SWIFT 一脉把任务图和缓存友好的粒子排序写成另一种并行口音。ChaNGa 把树码挂到可迁移对象上，负载再划分变成运行时的默认动作。Athena 一脉把 MHD 测试题文化带到更齐步的核。读这些名字，应看成河道改道：同一条物理河，河床被并行模型挖成不同断面。

空间填充曲线后来成为切域的常用模具：粒子按曲线排序，邻域在内存里更可能相邻，幽灵层更瘦。曲线不能消灭个别时间步的时钟分裂，但能降低“物质去了邻域却要跨节点取数”的税。星团直接码走另一条河：复制力或共享列表，因为空间切分会把硬双星拆到两个进程上，时钟更乱。切法跟着对象走，这句话在谱系树上能被点名验证。

社区对照项目是谱系的体检。对照失败时，先查默认参数词典，再查守恒账，最后才谈物理。把对照差异直接写成宇宙学结论分歧，是把模具公差当成新定律。开源许可证救不了没有回归的分支；夜间激波管、静水和功率谱，才是方言能活过导师更替的保险。

要点串联

• 并行首先切空间：表面积就是消息税
• 个别时间步强迫动态负载：均匀公民权是幻觉
• 三层墙不要混在热循环里裸奔：方言的价值是把墙藏好
• 开源码是方言不是勋章：读切域与流体刀
• 同一反馈句子在不同方言里不是同一操作
• 可重复靠测试题与钉死的分析版本

下一节面对拍字节：快照策略、就地统计与可视化，如何决定模拟结果究竟能不能被看见、被审计。