DCL:MATH2:SHTLIB : 球面調和関数:サブルーチンの説明
スペクトルデータに対してラプラシアンを演算する.
球面調和関数展開
(4.13)
に対して、水平Laplacian
(4.14)
を作用させると, 球面調和関数の性質から,
(4.15)
となる. そこで,
(4.16)
を導入すると,
(4.17)
と表せる. また, 逆に
(4.18)
であるとき,
(4.19)
を導入すると,
と表せる.
(4.20)
本サブルーチンは, IND=1の場合はAmnから を,
IND=-1の場合はAmnから を計算するものである。
SHTLAP(MM,IND,A,B)
MM (I) 入力. 切断波数(M). IND (I) 入力. ラプラシアンの演算形式を指定する(上記定 義を参照). A (R) 入力. Amnが格納されている配列(長さおよび並び 方は SHTNMLを参照). B (R) 出力. Bmnが格納る配列(長さおよび並び 方は SHTNMLを参照).