[Top][All Lists]

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

Re: [Help-gsl] gsl_integration_qawf

From: Patrick Alken
Subject: Re: [Help-gsl] gsl_integration_qawf
Date: Mon, 15 Jan 2018 11:29:20 -0700
User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:52.0) Gecko/20100101 Thunderbird/52.5.0


  A return code of 0 is equivalent to GSL_SUCCESS, so yes that is good.
I'm not sure what you mean by inner/outer integration - can you provide
a minimal working source code to demonstrate the problem you're having?

  The adaptive integration routines in GSL are based on QUADPACK - you
can find further documentation on QAWF here:
There is also an example calling program for QAWF on that page (in
fortran), perhaps it will help you.

  Thanks for the whitespace issue - I have fixed it on the git.


On 01/15/2018 01:35 AM, Richard Hartmann wrote:
> Hi
> I'm struggling with the numeric Fourier integration, in particular
> with the error handling of the 'gsl_integration_qawf' routine.
> As of the nested integration and if the error handling is default, my
> test aborts since the inner integration reports a roundoff error.
> (current Git 3bbe65ad5a07362e82103b311feda55d4355a2fe, qawf.c l. 133)
> So far so good, however, if the error handling is off, the outer
> integration seem to run through. The overall return code is 0,
> indicating success. Is that intensional?
> I tried to answer the question my self, but I could not understand
> what the 'gsl_integration_qawf' routine was doing.
> Secondly and very minor, there is a 'space' missing in the error
> message at qawo.c l. 104
>       GSL_ERROR ("cannot reach tolerance because of roundoff error"
>                  "on first attempt", GSL_EROUND);
> should be
>       GSL_ERROR ("cannot reach tolerance because of roundoff error "
>                  "on first attempt", GSL_EROUND);
> Thanks for this great library and thanks for you help.
> Richard

reply via email to

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