7 void main(
int argc,
char *argv[])
19 int foldvals[5] = {1,2,4,8,16};
35 float **power = malloc(nfolds *
sizeof(*power));
39 fp = fopen(argv[1],
"r");
41 while ( fgets(line, 256, fp) != NULL ) nf++;
43 fprintf(stderr,
"read %d lines of data\n",nf);
46 for(fold = 0; fold < nfolds; fold++){
47 power[fold] = malloc(nf *
sizeof(**power));
48 assert(power[fold] != NULL);
52 sp = (
float *) malloc(nf*4);
56 fseek(fp, 0, SEEK_SET);
60 while ( fgets(line, 250, fp) != NULL ) {
61 sp[i] = strtof(line, NULL);
64 fprintf(stderr,
"converted %d lines of data\n",i);
68 for (fold=0; fold< nfolds; fold++){
69 fval = foldvals[fold];
71 if (fval == 1)
for (n=0;n<nf;n++) power[0][n] = sp[n];
73 NFS = fval*nf1 - fval/2;
75 for (n=NFS; n< nf; n++){
78 for (m=0;m<foldvals[fold];m++){
79 x = n * m * xdiv + 0.5;
81 if ( ii > 1 && ii != lstidx) power[fold][n] = power[fold][n] + sp[ii];
89 for (ii=0;ii<20000;ii++){
90 fprintf(stdout,
"%07d %10.5f %10.5f %10.5f %10.5f %10.5f %10.5f",
91 ii, sp[ii],power[0][ii],power[1][ii],power[2][ii],power[3][ii],power[4][ii]);
96 for(fold = 0; fold < nfolds; fold++) free(power[fold]);