Skip to content

Commit

Permalink
Fix FGQ. (#6)
Browse files Browse the repository at this point in the history
  • Loading branch information
pkofod authored Nov 15, 2017
1 parent e1b6132 commit 5aff52d
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion src/DistQuads.jl
Original file line number Diff line number Diff line change
Expand Up @@ -17,7 +17,7 @@ distribution(dq::DistQuad) = dq.d

DistQuad(d; N = 32) = DistQuad(d, N)
function DistQuad(d::Distributions.Beta, N)
gj = gaussjacobi(N, d.α-1, d.β-1)
gj = FastGaussQuadrature.JacobiRec(N, d.α-1, d.β-1)
G = gamma(d.α)*gamma(d.β)/gamma(d.α+d.β)
w = gj[2]/((2.0^(d.α+d.β-1.0))*G)
n = (1.-gj[1])./2
Expand Down

0 comments on commit 5aff52d

Please sign in to comment.