コアを64倍にしたら効率が16.3%まで落ちた — 格子を細かくして測るweak scalingの罠
·8 min read
格子を細かくするweak scalingでKrylov反復が6.1倍になった理由と、マルチグリッドのレベル1段がその線を平らに戻す場所。
音速CFL 10で前処理器が誤差をsweepごとに100倍にした — ブロック前処理器が捨てた矢印一本
·9 min read
ブロックGauss-Seidelが下三角だけで捨てる圧力↔速度の帰還利得は、ちょうど音速CFLの二乗です。Schur補元が何をするかを1D音響系で測ります。
[論文レビュー] ヤコビアンを手で微分するのに疲れたら — 演算子オーバーロード自動微分（ADOO）
·6 min read
Fraysse（2019）のADOOは流束コードをそのままに型だけ替え、機械精度のヤコビアンを手計算なしで取り出す手法です。
[論文レビュー] GMRESが1000回回っても止まらなかった — 全速度 melt pool のための物理ベース・ブロック前処理
·6 min read
M=1e-3でGMRESが1000回、残差は1e-2で止まった。速度–圧力のSchur補元が、その反復を二、三回に固定する。