用例とそのコンピューター上での実行に重点を置く テイラー写像の演算と 加速器物理への応用 全6回の講義は 用例とそのコンピューター上での実行に重点を置く Etienne Forest 総合研究大学院大学 2004年12月/2005年1月
“微分代数(DA)”ソフトウェアの原型とは? その働きは? なぜ有用なのか? “演算子の多重定義”とは? テイラー多態とは? “リー/微分代数”とは? 正規形とは? 積分ベースのプログラムとは? 行列ベースのプログラム、Cosy-Infinityとは? コンパイルしたものと解釈したプログラム(Cosy-Infinity)に おける多態の違いは? 以上のもの全てを自分で作成したプログラムで使うことが できるのか? 打切りべき級数代数(TPSA) ソフトウェア テイラー写像の近似値を計算 テイラーを使ったプログラムが 普通のFortranプログラムのようになる テイラーか実数かの選択を実行時まで延ばすことができる 摂動論ではベクトル場を利用する 写像の様々な正規化形式は繰り返しにおいて有用である これらの質問については6回の講義の中で理論よりもむしろできるだけ多くの例を挙げて考えていく。
自動微分(AD)ー微分代数(DA)について 第1回講義 自動微分(AD)は微分代数(DA)ではなく 打切りべき級数代数(TPSA)に基づいている 実用的なベクトル場に関連する微分代数(DA)に ついてはのちの講義で述べることとする
第1回講義における基本的事項: 打切りべき級数代数(TPSA) (すなわち微分代数 (DA)) ソフトウェアの原型はどのようなものか? 第1回講義における基本的事項: 第1回講義における基本的事項: 打切りべき級数代数(TPSA) (すなわち微分代数 (DA)) ソフトウェアの原型はどのようなものか? その働きは? なぜ有用なのか? “演算子の多重定義”とは? これらの問題は Berz が好んだ抽象的な数学的専門用語を使わなければ、難しくはない
無駄を省くために 簡単: この講義のために 十分機能的な小型TPSA ソフトウェア を作成! 空白を含めたった400行のプログラム! 演算子が完全に多重定義され、自動的 2変数 3次数に限定 もちろんこの小型ソフトは Berz のソフトウェアではないので、Berz のソフトウェアで演算子の多重定義をする際の様々な事項を考える必要はない。 ここでの目的はこれらのテクニックをわかりやすく説明することであり、最終的にはほとんど関係のない専門用語でみなさんを煩わせることではない。
小型TPSAソフトウェアを通してTPSA を見てみよう 用い、コンピューター上で計算を展開するための道具である。 光学や加速器においては摂動演算の写像の計算、また時には シミュレーションを可能にする。 演算子の多重定義は自動微分、すなわちAutomatic Differentiationに近いものになりつつあるが、最近は Computational Differentiation という単語が使われている: テイラー展開を計算するプログラムは、もとの数量を計算した プログラムとほぼ同じになる。 摂動演算 : 転送行列、 ラティス関数(トゥイス)、ベータトロン振動数、 減衰、振幅依存チューンシフト、光学システムにおける収差、 共鳴、 シンクロトロン積分など。これらにおいては微分演算子とリー/微分代数(DA)が関係してくる。
小型 TPSAソフトと Gun/FODOセル アプリケーションのメインプログラムの ダウンロードは こちら から! なぜテイラー級数が必要なのか? 小型 TPSAソフトと Gun/FODOセル アプリケーションのメインプログラムの ダウンロードは こちら から! なぜテイラー級数が必要なのか?
一次元四極 FODOセル 一次元四極 FODOセル
TPSAの加法/減法
乗法
逆数 : 冪零の導入
関数: 指数関数 ここでもまた冪零DTを作るために定数部分をまとめる
冪零を無視できる? 時にはOK これはCOSY-Infinity テクニックの中で使われており、 DT は多項式のベクトル場である。 従って、なぜこの級数が使えるのかよく考えてほしい! 小型TPSAの中では TPSA_EXP として使われており、 N=n_tpsa_exp である。
結論 コンピュータ上に打切り級数代数(a Truncated Series Algebra)を導入した。 Berzの “微分代数(DA)” ソフトウェアの原型とはこれだけのものである。 数学的には、この打切りべき級数代数(TPS Algebra)は単なる環である。 微分代数の概念はTPSAに立入るものではない。現段階では何ら関連性はない。
今後の講義 摂動論について幾何学的観点から論じていく : 不変条件、 トゥイス関数およびその非線形歪み 不変条件、 トゥイス関数およびその非線形歪み 幾何学的な図はテイラー写像から導かれた演算に 翻訳される。 自分で作成したFPPライブラリを使用していく。ここでは数学的な詳細には立ち入らないが、このライブラリは公表どおりの働きをすると考えてほしい。 FPPはテイラー写像に関連したベクトル場を扱うことがわかると思う。これらのベクトル場は微分リー代数を 形成する。しかしそれはさほど重要なことではない : 幾何学的側面に注目してほしい。 第2回講義は従来型の計算から始まり、幾何学的観点の導入で終わる予定である。
問題 関数のうち ArcSin(x) と asin(x) は簡単には定数部分と冪零部分とに分けることができない。 x が多項式で表されるとき、 asin(x) の値を求めるアルゴリズムを作成せよ。 それをもとに、 小型TPSA ルーチンを用いてプログラムを作成せよ。 アルゴリズムは3次多項式に限定されてはならない。