Orz日記 by Akio Morita

_ [SAD]SAD->MADX変換

日曜日にinjection studyの様子を見に行ったコントロールルームで Rama達とデバックした結果をまとめておく

• QC1/2のスライスされたMULT(Thick element)をMADXのmultipole(Thin element)へ展開する時のサブスライスの枚数が不足していた[Fixed]
• QCS(MULT in SOL)内部での central orbitが一致しない。結果として、SOL区間の出口での座標変換後の軌道に角度が付いてしまう
• LER wiggler sectionの$\beta_y$が再現しない。ANGLE無しBENDをkickerエレメントに展開する際に linear fringeが落ちてるのが原因かと思ったが、実際はMADXのkickerエレメントが thickなdipole magnetic fieldで実装されていない(エッジから来る vertical focusingが未実装)のが原因であった。[Fixed but ugly hack]

sbendのマニュアルにはk0パラメータは無視される的な記述があるので design orbitを曲げずにビームを蹴るエレメントがkickerか multipoleぐらいしか無いのだが、thickな記述はkickerでしか出来な い Orz

まともに変換出来るようにするたびに、エレメントの分割と matrixエレメントの挿入が増えている...

MADX変数経由のパラメータ埋め込みにすれば、一応マッチング等は可能だと 思うが一部に超越関数が含まれているので、この部分は MADXの実行エンジン にネイティブに実装されて無いとかなり無理っぽい(指数関数とか三角関数系 なので原理的にはマクローリン展開とかで書けるけど...)

こうした変換を行った latticeにどこまで価値があるかどうかは 正直疑問である。

_ [SAD]Long List Support

seq_utils.hを使ったJoinの再実装を作成

インタープリタスタックと関係なく、メモリがある限り長いListを組み立てられるようになった。

実用上、次に必要なのはMap1/Scan1の再実装

多分、コードを読むとイテレーション動作の実装が微妙に違うのだが、定義的にはイテレーションのカーネルが違うだけなので、for_allイテレータとカーネルの形で実装するのが素直か？

実行効率優先なら、現行実装と同様にカーネルをイテレーションループ内にinline展開した方がいいのだが…

テンプレートライブラリ使いたい…Orz

_ [雑記]GPS/1PPSトレンド測定

GPSレシーバ1号機+Intel N305とGPSレシーバ3号機+Ryzen 5950Xのとある1日のGPS/1PPSによるNTP PLLのトレンドを測定した。

• 各レシーバのアンテナは窓に並べて固定
• PC筐体は同一の部屋に設置
• チャートはmuninによる5分毎のサンプリングによる

Ryzen5950X環境だと、散発的に発生する大きなPLL offsetスパイクに引きずられて、PLL freqにスパイクで見える

PLL offsetに見えるノイズ成分は、Ryzen5950X環境側はノイズフロアは低いが、スパイク状のノイズが乗り安定性を損なっているように見える (N305環境は暗雑音多めだが、全般的な安定性は優秀)

GPSレシーバ1号機+Intel N305

GPSレシーバ3号機+Ryzen 5950X