On $hp$-Convergence of PSWFs and A New Well-Conditioned Prolate-Collocation Scheme

技術要約：PSWFのhp収束性と、新しい良条件化されたプロレート・コロケーション・スキームについて

問題提起
プロレート・スフェロイダル波関数（PSWF）は、帯域制限関数の近似にスペクトル基底を提供し、有限および無限区間の両方における直交性や、最適な時間・周波数集中性といった特有の性質を備えている。これらは、擬スペクトル法やスペクトル要素法の優れた代替案として提案されてきたが、その適用には重大な理論的および数値的な障壁が存在する。具体的には、PSWFは非多項式基底であるため、標準的な多項式に基づく数値的特性が成立しない。最近、Boydらによって指摘された重要な問題は、プロレート要素法における$h$細分化（$h$-refinement）の非収束性であり、これは標準的なスペクトル要素法の収束理論と矛盾する挙動である。さらに、PSWFの積に対する厳密な求積則が存在しないため、ガラーキン定式化の魅力が低下しており、標準的なコロケーション・スキームはしばしば悪条件（ill-conditioning）に陥る。

手法
著者らは、以下の3つの主要な手法開発を通じてこれらの課題に対処している。

1. hp収束性の理論的解析： 本論文では、PSWF投影に関する$hp$誤差評価を確立している。分解されたソボレフ空間上の投影作用素を分析することにより、著者らはメッシュサイズ$h$と多項式次数$N$への依存性を分離した厳密な誤差評価を導出している。
2. Kong-Rokhlin則の実装： コロケーション点の数$N$に対する帯域幅パラメータ$c$の選択は、精度において極めて重要である。著者らは、積分および微分における高精度を保証するために$(c, N)$を組み合わせる、Kong-Rokhlin則の実用的な近似を提案している。この近似は、Stirlingの近似から導かれた上界を利用することで、計算コストの高いPSWF固有値$\lambda_N$の計算を回避する。
3. 良条件化のための新しい基底の構築： 標準的なプロレート・コロケーション系の悪条件性を克服するため、著者らは「二重の性質」を持つ新しい基底$\{\beta_k\}$を構築している。
   • この基底は、第2次導関数が基底関数のカーディナル基底特性を近似するように設計された、一般化されたBirkhoff補間問題を通じて定義される。
   • 決定的な点として、内部点の基底関数$\beta_k$は、PSWFそのものではなく、PSWFの二重積分によって生成される空間を用いて構築されている。これにより、悪条件な微分行列の逆行列を計算することなく、安定した計算が可能となる。
   • 結果として得られる行列$B_{in}$は、2次微分行列$D^{(2)}_{in}$の近逆行列として機能し、システムを前処理（プリコンディショニング）する。

主な貢献

• 非収束性の厳密な証明： 本論文は、プロレート要素法における$h$細分化の非収束性に対する理論的な正当化を提供している。導出された誤差評価（定理4.1）は、$N$を固定した場合、誤差項のうち$h$に依存しない項が支配的になるため、$h$を細分化しても$L^2$ノルムにおいて収束しないことを示している。これは、標準的なルジャンドル多項式の近似が示す$h$収束とは対照的である。
• 良条件化されたコロケーション・スキーム： 著者らは、条件数が帯域幅パラメータ$c$および点数$N$の両方に依存しない新しいPSWF基底を導入している。2つのアプローチが提示されている：
  1. 新しい基底行列$B_{in}$を用いて標準的なコロケーション系を前処理する手法。
  2. 新しい基底を直接用いてコロケーション形式を構成する手法（この場合、最高次微分の行列が単位行列となる）。
• 最適化されたパラメータ・ペアリング： $(c, N)$をペアリングするための実用的なアルゴリズムが提案されており、これにより、プロレート微分行列の離散固有値が連続作用素の固有値を高い精度で近似することを保証している。これは、ルジャンドル法と比較して、スペクトルの大幅に広い範囲において、最大12桁の精度を実現する。

結果

• 固有値の解像度： 数値実験によれば、提案された$(c, N)$ペアリングを使用した場合、プロレート微分行列の離散固有値の約87%が、12桁の精度で連続固有値を近似している。対照的に、ルジャンドル法では、この精度を達成できるのは固有値の約25%のみである。
• 条件数とソルバーの性能： 新しい良条件化されたスキーム（P-PCOLおよびN-PCOL）は、$N$が増加しても安定した条件数（1.3～1.9程度）を示す。一方、標準的なプロレート・コロケーション・スキーム（PCOL）は、条件数が$O(N^4)$で増大する。その結果、高次の$N$において、標準的なスキームでは反復回数が数千ステップに及ぶのに対し、新しいスキームを用いた反復ソルバー（BiCGSTAB）は、一定のステップ数（6～7ステップ）で収束する。
• ヘルムホルツ方程式： 不均質媒体における高波数領域のヘルムホルツ方程式の解決において、新しい$p$版プロレート要素法は、ルジャンドル基底を用いたスペクトル要素法を大幅に凌駕する。この新しいアプローチは、より粗い格子で高い精度を達成し、高度に振動する解に対して急速な収束率を示す。

意義
本論文は、PSWFが帯域制限関数に対して独自の利点を持つ一方で、その単純な適用は、非多項式基底であるために標準的なスペクトル要素法の収束理論が直接適用できないという点で欠陥があることを主張している。本研究の意義は以下の通りである：

1. 限界の明確化： プロレート要素における$h$細分化の失敗について厳密な理論的説明を提供し、標準的なスペクトル要素の収束理論が非多項式基底に直接適用できるという仮定を修正した。
2. 実用的な利用の実現： 堅牢で良条件化されたコロケーション・フレームワークを提供することで、PSWFを帯域制限関数や高度に振動する解の問題に対する、有力かつ優れた代替手段とした。
3. 優れた性能： 適切な$(c, N)$ペアリングと新しい基底を組み合わせることで、提案手法が連続スペクトルの解明および不均質媒体における高波数ヘルムホルツ方程式の解決において、ルジャンドル多項式ベースの手法を大幅に上回ることを示した。

著者らは、PSWFの積に対する厳密な求積則が存在しないことから、PSWFを用いたスペクトルアルゴリズムには、ガラーキン法よりも提案された良条件化コロケーション法を用いることが最適であると結論付けている。