- //calculate normalized read fraction
- double *nrf = new double[M + 1];
- memset(nrf, 0, sizeof(double) * (M + 1));
- denom = 0.0;
- for (int i = 1; i <= M; i++)
- if (eel[i] > EPSILON) {
- nrf[i] = counts[i];
- denom += nrf[i];
- }
- else {
- if (counts[i] > EPSILON) { printf("Warning: An isoform which EEL is less than %.6g gets sampled!\n", MINEEL); }
- }
- assert(denom > 0.0);
- for (int i = 1; i <= M; i++) nrf[i] /= denom;
-