pro->read(fi);
npro->read(fi);
- if (fscanf(fi, "%d", &M) == 1) {
- mw = new double[M + 1];
- for (int i = 0; i <= M; i++) fscanf(fi, "%lf", &mw[i]);
+ if (fscanf(fi, "%d", &val) == 1) {
+ if (M == 0) M = val;
+ if (M == val) {
+ mw = new double[M + 1];
+ for (int i = 0; i <= M; i++) fscanf(fi, "%lf", &mw[i]);
+ }
}
fclose(fi);
}
-//Only master node can call
+//Only master node can call. Only be called at EM.cpp
void PairedEndModel::write(const char* outF) {
FILE *fo = fopen(outF, "w");
qpro->read(fi);
nqpro->read(fi);
- if (fscanf(fi, "%d", &M) == 1) {
- mw = new double[M + 1];
- for (int i = 0; i <= M; i++) fscanf(fi, "%lf", &mw[i]);
+ if (fscanf(fi, "%d", &val) == 1) {
+ if (M == 0) M = val;
+ if (M == val) {
+ mw = new double[M + 1];
+ for (int i = 0; i <= M; i++) fscanf(fi, "%lf", &mw[i]);
+ }
}
+
fclose(fi);
}
-//Only master node can call
+//Only master node can call. Only be called at EM.cpp
void PairedEndQModel::write(const char* outF) {
FILE *fo = fopen(outF, "w");
npro->collect(*(o.npro));
}
+//Only master node can call
void SingleModel::read(const char* inpF) {
int val;
FILE *fi = fopen(inpF, "r");
pro->read(fi);
npro->read(fi);
- if (fscanf(fi, "%d", &M) == 1) {
- mw = new double[M + 1];
- for (int i = 0; i <= M; i++) fscanf(fi, "%lf", &mw[i]);
+ if (fscanf(fi, "%d", &val) == 1) {
+ if (M == 0) M = val;
+ if (M == val) {
+ mw = new double[M + 1];
+ for (int i = 0; i <= M; i++) fscanf(fi, "%lf", &mw[i]);
+ }
}
fclose(fi);
}
+//Only master node can call. Only be called at EM.cpp
void SingleModel::write(const char* outF) {
FILE *fo = fopen(outF, "w");
qpro->read(fi);
nqpro->read(fi);
- if (fscanf(fi, "%d", &M) == 1) {
- mw = new double[M + 1];
- for (int i = 0; i <= M; i++) fscanf(fi, "%lf", &mw[i]);
+ if (fscanf(fi, "%d", &val) == 1) {
+ if (M == 0) M = val;
+ if (M == val) {
+ mw = new double[M + 1];
+ for (int i = 0; i <= M; i++) fscanf(fi, "%lf", &mw[i]);
+ }
}
fclose(fi);
}
-//Only master node can call
+//Only master node can call. Only be called at EM.cpp
void SingleQModel::write(const char* outF) {
FILE *fo = fopen(outF, "w");
for (int i = 1; i <= M; i++) {
getline(fin, line);
size_t pos = 0;
- for (int j = 0; j < 3; j++) pos = line.find_first_of('\t', pos) + 1;
+ for (int j = 0; j < 2; j++) pos = line.find_first_of('\t', pos) + 1;
size_t pos2 = line.find_first_of('\t', pos);
if (pos2 == string::npos) pos2 = line.length();
tau = atof(line.substr(pos, pos2 - pos).c_str());