[論文レビュー] Optimal Quadrature Subsampling for Least Squares Polynomial Approximations

Pranay Seshadri, Akil Narayan+1 more|arXiv (Cornell University)|Jan 20, 2016

Probabilistic and Robust Engineering Design参考文献 26被引用数 4

ひとこと要約

本稿では、QR列ピボットと重み付き直交設計行列を用いて、テンソルグリッドから疎に点を選択する最適な四則サブサンプリング手法を導入する。この手法により、多項式 chaotic 近似におけるwell-conditionedな最小二乗推定が可能になる。ランダム化四則と同等の計算コストで、正規のテンソルグリッドやランダムサンプリングよりもより安定的かつ正確な係数推定を達成する。

ABSTRACT

This paper proposes a new sampling strategy for constructing polynomial chaos approximations for expensive physics simulation models. The proposed approach, optimal quadrature subsampling involves sparsely subsampling an existing tensor grid such that well-conditioned least squares estimates can be computed. By implementing QR column pivoting on our weighted orthogonal design matrix along with a simple total order basis selection strategy, we demonstrate the efficacy of our approach compared with regular tensor grid quadratures and randomized subsampling. We include numerical experiments on an analytical problem and an aerodynamic computational fluid dynamics model. Essentially for the same computational cost as randomized quadrature, the proposed method is able to offer more stable and accurate polynomial coefficient estimates.

研究の動機と目的

• 高価な物理シミュレーションにおける多項式 chaotic 近似のための完全なテンソルグリッド四則の高い計算コストに対処すること。
• 多項式 chaotic 展開における最小二乗係数推定の条件数と正確さを向上させること。
• 計算コストがランダム化サブサンプリングと同等である一方で、より優れた安定性と正確さを達成するサンプリング戦略を開発すること。
• 既存のテンソルグリッドからの疎でwell-conditionedなサンプリングを用いて、信頼性の高い多項式 chaotic 近似を可能にすること。

提案手法

• 既存のテンソルグリッドから疎にサブサンプリングを選択することで、計算コストを低減しつつも正確さを保持する。
• 数値的安定性と条件数の向上を図るために、重み付き直交設計行列にQR列ピボットを適用する。
• 全順序ベース選択戦略を用いて多項式基底を定義し、階層的かつ体系的な展開を保証する。
• 新しいグリッドを生成するのではなく、既存のテンソルグリッドを活用することで、オーバーヘッドを低減し、効率を向上させる。
• 設計行列に重みを付けることで、直交性を維持し、最小二乗問題の条件数を改善する。

実験結果

リサーチクエスチョン

• RQ1 テンソルグリッドからのサブサンプリング戦略は、完全なテンソルグリッドよりもより安定的かつ正確な多項式 chaotic 近似を実現できるか？
• RQ2 同じ計算コストで、最適な四則サブサンプリングはランダム化サブサンプリングよりも条件数と正確さが優れているか？
• RQ3 重み付き直交設計行列を用いたQR列ピボットは、多項式 chaotic 展開における最小二乗推定の安定性をどのように向上させるか？
• RQ4 疎なサンプリング戦略は、必要なシミュレーション評価回数を大幅に削減しながらも正確さを維持できるか？

主な発見

• 提案手法は、完全なテンソルグリッドおよびランダム化サブサンプリングと比較して、より良い条件数の最小二乗問題を達成する。
• ランダム化四則と同等の計算コストで、多項式 chaotic 近似における係数推定がより正確かつ安定的である。
• 解析的問題およびCFDモデルにおける数値実験により、ベースライン手法よりも安定性と正確さが向上していることが確認された。
• QR列ピボットは設計行列の条件数を効果的に向上させ、係数推定における数値誤差を低減した。