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の半分しか無いぞ... 変換行列の展開ミスってる？それとも、非線形項からの寄与が大きい？

_ [FreeBSD]8.0-BETA3

Subversion - CVS間の同期問題等で遅延していたBETA3が出たらしい。

そろそろ、RELENG_8への移行準備を始めるべきかな...

_ [SAD]revision 3200

idlist common blockを修正するために MAC*.incを整理し始めたら、 作業が膨れ上がってキリ番に到達...

-m64を使った darwin buildは現時点では動きませんが、 関数ポインタを格納するのに integer型の idval(*)ではなく TYPE(C_PTR)型の idptr(*)を導入して根本的な解決を行うことを 目標に作業中

_ [SAD]インデックスマジックナンバー

Elementパラメータ参照コードからマジックナンバーを除去するには、 一般にはkytblを参照する必要があるが、特定のElement実装に 対するkytblの内容はコンパイル時もしくはElement実装の 登録時点で確定しているので、呼び出し毎に間接参照を行なう必要は無い。

Cだと、enum型でインデックス定数を生成する若しくは、 関数スコープの static const変数の初期化式が初回呼び出し時に 排他的に評価されることを使って、サブルーチン内のテーブル参照を ローカルなstatic変数の参照に置換する(インデックス計算コードが コンパイル時に畳み込めない分不利)などの手法を使って、 可搬性と効率を両立させる訳だが、Fortranの場合はどうするのが 良いだろう。

• Fortranでは、シングルトンパターンのThread Safeを保証するのが難しい。
• コンパイル時に定数シンボルを割り当てる場合、enum型が無いのでインデックスの重複管理が人力となる。

_ [FreeBSD]stable/14がbranchした

currentからstable/14の枝が生えたので、14.0向けた作業が本格化する模様

直前のcurrent mailing listだとZFSのdeadlockが報告されているので、予定通りにリリースされるかは微妙な気が…