jacobi-theta 0.1.1.1 → 0.1.2.0
raw patch · 3 files changed
+15/−7 lines, 3 filesPVP ok
version bump matches the API change (PVP)
API changes (from Hackage documentation)
Files
- CHANGELOG.md +5/−0
- jacobi-theta.cabal +1/−1
- src/Math/JacobiTheta.hs +9/−6
CHANGELOG.md view
@@ -1,3 +1,8 @@+## 0.1.2.0 - 2023-03-02++Some values of the Jacobi theta functions were wrong.++ ## 0.1.1.1 - 2023-03-01 Removed a useless type alias.
jacobi-theta.cabal view
@@ -1,5 +1,5 @@ name: jacobi-theta-version: 0.1.1.1+version: 0.1.2.0 synopsis: Jacobi Theta Functions description: Evaluation of the Jacobi theta functions. homepage: https://github.com/stla/jacobi-theta#readme
src/Math/JacobiTheta.hs view
@@ -64,18 +64,21 @@ termp = altp * exp (- square (nplus' - 0.5 + z') / t') newseries = series + termm + termp -falpha :: Cplx -> Cplx -> Cplx-falpha z tau = - sqrt (-i_ * tau) * exp (i_ / tau * z * z / pi)+--falpha :: Cplx -> Cplx -> Cplx+--falpha z tau = +-- sqrt (-i_ * tau) * exp (i_ / tau * z * z / pi) jtheta1Alt :: Cplx -> Cplx -> Cplx jtheta1Alt z tau = if imagPart tau > 1.3 then- let w = pi * tau in - i_ * jtheta1Alt2 (z / w) (i_ / w) / falpha z tau+ jtheta1Alt1 z (exp (i_ * pi * tau))+-- let w = pi * tau in +-- i_ * jtheta1Alt2 (z / w) (i_ / w) / falpha z tau else- i_ * jtheta1Alt1 (z / tau) (exp (-i_ * pi / tau)) / falpha z tau+ let t = - i_ * tau / pi in+ jtheta1Alt2 (z / pi) t+-- i_ * jtheta1Alt1 (z / tau) (exp (-i_ * pi / tau)) / falpha z tau tauFromQ :: Cplx -> Cplx tauFromQ q = -i_ * log q / pi