void initialize() { p4423[0] = 1; for (index_t i = 1; i < max_len; i++) { p4423[i] = p4423[i - 1] * 4423; } gamma_init(); }
void gamma_save(int record, int cpu) { char name[128]; sprintf(name, "binData/gamma_%06d_%03d.txt", record, cpu); FILE *fp = fopen(name, "wt"); MF_ENSURE(fp, "Cannot open file '%s' for writing.", name); fprintf(fp, "%lg %d\n", gamma_step, gamma_multiplier); for (int i = 0; i < 100; i++) fprintf(fp, "%lg\n", gamma[i]); fclose(fp); gamma_init(); }