最近の計算機事情

空間 ずいぶんとご無沙汰しました。5年ぶりの再会ということになります。
先生 それ少し変だよ。君は大学4年生から年を取らないことにしてあるのだから。
空間 そんなこと言ったてこのホームページは1996年から、ぜんぜんさわられていませんよ。
先生 その通りだね。
空間 そんなことはしませんよ。空間群を使えば,何とかなると思って相談にきたのですよ。大体のことはTSPWAMとTSLCLAの中にある。射影演算子を使っているところをまねすればと分かるので自分で作ってみました。平面波はそっくりまねで問題がありませんが,球関数で表わされているところが分かりません。
先生 平面波はそっくりまねと言っても,平面波をしまっているCOMMON/SPW/をそれぞれのkで作らないといけないだろう。球関数の部分はTSLCLAが立方調和関数を使っているのに,君のバンド計算のプログラムは球関数をそのまま使っているから,困ると言うことだね。
空間 そうなんですよ。おそらく少し面倒でもTSTRLMを使って一度立方調和関数にしてTSLCLAのまねをして射影したあとで,もう一度もどすというのが手近な方法ですよね。
先生 それが良さそうだね。もう一つはEULERCがオイラー角をあたえるので,TSRMIで回転群の表現行列を作ってから射影するというのがあるが,こちらは今度の本で説明を省略した回転群の勉強をしてからということになる。
空間 どれどれ,あっP252の表題がまちがっていますよ。このTSIRMIの三番目のIは余分ですよ。
先生 ここに球関数を点群の操作で回転するサブルーチンTRSHHMが作ってあるから持っていってもいいよ。
1996/05/03