Rhys,
Unfortunately our 'test' code is 400-500 lines. We may be able to reduce
it in the coming weeks. Apart from all the bells and whistles, the
original code allocated a single workspace and then used that in a loop to
compute billions of integrals (calls to gsl_integration_qags()) with
stochastic parameters. The results were really strange (and wrong) until
we changed the code to allocate a workspace for every single integration.
I.e., all we did for the fix was adding one line before and after each
gsl_integration_qags() call to allocate/free a workspace (maybe a dozen
lines in total).
Denes
Email from Rhys Ulerich on Dec 1, 2011 at 9:59am:
Hi Denes,
On Wed, Nov 30, 2011 at 3:09 PM, Denes Molnar <mol...@physics.purdue.edu> wrote:
The documentation does not discuss whether gsl_integration_workspace is
reusable. Is it? For example, can one allocate a single workspace, and use
that to do two different integrations afterwards? Some of our tests with
gsl-1.14 indicate that the answer is 'no'.
Do you have some test code demonstrating that 'no' result?
- Rhys
_______________________________________________
Help-gsl mailing list
Help-gsl@gnu.org
https://lists.gnu.org/mailman/listinfo/help-gsl