Orz日記 by Akio Morita

_ [SAD]BENDの実装

BENDのパラメータが同じで、入射CODがほぼ同じケースで、 転送行列が妙にずれているとの話があり、内部実装を調査してみる。 (範囲的には、tturne.f, tbende.f, initb1.f辺り)

まず、内部での有効長　alは、 DISRAD == 0 && FRINGE != 0 && L > 0の場合には以下の補正を受ける $al = L - \frac{1}{24}\frac{(ROTATE * F1)^2}{L}\frac{\sin{\frac{1}{2}(ROTATE - ROTATE * (E1 + E2) - (AE1 + AE2))}}{\sin{\frac{1}{2}ROTATE}}$