We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
gaussjacobi
α, β
This works:
julia> gaussjacobi(500, 0.0, -1/2);
but changing 0.0 to an Int gives an error, probably indicating a missing call to float(α) or similar:
0.0
Int
float(α)
julia> gaussjacobi(500, 0, -1/2) ERROR: MethodError: no method matching asy1(::Int64, ::Int64, ::Float64, ::Int64) Closest candidates are: asy1(::Integer, ::Float64, ::Float64, ::Integer) at ~/.julia/packages/FastGaussQuadrature/BRLTf/src/gaussjacobi.jl:200 Stacktrace: [1] jacobi_asy(n::Int64, α::Int64, β::Float64) @ FastGaussQuadrature ~/.julia/packages/FastGaussQuadrature/BRLTf/src/gaussjacobi.jl:183 [2] gaussjacobi(n::Int64, α::Int64, β::Float64) @ FastGaussQuadrature ~/.julia/packages/FastGaussQuadrature/BRLTf/src/gaussjacobi.jl:45 [3] top-level scope @ REPL[98]:1
The text was updated successfully, but these errors were encountered:
jacobi_asy should be restricted to Float64 as the choice of boundary terms is precision dependent:
jacobi_asy
Float64
FastGaussQuadrature.jl/src/gaussjacobi.jl
Line 175 in 8b964b7
It would be easy to add a promote for the case of differing types of course.
promote
Sorry, something went wrong.
Maybe promote(float(α), float(β))... since passing integers should be allowed.
promote(float(α), float(β))...
Successfully merging a pull request may close this issue.
This works:
but changing
0.0
to anInt
gives an error, probably indicating a missing call tofloat(α)
or similar:The text was updated successfully, but these errors were encountered: