[論文レビュー] A Robust Approximation to a Lambert-Type Function

K. Roberts|arXiv (Cornell University)|Apr 8, 2015

Sports Dynamics and Biomechanics参考文献 3被引用数 8

ひとこと要約

この論文では、ハリーエの方法を用い、適応的初期推定値を採用することで、$ y = \log(W(e^x)) $ の関数を、言語に依存しない堅牢なアルゴリズムで計算する手法を提示する。この方法は数値的アンダーフロー・オーバーフローを回避し、すべての表現可能な実数 $ x $ に対して信頼性高く動作し、数回の反復で高い精度を達成する。

ABSTRACT

The function $y = g(x) = \mathrm{log}\big(W(e^x)\big)$, where $W()$ denotes the Lambert W function, is the solution to the equation $y + e^y = x$. It appears in various problem situations, for instance the calculation of current-voltage curves for solar cells. A direct calculation of $g(x)$ may be inaccurate because of arithmetic underflow or overflow. We present a simple algorithm for calculating $g(x)$ that is robust, in that it will work for almost all $x$ values which are representable in the arithmetic of one's chosen computer language. The algorithm does not assume that the chosen computer language implements the Lambert W function.

研究の動機と目的

• 浮動小数点演算におけるアンダーフローまたはオーバーフローのため、$ y = \log(W(e^x)) $ の計算に生じる数値的不安定性を解消すること。
• ネイティブなラムベルト W 関数の実装がない任意のプログラミング言語でも、この関数を信頼性高く評価できる方法を提供すること。
• 極端な値を含む、$ x $ の全実数領域において、高い精度と安定性を確保すること。
• 太陽電池のモデル化などの応用に適した、計算的に効率的かつ数値的に堅牢なアルゴリズムを開発すること。

提案手法

• 関数 $ y = \log(W(e^x)) $ は、$ y + e^y = x $ を解く形に再定式化され、ハリーエの方法を用いて反復的に解かれる。
• 初期推定値は、$ x \leq -e $ の場合 $ y_0 = x $、$ x \geq e $ の場合 $ y_0 = \log(x) $、$ -e < x < e $ の場合線形補間を用いる。
• ハリーエの方法は、導関数 $ h'(y) = 1 + e^y $ および $ h''(y) = e^y $ を用い、一般のハリーエ更新式から導かれた反復式を適用する。
• 各反復で $ e^y $ を一度だけ評価し、収束基準または固定反復回数を用いて精度を保証する。
• 更新式の分母において、最初の項が支配的になるように設計されており、キャンセルエラーを最小限に抑える。
• この方法は $ x $ の値が $ \pm 10^6 $ に達する場合でも堅牢であり、4回の反復で $ 10^{-80} $ 未満の誤差を達成する。

実験結果

リサーチクエスチョン

• RQ1 直接評価がアンダーフローまたはオーバーフローのリスクを伴う場合、$ y = \log(W(e^x)) $ をどのようにして堅牢に計算できるか？
• RQ2 全実数直線にわたり、この関数の反復ソルバーの収束を速く保証する初期推定戦略は何か？
• RQ3 ハリーエの方法は、$ y + e^y = x $ を高精度かつ最小限の計算コストで効果的に解くことができるか？
• RQ4 組み込みのラムベルト W 関数がなくても、任意のプログラミング言語でこのアルゴリズムをどのように実装できるか？

主な発見

• ハリーエの方法を4回反復するだけで、$ x \in [-10^6, 10^6] $ の範囲で絶対誤差が $ 10^{-80} $ 未満に抑えられる。
• 関数 $ g(x) = \log(W(e^x)) $ は厳密に単調増加かつ凹関数であり、大きな負の $ x $ では $ x $ に漸近し、大きな正の $ x $ では $ \log(x) $ に漸近する。
• 関数 $ g(x) $ が収縮的であるため、数値的に安定しており、$ |g'(x)| \leq 1 $ であるため、入力ノイズへの感度が小さい。
• ハリーエの方法を用いることで、立方収束が達成され、粗い初期推定値からも迅速に精度が向上する。
• 更新式の分母におけるキャンセルエラーは、式の中で最初の項が支配的になることで回避される。
• この方法は、$ \log(W(e^x)) $ の高精度評価が不可欠な太陽電池の電流-電圧モデル化などの実世界の問題に適用可能である。