*gicosy [#n768f4b1]
光学計算のコード。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の輸送行列 [#ma3803aa]
例えばドリフトの輸送行列は以下のように吐かれている。
-----------------------------------------------------------------------------
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。
* ??? [#t641caa0]
A S E 0.112 0.12;-------------F13
がわからない。コメントにしても結果同じっぽい。GIOSの名残か?