[teval] Assert s_i == s_f in configuration tree traversal

krivenko · Oct 2, 2024 · 8f2a675 · 8f2a675
1 parent f22dceb
commit 8f2a675
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/src/topology_eval.jl b/src/topology_eval.jl
@@ -459,9 +459,10 @@ function _traverse_configuration_tree!(eval::TopologyEvaluator,
 
     # Are we at a leaf?
     if pos > length(eval.conf)
-        if s_i == s_f # Is the resulting configuration block-diagonal?
-            eval.top_result_mats[s_i] .+= pair_int_weight .* eval!(eval.matrix_prods[s_f])
-        end
+        @assert(s_i == s_f, "A block off-diagonal contribution to a pseudo-particle \
+                propagator detected, try passing more `symmetry_breakers` when \
+                constructing the local exact diagonalization object")
+        eval.top_result_mats[s_i] .+= pair_int_weight .* eval!(eval.matrix_prods[s_f])
         return
     end