int main(int argc, char *const argv[])
{
int ret = EXIT_SUCCESS;
unsigned int err;
unsigned int n;
double *ary;
if (argc < 2)
usage();
err = read_number_of_values(stdin, &n);
if (err)
return EXIT_SUCCESS;
ary = malloc(n * sizeof(*ary));
if (!ary) {
perror("malloc failed");
goto out;
}
err = read_doubles(stdin, ary, n);
if (err)
goto out;
err = do_sum(stdout, argv[1], ary, n);
if (err)
goto out;
out:
if (ary)
free(ary);
return ret;
}
int main(void) {
int i, n;
double arr[MAX];
n = read_doubles(stdin, arr, MAX);
for (i = 0; i < n; i++)
printf("%3d %g\n", i, arr[i]);
return 0;
}
//read a whole ape propagators
void read_ape_propagators(ccss_propagator **out,const char *path,int nmass)
{
FILE *fin=open_file(path,"r");
ape_propagator buf[2][nmass];
read_doubles((double*)buf,fin,2*nmass*sizeof(ape_propagator)/sizeof(double));
for(int imass=0;imass<nmass;imass++)
for(int r=0;r<2;r++)
for(int imom=0;imom<nmom;imom++)
for(int ic_so=0;ic_so<3;ic_so++)
for(int id_so=0;id_so<4;id_so++)
for(int ic_si=0;ic_si<3;ic_si++)
for(int id_si=0;id_si<4;id_si++)
{
memcpy(out[r][imass][imom][ic_si][ic_so][id_si][id_so],buf[r][imass][imom][ic_so][id_so][ic_si][id_si],sizeof(complex));
printf("%d %d %d %d %d %d %d %lg %lg %lg %lg\n",imass,r,imom,ic_so,id_so,ic_si,id_si,out[r][imass][imom][ic_si][ic_so][id_si][id_so][0],out[r][imass][imom][ic_si][ic_so][id_si][id_so][1],buf[r][imass][imom][ic_so][id_so][ic_si][id_si][0],buf[r][imass][imom][ic_so][id_so][ic_si][id_si][1]);
}
}
