|Re: [eigen] Module for orth. Polynomials|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] Module for orth. Polynomials
- From: Manuel Yguel <manuel.yguel@xxxxxxxxx>
- Date: Mon, 25 Aug 2014 10:19:05 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20120113; h=message-id:date:from:user-agent:mime-version:to:subject:references :in-reply-to:content-type:content-transfer-encoding; bh=WYPsXUQmUSIunTJEgEMzlVQki/psvTi23S0fUMhjgXo=; b=vahkrXvK5HsTb38Qmm2Tmxf6DXsC0tiyLf29aPYvEUZhTuC1d/Ald1yfcRAKhT8jK9 vx+VgvH3I+jc0Pb5ju7f71rZMBjH1YURLzeagKduE8y/49NhK8z55afDjWkmZZW8hL8x K8DckHLrpdrCV42CSPNATfAFr7xCBgo4XoDAc1uox4tAx5UilL2BRWcjqbOwU13Ifhay SvpZrNNVPHLNPsO1DhPusKx5yNWagm9oTpPRa6ffi8OF7+DJxf0Rr7LRyIRwpbgY2/6T WOangntClu2pRP9tcLOtTwBaCBsgyujqcRWEg/qK1mdqGJFm5q4rEXZyi163LhOShHP7 qm4A==
I think it is a nice thing to have, I am interested to take a look at it.
Dr. Ing. Manuel Yguel
Porteur du projet StraTagGem
36, rue de l'Université
Tel: +33 9 73 52 86 75
Mobile: +33 6 59 59 17 30
On 08/25/2014 10:16 AM, Roman Pascal Schärer wrote:
Dear Eigen developers and users,
we developed a small C++ package for the computation of orthogonal polynomials (OPQ++) as an unsupported module inside Eigen for our own needs.
The main goal of this package is to provide methods to compute an orthogonal set of polynomials for a given weight function. This has been realised with the modified Chebyshev algorithm, which computes the recurrence coefficients of the polynomials. This algorithm makes use of the moments of the weight function, which in general need to be computed using some numerical integration method. For now we provide our own implementations of quadrature formulas (e.q. Gaussian quadratures). With the obtained recurrence coefficients the nodes and weights of a Gaussian quadrature formula for the weight function can be easily computed using e.g. the Golub-Welsch algorithm.
Since this package could be easily integrated in the Eigen library, I’d like to ask if there is interest for such a module to be made public. Since it is a rather small module, it could e.g. also be integrated in the “Polynomials Module”.
If such interest exists, we would be open for suggestions regarding the API and the implementation of more features such as multi-precision floating point number support etc.