Comments (4)
Thanks for the bug report! I made a PR to fix this bug.
from fastgaussquadrature.jl.
Thanks for the fast fix! Could you please be so kind to tag a new release with new fix?
This bug breaks our package for Bayesian inference on precompilation stage, because apparently we used 21 number of points as a default value 😅😅.
from fastgaussquadrature.jl.
done!
from fastgaussquadrature.jl.
This bug was introduced in release v0.4.8
.
julia> using FastGaussQuadrature # On v0.4.8
julia> gausshermite(21)
ERROR: BoundsError: attempt to access 10-element Vector{Float64} at index [11]
Stacktrace:
[1] setindex!
@ ./array.jl:843 [inlined]
[2] setindex!
@ ./array.jl:857 [inlined]
[3] hermite_initialguess(n::Int64)
@ FastGaussQuadrature ~/.julia/dev/FastGaussQuadrature/src/gausshermite.jl:271
[4] hermite_rec(n::Int64)
@ FastGaussQuadrature ~/.julia/dev/FastGaussQuadrature/src/gausshermite.jl:86
[5] unweightedgausshermite(n::Int64)
@ FastGaussQuadrature ~/.julia/dev/FastGaussQuadrature/src/gausshermite.jl:41
[6] gausshermite(n::Int64)
@ FastGaussQuadrature ~/.julia/dev/FastGaussQuadrature/src/gausshermite.jl:23
[7] top-level scope
@ REPL[2]:1
julia> using FastGaussQuadrature # On v0.4.7
julia> gausshermite(21)
([-5.550351873264678, -4.773992343411219, -4.1219955474918395, -3.5319728771376777, -2.979991207704598, -2.453552124512838, -1.9449629491862535, -1.448934250650732, -0.9614996344183689, -0.4794507070791075 … 0.4794507070791075, 0.9614996344183689, 1.448934250650732, 1.9449629491862535, 2.453552124512838, 2.979991207704598, 3.5319728771376777, 4.1219955474918395, 4.773992343411219, 5.550351873264678], [3.7203650688183795e-14, 8.818611242049893e-11, 2.5712301800593306e-8, 2.1718848980566762e-6, 7.478398867310059e-5, 0.0012549820417264083, 0.011414065837434388, 0.060179646658912324, 0.1921203240669979, 0.381669073613502 … 0.381669073613502, 0.1921203240669979, 0.060179646658912324, 0.011414065837434388, 0.0012549820417264083, 7.478398867310059e-5, 2.1718848980566762e-6, 2.5712301800593306e-8, 8.818611242049893e-11, 3.7203650688183795e-14])
from fastgaussquadrature.jl.
Related Issues (20)
- docstring format is not unified
- Question on unweighted Gauss-Hermite HOT 1
- gausslaguerre(39,44.8) doesn't converge, document these conditions for convergence
- Gauss-Chebyshev documentation has incorrect weights for third and fourth kinds
- Add mutating versions of functions that reuse external storage for sigma points and weights HOT 5
- missing promotion in `gaussjacobi` for integer `α, β` HOT 2
- Gauss–Kronrod points and weights? HOT 5
- export/document 3-term recurrence coefficients HOT 3
- The precision of `gausslegendre` can be improved HOT 1
- PLEASE DO NOT COMMIT DIRECTLY ON THE `master` BRANCH! HOT 3
- Road to v2.0.0 HOT 3
- Loss of accuracy in weights of large Gauss-Jacobi rules HOT 9
- Migrate to travis-ci.com HOT 3
- TagBot trigger issue HOT 12
- `gausslaguerre(n, α)` produces error if α is Int HOT 2
- `gaussjacobi( 100000, .9, -.1 )` takes much more time than before
- Piessens's Chebyshev series approximations doesn't have enough accurate HOT 4
- Time to switch to GitHub Actions? HOT 18
- Quadrature rules for the weight function w(x) = 1/(1+x^2) HOT 8
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from fastgaussquadrature.jl.