HCubature.jl
HCubature.jl copied to clipboard
JuliaMath
pure-Julia multidimensional h-adaptive integration
Hi, I made some benchmarks and noted that for 1-d integral `quadgk` is way more fast than `hquadrature`. Here is a simple example ```julia julia> using HCubature, BenchmarkTools julia> f()...
Dear @stevengj I would like to ask you a question when computing `∫∫g(x)*g(y)f(x,y)dxdy`. I hope it is fine to post my question here. Below is my problem. ``` using QuadGK,...
I started to put a PR together to allow for convergence testing based on the individual integrands, like is default in `Cubature.jl`. However, I quickly realized it may get more...
I am integrating `L(f) = \int K.(f .- g(x)) dx` where f is a vector of 100-1000 elements, `g` and `K` are scalar, the integration is in R^3. The dimension...
I am looking at some SciML problems where, depending on the stability of the ODE, the integrand used w/ HCubature.jl may return Inf. I am noticing inconsistent behavior across `hquadrature`...
I need a pure-julia integration library that support multi-dimensional, vector valued integrands **and** batch integrand evaluation. HCubature.jl hits all these criteria except for batch integrand evaluation as in [Cubature.jl](https://github.com/JuliaMath/Cubature.jl). Is...
Hi, I noticed some weird behaviour while using `hquadrature`: `hquadrature(t->1, 0, -1)[1]` outputs 1, while, mathematically speaking, it should be -1. Is this intended, or does anybody have an idea...
One of the code using HCubature.jl got internal error on Julia1.2.0 . Here is an example how to reproduce the error: ``` $ docker run --rm -it julia:1.2.0 \ julia...
In the previous Cubature.jl package, it is possible to evaluate the integrand in parallel by requiring the function `f` to accept matrix arguments. Is it possible such feature can be...
