Spheroidal wave functions: which angular normalization?
メンテナーはふだん 1 日以内に返信
まだ誰も着手していません。
評価
調査の方向性
The issue proposes prolate and oblate spheroidal wave functions and asks maintainers to choose an angular normalization. Read the linked discussion in #303 and compare the stated Meixner–Schäfke and Flammer conventions, including the SciPy conversion and reference constraints. Done means a convention is agreed on; implementation and high-precision reference tests are described as follow-up work.
索引モデルが issue の本文から書いたものです。
説明
Following @dschmitz89's request in #303: we're going to add prolate and oblate spheroidal wave functions (characteristic value, angular function of the first kind, radial functions of the first and second kind, and derivatives), ported from Van Buren's MIT-licensed codes and tested against high-precision Arb/MPFR references.
The radial functions are normalized the same way in every convention. The angular function $S_{mn}(c,\eta)$ is not, and we need to pick one:
Meixner–Schäfke (DLMF 30.4, Mathematica's SpheroidalPS, Van Buren's default)
- Includes the Condon–Shortley phase, consistent with
boost::math::legendre_p. - $|S_{mn}| = |P_n^m|$ for every $c$, so the function stays bounded as $c$ grows.
Flammer (scipy.special.pro_ang1/obl_ang1)
- No Condon–Shortley phase; $S_{mn}(c,0) = P_n^m(0)$ for even $n-m$, $S_{mn}'(c,0) = P_n^{m\prime}(0)$ for odd.
- Drop-in for SciPy.
- For oblate functions the norm grows quickly with $c$: $|S_{00}|/|P_0^0|$ is about 1.2e3 at $c=10$, 3.4e11 at $c=30$ and 2.6e24 at $c=60$.
I lean towards Meixner–Schäfke. SciPy can convert with $S^{\text{Flammer}} = \frac{P_n^m(0)}{\mathrm{Ps}(0)},\mathrm{Ps}$ (and use derivatives at 0 for odd $n-m$), with the phase factor $(-1)^m$. The other option is to ship both. Opinions welcome, especially from the SciPy side.
Note that SciPy's current values can't serve as test references: the reference pro_rad2 satisfies the Wronskian $R^{(1)}R^{(2)\prime} - R^{(1)\prime}R^{(2)} = 1/(c(\xi^2-1))$ to 1e-47, while SciPy's differs from it by 1.1% at $(m,n,c,\xi) = (1,2,2,1.5)$ and by a factor of about 7 at $(3,7,10,1.2)$.
- 主要言語
- C++
- スター
- 378
- フォーク
- 271
- 平均マージ
- 12時間 37分
- マージ済み PR（30日）
- 66
環境構築
このプロジェクトには開発コンテナ、Dockerfile、コントリビューションガイドがありません。まず README を読み、一般的な手順ははじめてのコントリビューションガイドを参照してください。
はじめの一歩
- issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
- 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
- リポジトリをフォークし、ブランチを切って変更します。
- issue 番号を参照したプルリクエストを送ります。
boostorg/math のほかの issue
-
難易度 5/5 1週間以上 初心者へのやさしさ 12/100
メンテナーはふだん 1 日以内に返信
-
Reverse-mode autodiff: order-2 tape slows order-1 gradients, rewound tapes keep allocating blocks, rvar copies silently share a node対応中かも @NAThompson が 2 日前に担当しました。 オープン
難易度 4/5 3〜5日 初心者へのやさしさ 20/100
メンテナーはふだん 1 日以内に返信
-
Feature Wishlistオープンenhancement
難易度 5/5 1週間以上 初心者へのやさしさ 20/100
boostorg/math#303 · コメント 20 件 ·
メンテナーはふだん 1 日以内に返信
似ている issue
-
難易度 2/5 1〜3時間 初心者へのやさしさ 78/100
NVIDIA-RTX/NRI#312 ·
メンテナーはふだん 1 日以内に返信
-
permute_q_b_proj / permute_wv_proj accept scale factors from another device対応中かも @0z5a が今日担当しました。 オープン
難易度 2/5 1〜3時間 初心者へのやさしさ 78/100
deepseek-ai/FlashMLA#247 ·
メンテナーはふだん 1 日以内に返信
-
Linux npm packaging leaves workspace files that block subsequent CI checkout対応中かも @Spockkk0225 が今日担当しました。 オープンci
難易度 2/5 1〜3時間 初心者へのやさしさ 73/100
alibaba/neug#1195 · コメント 1 件 ·
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 82/100
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 78/100
メンテナーはふだん 3 日以内に返信