30 using namespace pwiz::math;
31 using namespace pwiz::util;
32 using namespace pwiz::analysis;
33 using namespace pwiz::data::peakdata;
41 if (
os_) *
os_ <<
"testParabola()\n";
43 const double center = 2.1;
44 const double height = 5;
46 vector<OrderedPair> pairs;
47 for (
double i=0; i<5; i++)
48 pairs.push_back(
OrderedPair(i, height-(i-center)*(i-center)));
61 copy(peak.data.begin(), peak.data.end(), ostream_iterator<OrderedPair>(*
os_,
" "));
77 if (
os_) *
os_ <<
"testMultiplePeaks()\n";
79 const double center = 2.1;
80 const double height = 5;
82 vector<OrderedPair> pairs;
83 for (
double i=0; i<5; i++)
84 pairs.push_back(
OrderedPair(i, height-(i-center)*(i-center)));
85 for (
double i=0; i<5; i++)
86 pairs.push_back(
OrderedPair(i+5, height-(i-center)*(i-center)));
87 for (
double i=0; i<5; i++)
88 pairs.push_back(
OrderedPair(i+10, height-(i-center)*(i-center)));
90 vector<size_t> indices;
93 indices.push_back(12);
98 fitter.
fitPeaks(pairs, indices, peaks);
102 for (vector<Peak>::const_iterator it=peaks.begin(); it!=peaks.end(); ++it)
105 copy(it->data.begin(), it->data.end(), ostream_iterator<OrderedPair>(*
os_,
" "));
125 int main(
int argc,
char* argv[])
131 if (argc>1 && !strcmp(argv[1],
"-v"))
os_ = &cout;