+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");
+}
+