[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Help-gsl] looking for gsl_integration_long_double

From: Tobias Munk
Subject: Re: [Help-gsl] looking for gsl_integration_long_double
Date: Thu, 31 May 2007 13:42:08 +0200 (CEST)

The problem I have to solve is a summation of integrals of Mathieu function 
(some special
periodic functions, in gsl since v 1.9). The integral solver for oscillatory 
does not work in this case, since the integrand is known only numerically, thus 
cannot factor out a trigonometric term. In any case, a long double integration 
would not be sufficient to increase accuracy for this problem, since the Mathieu
functions must be available at the same precision. Currently, they are only 
present as
double, too.

Mathematica is able to approximate the series much better than the code I wrote 
gsl, presumably since Mathematica uses arbitrary precision arithmetics. I am 
not sure
which approach is more promising, implementing everything in long double, or 
arbitrary precision arithmetics, like gmp (


--- Evgeny Kurbatov <address@hidden> schrieb:

> Please tell, what the task you solve with such requirement to precision?
> Tobias Munk wrote:
> > I am looking for a integration routine with "long double" type, and cannot 
> > find
> anything
> > in the gsl. Has anyone already thought of implementing something like this?
> > 
> > Cheers, Tobias

      __________________________________  Yahoo! Clever: Sie haben Fragen? 
Yahoo! Nutzer antworten Ihnen.

reply via email to

[Prev in Thread] Current Thread [Next in Thread]