Orz日記 by Akio Morita

_ [SAD]*frac系の継続調査

• tfbeamfrac - gammab配列の線形補間に依存している
• qtwissfrac1 - gammab配列の線形補間に依存している

tfgeoによるpos/gammab/geo配列の構成は、tfgeo1による計算直後は、物理並びであるが、tfgeoがMARKer要素に関しては、OFFSET参照先の値を計算してpos/geoを更新しているが、gammabは明示的に更新していないので、tfgeoの設計意図的にはgammab配列上のMARKer要素の値は一つ前の物理element出口の値を保持している状態を意図しているものと思われる

しかし、fractionalなgeometryを計算するtfgeo1の読み出しで、再計算の対象となったelement出口側のgammabは作業領域として使われたあと復元されていない また、MARKerからのoffset計算に置いて、再帰的なMARKer列が考慮されていない (qcell1やtturneでは、tffselmoffsetによる再帰検索が行われている)

また、*frac関数系の動作としては、有限のfracional lengthかつqfraccompが偽(長さの更新が光学系を変更しない)場合、次の要素の値を返している

従って、MARKer要素の有限fractinal positionは、次のelementを指すことになる

qfraccompが真になり得るElement Type

• DRIFT
• BEND
• QUAD
• SEXT
• OCT
• DECA
• DODECA
• ST
• MULT
• CAVI
• TCAVI

他のTypeはすべて偽となる

例えば、

• UND
• WIG
• SOL
• MAP
• INS
• COORD
• BEAMBEAM
• PROT
• APERT

本質的に局所的に不連続な写像を含むSOL, COORD, APERTに関しては厚みが無いものとして扱うほかないと思われる

BEAMBEAMに対しては意味づけの議論が必要

MAPに関しては、非物理的な任意写像を含むので、一般に分割を定義するのは困難と思われる

一方、進行方向に広がった磁場分布をモデル化するUND, WIG, INS, PROTには中点を意味づけ可能

例えば、PROTは、逆写像が存在すればMatrixLog/Expベースで定義可能と思われる