Skip to content

Sttempt to implement elliptic filter

Marco Meyer requested to merge elliptic-filter into master
  • The first attempt implemented (https://people.math.sc.edu/Burkardt/cpp_src/elliptic_integral/elliptic_integral.html), but no special case for (m near 1)
  • The elliptic filter has been implemented based on python scipy ellipk methods
  • The second attempt used CEPHES library routines (calculation of ellipk, ellipkm1 are safe)
  • But there is stil a double precision issue in the calculation of the inverse jacobian (inv_jacobi_sn1 method) when m approches 1.. NB: test example in ./examples/elliptic.cpp

Merge request reports