Output theta and model even when no alignable reads

deweylab · Apr 6, 2019 · e6e14db · e6e14db
1 parent aa8bbc3
commit e6e14db
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/EM.cpp b/EM.cpp
@@ -615,6 +615,13 @@ int main(int argc, char* argv[]) {
 	if (N1 == 0) {
 		printf("Warning: There are no alignable reads!\n");
 		theta.resize(M + 1, 0.0);
+		FILE *fo = NULL;
+		sprintf(thetaF, "%s.theta", statName);
+		fo = fopen(thetaF, "w");
+		fclose(fo);
+		sprintf(modelF, "%s.model", statName);
+		fo = fopen(modelF, "w");
+		fclose(fo);
 		eel.resize(M + 1, 0.0);
 		for (int i = 1; i <= M; ++i) eel[i] = transcripts.getTranscriptAt(i).getLength();
 		double *countv = new double[M + 1];