[論文レビュー] FaVeST: Fast Vector Spherical Harmonic Transforms

Quôc Thông Lê Gia, Ming Li+1 more|arXiv (Cornell University)|Jul 31, 2019

Geophysics and Gravity Measurements参考文献 68被引用数 4

ひとこと要約

本稿では、球面上の接ベクトル場に適した高速な前向きおよび随伴ベクトル球面調和変換（FwdVSHT および AdjVSHT）を計算するための FaVeST というアルゴリズムを紹介する。この手法は、スカラー球面調和変換とクレブシュ＝ゴルダン係数を活用することで、近似的に線形の計算複雑度を達成する。FwdVSHT の計算コストは 𝒪(N log √N)、AdjVSHT は 𝒪(M log √M) に低減され、シミュレートされた接ベクトル場を用いた検証により、高い精度、効率性、数値的安定性が確認された。

ABSTRACT

Vector spherical harmonics on the unit sphere of $\mathbb{R}^3$ have broad applications in geophysics, quantum mechanics and astrophysics. In the representation of a tangent vector field, one needs to evaluate the expansion and the Fourier coefficients of vector spherical harmonics. In this paper, we develop fast algorithms (FaVeST) for vector spherical harmonic transforms on these evaluations. The forward FaVeST evaluates the Fourier coefficients and has a computational cost proportional to $N\log \sqrt{N}$ for $N$ number of evaluation points. The adjoint FaVeST which evaluates a linear combination of vector spherical harmonics with a degree up to $\sqrt{M}$ for $M$ evaluation points has cost proportional to $M\log\sqrt{M}$. Numerical examples of simulated tangent fields illustrate the accuracy, efficiency and stability of FaVeST.

研究の動機と目的

• 球面上の接ベクトル場を扱う応用分野における、効率的なベクトル球面調和変換の計算アルゴリズムの不足に対処すること。
• FwdVSHT および AdjVSHT の計算コストを、従来の 𝒪(N²) および 𝒪(M²) から近似的に線形の複雑度に低減すること。
• 地球物理学、天体物理学、流体力学への応用を想定した、実用的でソフトウェア実装済みの高速変換ソリューションを開発すること。
• さまざまなノード数を用いたシミュレートされた接ベクトル場における数値実験を通じて、提案手法の精度、安定性、効率性を検証すること。

提案手法

• クレブシュ＝ゴルダン係数を用いて、発散なしおよび渇きなしのベクトル球面調和関数をスカラー球面調和関数の明示的表現で表すこと。
• これらの表現を基に、FwdVSHT および AdjVSHT をスカラー球面調和変換の合成として定式化すること。
• スカラー調和係数の評価を高速化するために、球面上の高速フーリエ変換（FFT）技術を適用すること。
• 変換計算における積分の高精度離散化を保証するために、多項式正確な数値積分則を用いること。
• 実用的利用および既存のソルバへの統合を想定し、MATLAB ソフトウェアパッケージとしてアルゴリズムを実装すること。
• さまざまなノード数を用いた、シミュレートされた接ベクトル場における数値実験を通じて、手法を検証すること。

実験結果

リサーチクエスチョン

• RQ1 既存の高速スカラー調和変換を用いて、接ベクトル場のためのベクトル球面調和変換を近似的に線形時間で計算できるか？
• RQ2 クレブシュ＝ゴルダン係数をどのように活用することで、計算の高速化が図れるスカラー調和関数によるベクトル球面調和関数の表現が可能になるか？
• RQ3 評価点の数が異なる条件下で、提案された高速アルゴリズムの達成可能な精度と安定性はどの程度か？
• RQ4 このアルゴリズムはソフトウェアとして効率的に実装可能であり、球面上の偏微分方程式ソルバに統合可能か？

主な発見

• 前向き FaVeST は、N 個の評価点に対して 𝒪(N log √N) の計算コストで実行され、従来の 𝒪(N²) の複雑度と比べて顕著に低減されている。
• 随伴 FaVeST は、M 個の評価点に対して 𝒪(M log √M) の時間で動作し、調和係数からベクトル場を迅速に再構成可能である。
• 数値実験の結果、前向き変換における相対誤差は N と共にゆっくりと増加し、正規化誤差が 𝒪(1/N) であることが示され、良好な収束性と安定性が裏付けられた。
• 前向きおよび随伴変換において、N の増加に伴い正規化誤差の包絡線が減少するという事実から、高い数値的安定性が確認された。
• ソフトウェアパッケージは正常に動作しており、すでに球面上の接ベクトル場のテンソル・ニードレット変換の高速化に応用されている。
• 今後は、Python での HEALPix グリッド対応を含め、球面上の Stokes やナビエ＝ストークス方程式を解くソルバへの統合を想定した今後の研究が行われる。