#include<algorithm>
#include<sstream>
#include<iostream>
+#include<vector>
#include "utils.h"
#include "my_assert.h"
const LenDist& getGLD() { return *gld; }
- void startSimulation(simul*, double*);
+ void startSimulation(simul*, const std::vector<double>&);
bool simulate(READ_INT_TYPE, PairedEndRead&, int&);
void finishSimulation();
//Use it after function 'read' or 'estimateFromReads'
- double* getMW() {
+ const double* getMW() {
assert(mw != NULL);
return mw;
}
fclose(fo);
}
-void PairedEndModel::startSimulation(simul* sampler, double* theta) {
+void PairedEndModel::startSimulation(simul* sampler, const std::vector<double>& theta) {
this->sampler = sampler;
theta_cdf = new double[M + 1];