あと、Emittance[]の返すTransferMatricesとClosedOrbitだが、どうやらSOL領域は一個ずれているようだ。原因は、src/tsole.fとsrc/tturne.fでの転送行列(iatr)と軌道(iacod)の格納法の問題で、src/tturne.fではl番目のエレメントの入り口の状態をリストのl番目の要素に保存したあとにエレメントの転送行列を作用させているのに対して、src/tsole.fではl番目のエレメントの転送行列を作用させた結果をリストのl番目の要素に保存しているので、SOLで囲まれた境界内部では、TransferMatricesや ClosedOrbitのインデックスがズレている。これと同時に計算されるtwiss配列の更新では、 src/tsole.fではリストのl+1番目の要素へ格納するので整合する。多分、格納コードをコピー＆ペーストしたが保存と転送の順番を入れ替えに対応してコードを書き換えるのを忘れたと推測される。

修正は、Revision 2348で修正をコミットした。

アーキテクチャ的には、SOL領域内部の写像はソレノイド磁場が0となる極限ではSOL領域外と一致するので、単一のコードで処理すべきだと思う。

_ [SAD]SAD->MADX変換

昨日の問題だが、原因が判明

Emittance[]から得た転送行列をmatrixエレメントへ変換する際に、ROTATEパラメータによる不要な座標回転変換を加えていた
SOL内部のインデックスのズレにより、Emittance[]から得た転送行列と軌道との関係が、正しくないのでアフィン変換の定数項が正しく計算出来ていなかった
SOL内部のインデックスのズレに対する、Quick Hackの際にSOL境界を越えてNext Elementを探索したので、FBLX(何もしないMULT)を正しく単位行列に変換出来なかった

一つはこちら側のミスだが、残りはSADの実装バグに起因する Orz

Emittance[]で取り出した転送行列でエレメントを置換するルーチンは完成

SOL系のエレメントを転送行列近似に切り替えて様子を見ると

HERでは、分散関数と軌道の近似が劣化
LERでは、分散関数と軌道の近似が改善(主に垂直方向)

という、悩ましい結果に...

QCSL側で、まっすぐ通過している(LER) or 傾いている(HER)の違いが現れているのだろうか？

別の可能性としては、Emittance[]の返す転送行列が先頭から各エレメントまでのものなので、エレメント毎の行列を行列積で計算することによる誤差混入とかSAD-MADX間の座標系定義の差異とかが疑わしい。

座標系定義に関してはドキュメントをあたって定義を突き合わせてみる。

全部のエレメントを転送行列で展開しても軌道に残差が残る Orz

転送行列とCODをアフィン変換へ展開して、MADX内部で再度転送行列とCODを計算する過程の演算精度が原因か？

次のステップは、SAD上でアフィン変換展開して後でCODの再構築での再現精度を検証することかな？

[ツッコミを入れる]

2013-09-01

_ [雑記]気がつけば9月

ほんとにスケジュール通りにプロジェクトは進捗しているのだろうか orz

_ [雑記]Crab Waist

純粋なものは動作することが示せたが、実際に適用できるかで問題山積み

設置場所があるのか？
最終収束系の non-linear
- superconducting quadrupoleのnon-linearフリンジ
- 各種補正コイルのnon-linear磁場
Sextupole本体のnon-linearフリンジ
Sextupole本体の厚み

[ツッコミを入れる]

2019-09-01

_ [FreeBSD]SMTP-AUTH対応

プロバイダのOCNのOB25なメール中継サービスが、10月より順次SMTP-AUTH/暗号化必須となるとのことなので、メール送信のSMTP-AUTH対応作業を始めた

qmailのSMTP-AUTH/TLS拡張も検討したが、流石に保守があれな状況だし、そこまで設定すると全然簡単ではなくなるので、postfixに移行することに...

で、一応デフォルトのローカル配信系をpostfix化(ローカル配信系は同時立ち上げ可能、/etc/mail/mailer.confで/usr/bin/側のコマンドバックエンドを選択)実施

SMTP-AUTH over SMTPSでの送信テストを始めたが、TLSネゴシエーション後に RCPT TOコマンドへの応答で Sender address rejected: Access denied (in reply to RCPT TO command)で配信拒否される模様

どうやら、Envelope-FromがOCNのアドレスでないとお気に召さない様なのだが、何度か同一の設定でメールを送り出すと受け取って配信してくれるケースがあるのが謎

SMTPサーバーのv4アドレスは同じに見えるので、reversed proxyの先に微妙に設定の異なるSMTPサーバーでも並んでいるのかもしれない

Access deniedで蹴られると送信側のpostfixはそのままbounceしてしまうので、postfix側で送信時のEnvelope-From書き換えを仕込む

sender_canonical_classes = envelope_sender
sender_canonical_maps = regexp:/etc/postfix/envelope_sender_rewrite

当面、これで様子を見る

サーバー側の残作業は、SMTPDとPOPDの実装をqmailからpostfixへの切り替えることだが、Maildir監視用のmaildirwatchはどうしよう？

ports的にはmail/qmailにしか入っていない上にソースはdjb流なので切り出してコンパイルするのも面倒くさい…

[ツッコミを入れる]

カテゴリー: Admin | Emacs | EPICS | Fortran | FreeBSD | GCC | hgsubversion | IPv6 | KEKB | LHC | Lisp | LLVM | MADX | Ryzen | SAD | samba | tDiary | unix | WWW | YaSAI | お仕事 | イベント | 出張 | 宴会 | 数学 | 艦これ | 買いもの | 追記 | 雑記