Skip to content

Commit

Permalink
up
Browse files Browse the repository at this point in the history
  • Loading branch information
tinatorabi authored Sep 10, 2024
1 parent 2ffc4b2 commit 511ea5e
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion src/ellipjc.jl
Original file line number Diff line number Diff line change
Expand Up @@ -44,7 +44,7 @@ function ellipjc(u, L; flag=false)
else
kappa = ComplexElliptic.polyval(Complex{Float64}.([132.0, 42.0, 14.0, 5.0, 2.0, 1.0, 0.0]), Complex{Float64}(m / 4.0))
end
sn1, cn1, dn1 = Zygcompat_ellipjc(u / (1 + kappa), kappa ^ 2, flag=true)
sn1, cn1, dn1 = ellipjc(u / (1 + kappa), kappa ^ 2, flag=true)

denom = 1 .+ kappa .* sn1 .^ 2
sn = (1 .+ kappa) .* sn1 ./ denom
Expand Down

0 comments on commit 511ea5e

Please sign in to comment.