Computation of elliptic integrals with arbitrary precision

I have seen the issue Support for extended precision in Elliptic.jl. It references the pull request Jacobi elliptic functions and complete elliptic integral of the first kind in SpecialFunctions.jl which is year old and is still opened (could anyone explain why is it so?).

Besides that, I was not able to find anything.
So, is there a julian realisation of Elliptic Integrals with arbitrary precision?

If there is nothing like that, can someone advise a library in other language which would implement this feature and which would be efficient to call from Julia?

ArbNumerics.jl works.
see Elliptic Integrals (complete, incomplete, symmetric or Carlson).