101 ame solver gets wrong result for constant hamiltonian (#102)

* 🐛 add a 2π factor for in AME solver for constant Hamiltonian * 🔖 bump version to v0.7.4
USCqserver · Jan 14, 2023 · 2f66aad · 2f66aad
1 parent 06148c0
commit 2f66aad
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "OpenQuantumTools"
 uuid = "e429f160-8886-11e9-20cb-0dbe84e78965"
 authors = ["Huo Chen <[email protected]>"]
-version = "0.7.3"
+version = "0.7.4"
 
 [deps]
 DiffEqCallbacks = "459566f4-90b8-5000-8ac3-15dfb0a30def"

diff --git a/src/QSolver/ame_solver.jl b/src/QSolver/ame_solver.jl
@@ -128,7 +128,7 @@ function solve_ame(
     H = w |> Diagonal |> sparse |> Hamiltonian
     inters = rotate(A.interactions, v)
     # build gaps for AME
-    gap_idx = OpenQuantumBase.build_gap_indices(w, digits, sigdigits, cutoff, lvl)
+    gap_idx = OpenQuantumBase.build_gap_indices(2*π*w, digits, sigdigits, cutoff, lvl)
     # prepare for the initial state
     u0 = build_u0(A.u0, :m, vectorize=vectorize)
     u0 = v' * u0 * v
@@ -170,7 +170,7 @@ function build_ensemble_ame(
     H = w |> Diagonal |> sparse |> Hamiltonian
     inters = rotate(A.interactions, v)
     # build gaps for AME
-    gap_idx = OpenQuantumBase.build_gap_indices(w, digits, sigdigits, cutoff, lvl)
+    gap_idx = OpenQuantumBase.build_gap_indices(2*π*w, digits, sigdigits, cutoff, lvl)
     # prepare for the initial state
     u0 = build_u0(A.u0, :v)
     u0 = v' * u0