agerlach
agerlach
For `ODEProblem`, parameters are passed by reference, so parameter updates via callbacks can cause upstream issues. This occurs for both `VectorContinuousCallback` and `ContinuousCallback`. As expected, the `ODEProblem` object is also...
Interpolation of unitful types works for ℝ -> ℝ but fails for ℝ -> ℝⁿ ``` using Interpolations, Unitful x = collect(1:10) y = [[i, 2i] for i in x]...
I would like to use VexCL w/ ODEINT to do "event detection" or for finding "zero crossings", but I am struggling to determine the best way to do this. I...
Updating`transformation_if_inf` is sufficient to make HCubatureJL() type stable. Other methods are still type unstable.
Currently not working. See various commented out tests, e.g. https://github.com/SciML/Quadrature.jl/blob/2fa4a7d10e213a1206f759310264b411b304e731/test/derivative_tests.jl#L153
```julia using Quadrature function f(x,p,c) c[1]+=1 sum(sin.(x)) end algs = [HCubatureJL(), CubatureJLh(), CubatureJLp(), CubaSUAVE(), CubaDivonne(), CubaCuhre()] count_array = zeros(length(algs)) for (i,alg) ∈ enumerate(algs) try count = [0] prob = QuadratureProblem((x,p)->f(x,p,count),ones(2),3ones(2))...
I wasn't aware of this package previously. provides and interesting interface. https://github.com/JuliaApproximation/DomainIntegrals.jl
Related to #35 This initial PR commit introduces a `f(x,p)=Inf` integrand test for all quadalgs.
Many of the quadrature algorithms do not handle `Inf` in the integrand well. `QuadGK` throws a `DomainError`, 2D or higher HCubature returns `NaN`. Cubature has similar issues but at different...
For DiffEqUncertainty `expectation()` applications, distributions that are narrow relative to the support can lead to incorrect results via Quadrature as the integrand is not sampled at points w/ non-zero (numerically)...