Orz日記 by Akio Morita

_ [SAD]LINEバグ

LINE["S", position + frac]の動作が、Twiss[*, position + frac]系と異なるバグが昔から知られているが、少し深堀りしてみた

原因は、positionの次のbeamline-elementがOFFSET付きMARKerの場合、beamlineに沿った内部配列の次の要素が物理的に異なる場所を指しているにも関わらず、線形補間するためであるが、自体はもう少しややこしい模様

明らかに線形の変化で無い要素に関しては、qfraccompで一時的にbeamline-elementのパラメータをfracでスケールして部分的に再計算しているため、GEOなどは正しく扱われている模様

現時点で判明している正しくないLINEのkeywordは以下の通り

• S,LENG - itfpos配列の線形補間
• GAMMA,GAMMABETA - itfgamm配列の線形補間
• SIG* - 光学関数は再計算しているが、itfgamm配列の線形補間を結果のスケーリングに使っている

調査過程で、SIZE*keywordに関しては、itfsize配列の割り当て状況を検査していないので、SEGVで落ちるケースが発見された

Rev.6409時点では、S,LENG,GAMMA,GAMMABETAの補間は修正済み、SIZE*に関しては未割り当て時はエラーを返すように変更済み