int main(int argc, char* argv[]) {
string datafile = argv[1];
string treefile = argv[2];
string name = argv[3];
// MGModel* model = new MGModel(datafile,treefile);
GamMixModel* model = new GamMixModel(datafile,treefile);
cerr << "start\n";
ofstream os((name + ".trace").c_str());
ofstream ros((name + ".rate").c_str());
model->TraceHeader(os);
while (1) {
model->Move(1);
model->Move(0.1);
model->Move(0.01);
model->Trace(os);
model->rate->ToStream(ros);
ros << '\n';
ros.flush();
ofstream mon_os((name + ".monitor").c_str());
model->Monitor(mon_os);
}
}
virtual void write(std::ostream &os) const override
{
llvm::raw_os_ostream ros(os);
write(ros);
}
