bug fix.

wztzjhn · Dec 1, 2017 · ae5a5e8 · ae5a5e8
1 parent 8267e1e
commit ae5a5e8
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/model.cc b/src/model.cc
@@ -1850,7 +1850,7 @@ namespace qbasis {
             std::vector<std::pair<MKL_INT, T>> values;
             T value_gs = static_cast<T>(0.0);
             for (auto it = Bq.mats.begin(); it != Bq.mats.end(); it++) {
-                if (Bq.q_diagonal()) {                                           // only momentum changes
+                if (it->q_diagonal()) {                                           // only momentum changes
                     values.push_back(std::pair<MKL_INT, T>(j, sj * basis_vrnl[sec_old][j].diagonal_operator(props,*it)));
                 } else {
                     intermediate_states[tid].copy(basis_vrnl[sec_old][j]);