From: Dariusz Murakowski <murakdar@mit.edu>
Date: Tue, 23 Jun 2015 06:10:53 +0000 (-0400)
Subject: Define state2aa properly.
X-Git-Url: http://src.murakowski.org/?a=commitdiff_plain;h=7de0f423904530a45fb762c001404b3d83d5aa23;p=VirEvoDyn.git

Define state2aa properly.
---

diff --git a/ham_ss.cpp b/ham_ss.cpp
index df392f3..338d454 100644
--- a/ham_ss.cpp
+++ b/ham_ss.cpp
@@ -205,9 +205,10 @@ void PottsHamiltonian::getSeq2State(const std::string &FILENAME)
     
     for (size_t i=0; i<Schar.size(); i++) {
         state2aa.push_back(std::vector<aa>(Schar[i].size()));
-        std::vector<aa> last = state2aa.back();
+        std::vector<aa> &last = state2aa.back();
         for (size_t j=0; j<Schar[i].size(); j++) {
-            last.push_back(str2aa.at(Schar[i][j]));
+            last[j] = str2aa.at(Schar[i][j]);
+            // could have used reserve() above and then push_back() here but eh
         }
         // move the first element to last place, because that is how states are ordered
         std::rotate(last.begin(),++last.begin(),last.end());    // could have used std::next() but eh