Orz日記 by Akio Morita

_ [SAD]calc.y修正

先日の日記で書いてた calc.y回りの保守を実施。

新しいコードでは、src/yylex_.cへ

• Frotranから呼び出すyyparse_,yypushtoken_,yypoptoken_
• calc.yから呼び出すトークンリーダーyylex

を実装し、calc.yからyaccにてcalc.hとcalc.cを生成し そのまま使います。

yaccの出力を Cコンパイラに渡すために、calc.yのアクションは FortranからCへ書き換えて有ります。

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の際にトラッキングが並列化される。
  • 並列実行時は、エレメント毎でのプロセス間同期は保証されないので、粒子間相互作用や、粒子分布の統計量を使う計算は「結果が保証できない」はず