Cuba.jl
Cuba.jl copied to clipboard
giordano
Library for multidimensional numerical integration with four independent algorithms: Vegas, Suave, Divonne, and Cuhre.
Hi, Thank you for creating this fantastic package! I use it heavily for calculating Feynman diagrams in condensed matter physics. For my applications, I need to provide the integral with...
i would like to build a histogram of some distribution while integrating it with vegas. Is there a way to access the "vegas-weight" of a point at each step in...
Hello, the docs state: > nvec, minevals, maxevals, nnew, nmin, neval are passed to the Cuba library as 64-bit integers, so they are limited to be at most typemax(Int64) However,...
Hi Giordano, I found that when the number of components of f is greater than 1024, cuhre fails to give the correct answer. Is this the limitation of the algorithm?...
I'm not sure it's actually feasible, but it would be great if `Cuba.jl` could take advantage of parallelization capability of Cuba Library. Concurrency is achieved using `fork` and `wait`, but...
Cuba produces a lot of allocations, 4 per iteration in the example below. Is there a way to avoid this? ```julia julia> using Cuba julia> f!(x,ret ) = ret f!...
In the Mathematica interface to Cuba, I can obtain a list of the regions used by option Regions -> True. Is there any way to get this information in Julia?...
This gives a segmentation fault error: using Cuba vegas((x, f) -> f[1] = cos(x[1]), maxevals=100) It works fine at 1000, and also this works fine: vegas((x, f) -> f[1] =...
I think the interface would be much more intuitive if the users can directly specify the integration range, e.g. cuhre(f::Function, xmin, xmax), similar as the interface in Cubature.jl
On page 14 of the Cuba paper (https://arxiv.org/abs/hep-ph/0404043) there is discussion about having a phase argument on the integrand so that an approximation with similar peak structure can be passed...
