#include<cstdio>
#include<cstring>
#include<string>
+#include<cassert>
#include "utils.h"
#include "RefSeq.h"
logp = 0.0;
sum = 0.0;
for (int i = 0; i < NCODES; i++) sum += (p[i] + c[i]);
- if (sum <= 0.0) return;
+ if (sum <= EPSILON) return;
for (int i = 0; i < NCODES; i++) {
p[i] = (p[i] + c[i]) / sum;
if (c[i] > 0.0) { logp += c[i] * log(p[i]); }
int tmp_ncodes;
memset(c, 0, sizeof(c));
- fscanf(fi, "%d", &tmp_ncodes);
+ assert(fscanf(fi, "%d", &tmp_ncodes) == 1);
assert(tmp_ncodes == NCODES);
for (int i = 0; i < NCODES; i++)
- fscanf(fi, "%lf", &p[i]);
+ assert(fscanf(fi, "%lf", &p[i]) == 1);
}
void NoiseProfile::write(FILE *fo) {