核数增加64倍，效率却降到16.3% — 以细化网格方式测量weak scaling的陷阱
·6 min read
在细化网格的弱扩展测试中Krylov迭代增长6.1倍的原因，以及多重网格多加一层如何把效率曲线重新拉回平直。
声速CFL为10时预处理器每次sweep把误差放大100倍 — 块预处理器丢掉的那一支箭头
·7 min read
块Gauss-Seidel只用下三角，丢掉的压力↔速度反馈增益恰好等于声速CFL的平方。本文用一维声学系统实测Schur补在那个位置做了什么。
[论文评述] 手推雅可比矩阵推到心累？——算子重载自动微分（ADOO）
·5 min read
Fraysse（2019）的ADOO不改通量代码，只换数据类型，便取出机器精度的雅可比，有限差分的步长h从此不必再选。
[论文评述] GMRES 跑了 1000 次还停不下来 — 全速度 melt pool 的物理型块预条件
·5 min read
全速度熔池在M=1e-3时，GMRES转了1000次，残差停在1e-2。速度–压力Schur补把迭代次数压回到了两三步。