[Top][All Lists]

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

Re: [Help-gsl] Integration in a long while loop

From: Tommy Nordgren
Subject: Re: [Help-gsl] Integration in a long while loop
Date: Tue, 19 Feb 2008 21:32:30 +0100

On 16 feb 2008, at 19.32, Sam M. Vinko wrote:

I have noticed that using the gsl integration in a loop will work only a certain number of times. If I take an integral functional and try to integrate it in a loop (for a certain set of values of a parameter A - say from 1 to 1000) after some values are calculated (for example from A=1 to A=10) the error

gsl: qags.c:558: ERROR: roundoff error detected in the extrapolation table
Default GSL error handler invoked.
Abort trap

shows up. I thought that maybe the problem was with the integral for A = 10 or A = 11, but if I now start the same code at A = 10 I get all the numbers up to A = 20, and again the error. This makes me think this is a memory issue, so here is my question - is there any other way to clear out work space memory than with

gsl_integration_workspace_free ()

? Changing the size of workspace_alloc does not seem to do anything Is there something I'm missing?

Did you forget to re-initialize the gsl data structures used for integration at each new value of A?

Skinheads are so tired of immigration, that they are going to move to a country that don't accept immigrants!
Tommy Nordgren

reply via email to

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