fix a=2k instable case

Lightup1 · Sep 4, 2023 · 0d0fc20 · 0d0fc20
1 parent d9ccca9
commit 0d0fc20
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/char-coeff-new-indexing.jl b/src/char-coeff-new-indexing.jl
@@ -195,7 +195,7 @@ function MathieuExponent(a,q;ndet::Int=20)
     m=Tridiagonal(d[2:N], ones(N), d[1:N-1])
     delta=det(m)
     if x
-        alpha=acos(2delta-1)/pi
+        alpha=acos(2*Complex(delta)-1)/pi
     else
         alpha=2*asin(sqrt(delta*sin(pi*sqrt(Complex(a))/2)^2))/pi
     end