Re: [eigen] Numerical Integration module |

[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]

*To*: eigen@xxxxxxxxxxxxxxxxxxx*Subject*: Re: [eigen] Numerical Integration module*From*: Sreekumar Thaithara Balan <tbs1980@xxxxxxxxx>*Date*: Wed, 16 Jul 2014 14:24:19 +0100*Dkim-signature*: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20120113; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :content-type; bh=QNX44OZxTXUVhsbYv6gi6BawnaE70tFMPL2ASlLNcR8=; b=U82/XEPkDynXiUOJULTvV4JUJJOM17q+nKVq6Bwu4/TdQ8GtquN7co0c+vy/NNyZAW dLEC0VJVwcxUXthUJaQauHl4lRGN1kp0aBXEHO5wa3MQSY6Oc2m+dB0IznxyS3LMgpRw 1gHggtHbDWfL0mh2dtSZNHJ8afl/c8g43iyaiVffvN1Ha4khYJKdRflYbvFyAndPXSi0 dKzRlXxXBl3tq8qupQqU/bFXq4e52FAXI1BxPfYrl88aGTKxhUwbU6EuMAv4X1ve3p+o +GjhvLKbbvYI0VuAN2NUsn/SBgID8Kty3E9xun6hZSG+hx+qLhCRsNx7HNlayVCAH9WP sHIg==

Hi Jeff,

I am very much interested in this module as I use both Eigen and Quadpack++ in my code. I am happy to provide help in coding / documenting this module if I get some guidance.

Best,

Sree

On Mon, Jul 14, 2014 at 12:32 AM, Jeff J <complexzeros@xxxxxxxxx> wrote:

I am planning on making it precision-independent. I hadn't really looked into it though. It is good to know about mpreal. That will be convenient. Thanks for reaching out!

On Sat, Jul 12, 2014 at 10:32 PM, Pavel Holoborodko <pavel@xxxxxxxxxxxxxxx> wrote:Hi Jeff,I am interested in such extension.Do you plan to make it precision-independent, so that it can compute with any desired precision using multiple-precision scalar types (e.g. "mpreal")?(most of the code in Eigen is precision-independent and "mpreal" is included in distribution).The main obstacle for this is that QUADPACK uses pre-computed coefficients in "double" precision only.So that we have to implement additional algorithms to compute Gauss-Kronrod nodes & weights for (any) required precision on the fly.The are several algorithms to generate GK:D. Laurie (1997). Calculation of Gauss-Kronrod Quadrature Rules. Mathematics of Computation, 66(219).G. Monegato (1978).. Some remarks on the construction of extended Gaussian quadrature rules. Mathematics of Computation, 32(141).

Here is one of the (unfinished) attempts to implement Monegato method:I am ready to help with any "mpreal"-related questions (as I am its author).Thank you,Pavel.On Sat, Jul 12, 2014 at 2:38 AM, Jeff J <complexzeros@xxxxxxxxx> wrote:I'm planning on converting some QUADPACK code to the Eigen paradigm for my own purposes. Is there interest for this in the Eigen community? If so, I'm interested in hearing thoughts/requests on implementation and API.

Sreekumar Thaithara Balan

Research Associate in Cosmology

University College London

**Follow-Ups**:**Re: [eigen] Numerical Integration module***From:*Jeff J

**References**:**[eigen] Numerical Integration module***From:*Jeff J

**Re: [eigen] Numerical Integration module***From:*Pavel Holoborodko

**Re: [eigen] Numerical Integration module***From:*Jeff J

**Messages sorted by:**[ date | thread ]- Prev by Date:
**Re: [eigen] AVX integer and half-packet support (was Altivec port status)** - Next by Date:
**Re: [eigen] Numerical Integration module** - Previous by thread:
**Re: [eigen] Numerical Integration module** - Next by thread:
**Re: [eigen] Numerical Integration module**

Mail converted by MHonArc 2.6.19+ | http://listengine.tuxfamily.org/ |