Elliptic.jl
Elliptic.jl copied to clipboard
nolta
Elliptic integral and Jacobi elliptic special functions
This PR addresses https://github.com/nolta/Elliptic.jl/issues/28 by removing all type annotations for the elliptic integrals so that they are auto-differentiable. I have checked using Zygote.jl and ForwardDiff.jl and have added a suite...
Is there a reason all the function parameters are strictly typed to `Float64`? This causes issues when using auto-diff tools (which need a more general type: `T
As far as I understand, the elliptic functions in this package are defined only for real values of the first argument ``u``. At the same time, it is good to...
From the documentation of MATLAB on [ellipke](https://www.mathworks.com/help/matlab/ref/ellipke.html#d120e314065) and [ellipj](https://www.mathworks.com/help/matlab/ref/ellipj.html), we see that the inputs to these functions can also be arrays. In the current implementation, we have to use the...
I have seen that the complete elliptic integrals were extended for negative range which is great. Can the non complete elliptic integral be extended as well? In page 593, in...
~~Would be nice to have a catch-all SpecialFunctions package to put (ideally pure Julia) special functions in, maybe this package could be a start?~~ Just noticed that https://github.com/nolta/SpecialFunctions.jl already exists,...
Hi! I was trying to use Elliptic.jl in order to evaluate K(k) to high accuracy using BigFloats, but noticed that if given a BigFloat as an input, `Elliptic.K` returns a...
From @TSGut (CC @MikaelSlevinsky): ```julia julia> α = 4.453667231591123 4.453667231591123 julia> K = Elliptic.K(1-1/α^2) 2.9043516240612575 julia> K2f1 = convert(eltype(α),π)/2*HypergeometricFunctions._₂F₁(one(eltype(α))/2,one(eltype(α))/2,1,1-1/α^2) 2.904351624061259 julia> @test K ≈ K2f1 Test Passed ```
