Fast and accurate spectral treatment of coordinate singularities (Q808218)

The coordinate singularities along the polar axis can decrease the accuracy of classical spectral methods. For that reason, the authors consider two axisymmetric functions \(U_ k(r,\theta)\) (0\(\leq r\leq 1\); \(0\leq \theta \leq \pi\); \(k=1,2)\) with the Fourier expansions \(U_ k(r,\theta)=\sum^{\infty}_{m=0}U_{k,m}(r)r^ m \cos (m\theta),\) where every so-called Robert coefficient \(U_{k,m}(r)\) is bounded at \(r\to 0\). The goal is to compute the corresponding expansion of \(U_ 3:=U_ 1U_ 2\). By fast Fourier transforms, the Robert coefficients of \(U_ 3\) are calculated in a well-conditioned manner.