* 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA *
***************************************************************************/
-#include "geneorder.h"
-
+#include <algorithm>
+#include <stdexcept>
+#include <cstdlib>
using namespace std;
+#include "geneorder.h"
+
/**
* \file GeneOrder.cpp
* Implements the GeneOrder class
*/
-
GeneOrder::GeneOrder(const GeneOrder& go){
_geneorder = go._geneorder;
}
+void GeneOrder::pad(){
+ if (_geneorder[0] != 0)
+ _geneorder.insert(_geneorder.begin(),0);
+ if(_geneorder.back() != static_cast<Gene>(_geneorder.size() - 1))
+ _geneorder.push_back(_geneorder.size());
+}
+
+struct Abs{
+ Gene operator()(Gene x) const{
+ return abs(x);
+ }
+};
+
+void GeneOrder::verify(){
+ if(_geneorder[0] != 0)
+ throw invalid_argument("Permutation doesn't start with 0");
+ if ( _geneorder.back() != static_cast<Gene>(_geneorder.size() - 1))
+ throw invalid_argument("Permutation doesn't end with n+1");
+
+ GeneList genes(_geneorder);
+ transform(genes.begin(),genes.end(),genes.begin(),Abs());
+ sort(genes.begin(),genes.end());
+ if (unique(genes.begin(),genes.end()) != genes.end())
+ throw std::invalid_argument("Not all genes are present exactly 1 time");
+}
+
GeneOrder::~GeneOrder(){
}
*/
const Gene& GeneOrder::operator[](size_type i) const{
if (i < 0 || i >= _geneorder.size())
- throw out_of_range("Index is not in valid range");
+ throw out_of_range("Indexis not in valid range");
return _geneorder[i];
}
-int GeneOrder::size() const{
+GeneOrder::size_type GeneOrder::size() const{
return _geneorder.size();
}
const GeneOrder::GeneList& GeneOrder::list() const{
return _geneorder;
}
+
+
+GeneOrder::iterator GeneOrder::begin() const{
+ return _geneorder.begin();
+}
+
+GeneOrder::iterator GeneOrder::end() const{
+ return _geneorder.end();
+}