Multidomain spectral method for the Gauss hypergeometric function

Siegfried Crespo 1 Marco Fasondini 2 Christian Klein 1, * Nikola Stoilov 1 Corentin Vallée 1, 3
* Auteur correspondant
3 Empenn
INSERM - Institut National de la Santé et de la Recherche Médicale, Inria Rennes – Bretagne Atlantique , IRISA-D5 - SIGNAUX ET IMAGES NUMÉRIQUES, ROBOTIQUE
Abstract : We present a multidomain spectral approach for Fuchsian ordinary differential equations in the particular case of the hypergeometric equation. Our hybrid approach uses Frobenius’ method and Moebius transformations in the vicinity of each of the singular points of the hypergeometric equation, which leads to a natural decomposition of the real axis into domains. In each domain, solutions to the hypergeometric equation are constructed via the well-conditioned ultraspherical spectral method. The solutions are matched at the domain boundaries to lead to a solution which is analytic on the whole compactified real line R∪∞, except for the singular points and cuts of the Riemann surface on which the solution is defined. The solution is further extended to the whole Riemann sphere by using the same approach for ellipses enclosing the singularities. The hypergeometric equation is solved on the ellipses with the boundary data from the real axis. This solution is continued as a harmonic function to the interior of the disk by solving the Laplace equation in polar coordinates with an optimal complexity Fourier–ultraspherical spectral method. In cases where logarithms appear in the solution, a hybrid approach involving an analytical treatment of the logarithmic terms is applied. We show for several examples that machine precision can be reached for a wide class of parameters, but also discuss almost degenerate cases where this is not possible.
Type de document :
Article dans une revue
Liste complète des métadonnées
Contributeur : Imb - Université de Bourgogne <>
Soumis le : vendredi 26 juillet 2019 - 08:57:11
Dernière modification le : vendredi 13 septembre 2019 - 09:50:02

Lien texte intégral



Siegfried Crespo, Marco Fasondini, Christian Klein, Nikola Stoilov, Corentin Vallée. Multidomain spectral method for the Gauss hypergeometric function. Numerical Algorithms, Springer Verlag, 2019, pp.1-35. ⟨10.1007/s11075-019-00741-7⟩. ⟨hal-02194789⟩



Consultations de la notice