对于指数函数,GSL qag_integration非常慢

时间:2013-09-25 00:08:37

标签: c++ gsl

这个问题跟我上一个问题有关GSL integration behaves strange

由于缩放x -> (1-t)/t在无限积分方法gsl_integration_qagi中导致不受欢迎的答案,我现在使用有限支持上的集成。我现在要做的是:

我有一个(离散的)实数的S系列,我用给定的指数函数卷积(每个)

    exp(-t/T) ... T = decay constant

我选择对集成的支持为(min(S) - 10*T, max(S) + 10*T),以便我涵盖函数的大部分“重要”贡献。

使用gsl_integration_qag对此支持进行集成需要几秒钟,而gsl_integration_qagi`(无限支持)几乎不会产生几毫秒但会产生错误的结果。有谁知道原因吗?

如果卷积是高斯而不是指数,则相同的gsl_integration_qag效果很好。

提前致谢, NIKHIL

0 个答案:

没有答案