memo/gicosy

gicosy †

光学計算のコード。gicosyというシェルスクリプトで走らせる。これを改造して以下のようにした。

cp $1 GICOSYIN.DAT
./gicosy.x
./meta2ps META.DAT
sed s/"XX.XX.XXX"/"`date`"/g -i gicosyplot.ps
sed s/"YY.YY.YY"/""/g -i gicosyplot.ps

OUTFILE=$1
OUTFILE=`echo $OUTFILE | sed -e "s/inp/out/g"`
cp GICOSYOUT.DAT $OUTFILE

PSFILE=$1
PSFILE=`echo $PSFILE | sed -e "s/inp/ps/g"`
cp gicosyplot.ps $PSFILE

OUTPUTの輸送行列 †

例えばドリフトの輸送行列は以下のように吐かれている。

-----------------------------------------------------------------------------
                           U1 = 0.000000000E+00 M     V1 = 0.000000000E+00 M
 DRIFT LENGTH = 1.500000000E+00 M                     W  = 0.000000000E+00 DEG
                           U2 = 1.500000000E+00 M     V2 = 0.000000000E+00 M
 -----------------------------------------------------------------------------
Entrance of STQ15

 NON SYMPL. SYSTEM TRANSFER MATRIX AT PATH-LENGTH L= 1.500000000E+00 M
 *********************************************************************
                X             A             Y             B             L
  -----------------------------------------------------------------------------
  0        0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00  1.500000E+00
  -----------------------------------------------------------------------------
  1 X      1.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
  2 A      1.500000E+00  1.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
  3 Y      0.000000E+00  0.000000E+00  1.000000E+00  0.000000E+00  0.000000E+00
  4 B      0.000000E+00  0.000000E+00  1.500000E+00  1.000000E+00  0.000000E+00
  5 G      0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
  6 P      0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00

つまり転置された行列とみることができる。上の例では始点角度の終点位置への寄与は(x|a)=1.5、始点位置の終点角度への寄与は(a|x)=0。

??? †

A S E  0.112 0.12;-------------F13

がわからない。コメントにしても結果同じっぽい。GIOSの名残か？