Orz日記 by Akio Morita

_ [FreeBSD][数学]libmががが！

Fortran2008で$\Gamma$関数と対数$\Gamma$関数が標準に入ったので、 使おうと思って調べていたのだが、最終的には libmを呼び出す訳で 呼び出し先のコードを追跡していたら... gamma(3m)は$\Gamma$関数ではなく、 対数$\Gamma$関数の値を返すんだね... Orz 何のために、lgamma(3m)が有るんだか src/lib/msum/src/e_gamma_r.cの __ieee754_gamma_rは$\Gamma$関数の対数を返すと言う コメントを読んでる当たりで、頭痛がしてきた。

$\Gamma$関数を返すのは、tgamma(3m)です。

はい、皆さん準備はいいですか？

右見て、左見て、半顔、見上げて...

絶望した！libmの紛らわしい仕様に絶望した！

_ [SAD]SAD->MADX変換

厚み0のBENDエレメントの写像を tbthin()@src/tbendfから拾ってきて変換器に実装。 分散関数の一致が改善したが、依然としてSOL領域の 分散関数がSADモデルと一致しない Orz

SOL境界の座標変換の実装にバグが有るのかなぁ？

HERの変換結果を見る限り、SOLへの境界直後の分散関数の微分が 有っていないというかSADの半分しか無いぞ... 変換行列の展開ミスってる？それとも、非線形項からの寄与が大きい？