Hi, I'm testing the use of QAWO integration to perform my fourier integration because FFT is not suitable as I need more input points (which are computationally had to generate) to get the same accuracy.
In order to see if QAWO does the job, I test my code on the block function of the FFT example in the gsl manual. I then print the frequency space result which looks correct. Then when I try to invert the result to see how close I am to the original block function, I get the gsl error GSL_EDIVERGE. I've checked the code and it seems that I'm following the right approach i.e. call gsl_integration_qawo_table_set in a loop for each omega value then do the integration & save the result for the given omega. Has anyone tested this? The other issue that I have is that I do not fully understand (not clear from the manual) : what n should be set to in: gsl_integration_qawo_table_alloc what limit should be set to in int gsl_integration_qawo in order to get a certain accuracy. Say I have the region from [0, 100] if I think the minimum quadrature interval should be 1e-6, what do I set n and limit to. They seem to be correlated. Regards Serle Shuman ----------------------------------------- Please click on the following link to read the Old Mutual legal notice: <http://www.oldmutual.co.za/e-mail_legal_notice> Alternatively, send a "blank" e-mail to <mailto:[email protected]> _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
