| 58 | simul *sampler; |
| 59 | |
| 60 | void genOutReadStreams(int type, char *outFN) { |
| 61 | switch(type) { |
| 62 | case 0 : |
| 63 | n_os = 1; |
| 64 | sprintf(outReadF[0], "%s.fa", outFN); |
| 65 | break; |
| 66 | case 1 : |
| 67 | n_os = 1; |
| 68 | sprintf(outReadF[0], "%s.fq", outFN); |
| 69 | break; |
| 70 | case 2 : |
| 71 | n_os = 2; |
| 72 | for (int i = 0; i < n_os; i++) |
| 73 | sprintf(outReadF[i], "%s_%d.fa", outFN, i + 1); |
| 74 | break; |
| 75 | case 3 : |
| 76 | n_os = 2; |
| 77 | for (int i = 0; i < n_os; i++) |
| 78 | sprintf(outReadF[i], "%s_%d.fq", outFN, i + 1); |
| 79 | break; |
| 80 | } |
| 81 | |
| 82 | for (int i = 0; i < n_os; i++) |
| 83 | os[i] = new ofstream(outReadF[i]); |
| 84 | } |
| 85 | |
| 86 | template<class ReadType, class ModelType> |
| 87 | void simulate(char* modelF, char* resultsF) { |