Orz日記 by Akio Morita

_ [雑記]どたばた

朝10時半ごろから、Υ(5s)への移行作業が始まったはずなのだが...

14:15現在、未だに電磁石の初期化中...なぜだ〜〜〜

日が変わって、00:30現在...未だに大電流運転が安定せず、 六極と格闘を続けている状況 Orz

_ [SAD]SuperKEKB LER Tracking Profile

Google CPU Profilerによる測定結果

Total: 68153 samples
  22354  32.8%  32.8%    22354  32.8% tdrift_
   4828   7.1%  39.9%     4828   7.1% tquad_
   4801   7.0%  46.9%     4801   7.0% tsteer_
   4315   6.3%  53.3%     4315   6.3% ttfrin_
   4228   6.2%  59.5%     4228   6.2% tsolqu_
   4011   5.9%  65.3%     4011   5.9% tmulti_
   2267   3.3%  68.7%     2267   3.3% tbend_
   1745   2.6%  71.2%     1749   2.6% tturn1_
   1353   2.0%  73.2%     1353   2.0% fintrb_
    831   1.2%  74.4%      831   1.2% tmulbs_
    766   1.1%  75.6%      766   1.1% tthin_
    729   1.1%  76.6%      729   1.1% 0x000000020112dee6
    679   1.0%  77.6%      679   1.0% xsin_
    635   0.9%  78.6%      635   0.9% touckf_
    595   0.9%  79.4%      595   0.9% spdrift_
    500   0.7%  80.2%      500   0.7% tsfrin_
    395   0.6%  80.7%      395   0.6% tmultr5_
    303   0.4%  81.2%      303   0.4% spkick_
    302   0.4%  81.6%      302   0.4% tbende_
    287   0.4%  82.1%      287   0.4% bint_
    256   0.4%  82.4%      256   0.4% tintrb_
    238   0.3%  82.8%      238   0.3% 0x000000020112c1fc
    197   0.3%  83.1%      197   0.3% tsteee_
    168   0.2%  83.3%      168   0.2% tcav_
    153   0.2%  83.5%     1307   1.9% _init
    140   0.2%  83.7%      140   0.2% ttfrins_

tdrift(ドリフトスペース)がCPU時間の1/3を占めている

_ [SAD]kblogrdの読み出し速度

afsad1上で 2010/5/31一日分をSAD形式で読み出してみた。(ストア先は、NFSな home directory)

SAD形式(SAD EPOCHな時刻タグ＋レコード名＋データ列)だと格納効率はかなり 低く、5月31日分だけで22GiBの領域を必要とする

_ [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ベースで定義可能と思われる