26 #include "boost/random.hpp"
31 using namespace pwiz::proteome;
32 using namespace pwiz::analysis;
33 using namespace pwiz::util;
44 for (
size_t i=0; i < pl->size(); ++i)
47 ProteinPtr decoy = decoyList.protein(i + pl->size());
49 if (
os_) *
os_ << target->id <<
" " << target->sequence() << endl;
50 if (
os_) *
os_ << decoy->id <<
" " << decoy->sequence() << endl;
54 unit_assert(
string(target->sequence().rbegin(), target->sequence().rend()) == decoy->sequence());
65 boost::mt19937 engine(0);
66 boost::uniform_int<> distribution;
67 boost::variate_generator<boost::mt19937, boost::uniform_int<> > rng(engine, distribution);
69 for (
size_t i=0; i < pl->size(); ++i)
72 ProteinPtr decoy = decoyList.protein(i + pl->size());
74 if (
os_) *
os_ << target->id <<
" " << target->sequence() << endl;
75 if (
os_) *
os_ << decoy->id <<
" " << decoy->sequence() << endl;
79 string sequence = target->sequence();
80 random_shuffle(sequence.begin(), sequence.end(), rng);
96 int main(
int argc,
char* argv[])
102 if (argc>1 && !strcmp(argv[1],
"-v"))
os_ = &cout;