calc.yの文法に関しては、 ";",")",ID型トークンによる exprの終端を受け入れる(YYACCEPTする)構文に変更しています。これは、calc.yが解析する構文が、部分構文で次のトークンが構文定義に含まれないために、SADの MAINレベル的には正常な構文に対して構文解析器が必ずsyntax errorで終了するという問題への修正です。

_ [SAD]itopenbuf_のデフォルトを変更

Revision 1500から USE_NEW_FORTRAN_OPEN_SIM=YESな状態がデフォルトになりました。

古い実装を使いたい場合は、USE_OLD_FORTRAN_OPEN_SIM=YESを sad.confに追加してください。

_ [SAD]SPAC系フラグに関連するコード

なんか、SAPC系のフラグが乱立しているので、状況を整理してみる

PSPAC
- tfsetparam.f(PSPAC*変数の読み込み)
- tturn.f(tspac_ルーチンの呼び出し @ トラッキング)
WSPAC
- track.f(WSPACフラグの初期化)
- tturn.f(twspac_ルーチンの呼び出し @ トラッキング)
- qtwiss.f(Space Charge付きのOptics計算)
- temit.f(Space Charge付きのEmittance表示)
- tintrb.f(Space Charge付きのEmittance計算)
SPAC
- track.f(SPACフラグの初期化)
- tturn.f(DRIFTエレメントに対してspdrift_を呼び出す)
- tmulti.f(spkick_とspapert_を呼び出す)
- tsol.f(DRIFTエレメントに対してspdrift_を呼び出す)
- tftrack.f(SPACフラグ付きの場合は、NPARAを殺している)

なお、Space Charge関連と思われるtspac.fにある tspac_ルーチンは、呼び出されない。

ざっと眺めた感じでは、次の問題がありそう

PSPAC、WSPAC、SPACは同時にオンに出来る
PSPACとWSPACに関しては、NPARA>1の際にトラッキングが並列化される。
- 並列実行時は、エレメント毎でのプロセス間同期は保証されないので、粒子間相互作用や、粒子分布の統計量を使う計算は「結果が保証できない」はず

[ツッコミを入れる]

2023-03-20

_ [雑記]tDiary & Hiki

運用環境のtDiaryを最新版に、Hikiは 1.0.0ベースでRuby 3.1対応改修版に更新して、Ruby3.1での運用に切り替え実施

Hiki環境に関しては、Ruby2.7 EoL後にtaint回りのmethodを削除して廻る必要あり(Ruby 3.2で削除予定)

[ツッコミを入れる]

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