int cnt = 0;
while (fscanf(fi, "%d", &cvec[0]) == 1) {
for (int j = 1; j <= M; j++) assert(fscanf(fi, "%d", &cvec[j]) == 1);
int cnt = 0;
while (fscanf(fi, "%d", &cvec[0]) == 1) {
for (int j = 1; j <= M; j++) assert(fscanf(fi, "%d", &cvec[j]) == 1);
- gammas[j] = new gamma_dist(cvec[j]);
- rgs[j] = new gamma_generator(*(params->engine), *gammas[j]);
+ gammas[j] = NULL; rgs[j] = NULL;
+ if (cvec[j] > 0) {
+ gammas[j] = new gamma_dist(cvec[j]);
+ rgs[j] = new gamma_generator(*(params->engine), *gammas[j]);
+ }
}
if (verbose && cnt % 100 == 0) { printf("Thread %d, %d count vectors are processed!\n", params->no, cnt); }
}
if (verbose && cnt % 100 == 0) { printf("Thread %d, %d count vectors are processed!\n", params->no, cnt); }