MPS_search gives correct energy when max_states < N^2-1

euro-hpc-pl · Feb 23, 2021 · b3540d6 · b3540d6
1 parent 313917c
commit b3540d6
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/test/MPS_search.jl b/test/MPS_search.jl
@@ -107,9 +107,9 @@ states = all_states(get_prop(ig, :rank))
         end
 
         @testset "Results from solve agree with brute-force" begin
-            # The energy is wrong when max_states > N
+            # The energy is wrong when max_states > N^2-2
 
-            for max_states ∈ [1, N, 2*N, 3*N]#, N^2-1, N^2]  
+            for max_states ∈ [1, N, 2*N, 3*N, N^2-3, N^2-2]#, N^2-1, N^2]  
                 states, prob, pCut = solve(rψ, max_states)
                 sp = brute_force(ig, num_states = max_states)