NumericalIntegration.jl
NumericalIntegration.jl copied to clipboard
Basic numerical integration routines for presampled data.
Add support for latest version of Interpolations. Tests seem to run fine.
Added the methods ``` integrate(x::AbstractVector, y::AbstractVector, ::Simpson) integrate(x::AbstractVector, y::AbstractVector, ::SimpsonFast) ``` which work on irregular grids. Also added a `testset` to test the new methods (together with Trapezoidal/Trapezoidal Fast) Example:...
```julia julia> x = collect(1.0:5.0) 5-element Vector{Float64}: 1.0 2.0 3.0 4.0 5.0 julia> y = x.^2 5-element Vector{Float64}: 1.0 4.0 9.0 16.0 25.0 julia> cumul_integrate(x, y) 5-element Vector{Float64}: 2.5 5.0...
I am trying to compute the numerical integral from a time series. This fails due to a type conversion error that I do not really understand. ``` julia> using Dates,...
When I use culum_integrate, it appears that it integrates 1 index beyond the array. ``` using NumericalIntegration a = LinRange(0,1,100) print(integrate(a,a), ' ') print(last(cumul_integrate(a,a))) ``` returns: ``` 0.5 0.5000510152025299 ```...
I posted an implementation on discourse (https://discourse.julialang.org/t/romberg-jl-simple-romberg-integration-in-julia/39313/8?u=stevengj) using Richardson.jl that is much more general (is not restricted to power-of-two sizes), seems to be more accurate than your current code, and...
Hi, Is there a reason `dims` is ordered differently than the multi-dimensional arrays? I have to do: ```julia julia> x=1:3 julia> y=1:5 julia> Z = randn(3, 5) julia> integrate(x, Z;...
Many state-of-the-art Julia Frameworks rely on Interpolations.jl v0.14, but NumericalIntegration.jl forces v0.13. This results in many libraries aren't usabel together with NumericalIntegrations.jl. Would be nice to add compat entry like:...
Since version 0.14.2 the Interpolations package changed some of its notation from CamelCase to snake_case in accordance to Julia's recommendations. I simply adjusted the package so that it is now...
Now that the package has been moved to JuliaMath, would it be possible to register the current version `v0.3.4`? It solves the long-standing issue with `Interpolations.jl` compatibility - see https://github.com/JuliaMath/NumericalIntegration.jl/issues/47...