fix npdm mpo simpl

block-hczhai · Jun 24, 2024 · 6ae27a9 · 6ae27a9
1 parent 81e3c11
commit 6ae27a9
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/dmrg/mpo_simplification.hpp b/src/dmrg/mpo_simplification.hpp
@@ -335,6 +335,8 @@ template <typename S, typename FL> struct SimplifiedMPO : MPO<S, FL> {
                 set<shared_ptr<OpExpr<S>>, op_expr_less<S>> left_zero_ops,
                     right_zero_ops;
                 MPO<S, FL>::load_left_operators(i);
+                if (MPO<S, FL>::tag != mpo->tag)
+                    mpo->load_left_operators(i);
                 for (size_t j = 0;
                      j < MPO<S, FL>::left_operator_names[i]->data.size(); j++)
                     if (MPO<S, FL>::left_operator_names[i]
@@ -345,6 +347,8 @@ template <typename S, typename FL> struct SimplifiedMPO : MPO<S, FL> {
                 mpo->unload_left_operators(i);
                 MPO<S, FL>::unload_left_operators(i);
                 MPO<S, FL>::load_right_operators(i + 1);
+                if (MPO<S, FL>::tag != mpo->tag)
+                    mpo->load_right_operators(i + 1);
                 for (size_t j = 0;
                      j < MPO<S, FL>::right_operator_names[i + 1]->data.size();
                      j++)