我正在使用库中实现的 Vegas（重要性采样）蒙特卡洛方法gsl来评估形式的各种积分I_i = int f(x) g_i(x) dx，不同g_i但相同f > 0。特别是，其中一个g_i，说g_0是相同的1，所以它对应的积分就是 的积分f。可以将另一个视为分布下I_i的平均值，更明确地说：g_if
I_i/I_0 = (int f(x) g_i(x) dx)/(int f(x) dx) = <g_i>_f
由于重要性采样背后的想法是根据其绝对值对被积函数进行采样，我想知道是否有可能（并且有效）利用积分后获得的已经优化的网格，f并从该分布中采样所有的g_i，通过设置选项stage=3并将调用中的被积函数更改为gsl_monte_vegas_integrate.