Question

我编写了一个C代码，它从大小约为1 GB的二进制文件中提取数据。有101（0到100）个配置，C代码提取所选配置的数据，并将输出写入文件。要编译C代码，我在终端中给出用户定义的配置号：

gcc binary2textperconfig.c -o f.out
./f.out proton-p000-1.bin out1.txt

然后它要求配置号：

Enter the configuration number:

之后，数据写入文件“out0.txt”。现在我想为所有101个配置运行此代码并将这些数据写入文件“out0.txt”，“out1.txt”，....，“out100.txt”等。我不知道该怎么做没有在终端中输入101次配置号码。有人可以帮帮我吗？这是我的C代码：

#include <stdio.h>
#include<complex.h>
#include<math.h>
#include <stdlib.h>
#include <gsl/gsl_sf_gamma.h>
#include <gsl/gsl_matrix.h>


typedef double complex  dcomplex;


//Data is converted to Bigendian using io-general

double constructfloat(char bytes[sizeof(double)/2],  int order)
{
    double dRes;
    char *pc;
    int i;

    if (order == 0)
        for(i=0, pc = (char*) &dRes; i<=sizeof(double)-1 ; i++, pc++)
            (*pc) = bytes[i];
    else
        for(i=sizeof(double)-1, pc = (char*) &dRes; i>=0; i--, pc++)
            (*pc) = bytes[i];

    return dRes;
}


int main(int argc, char *argv[]){



    int configcount = 101;
    int mcount = 14;
    int tcount = 64;
    int d1count = 4;
    int d2count = 4;
    int pcount = 45;


    int userci;
    int usermi;
    int userpi;


    // number of complex numbers per configuration
    int unitcountperconfig =(mcount*tcount*d1count*d2count*pcount);


    // initialize loop index variables
    int ci = 0; //config
    int mi = 0; //mass
    int ti = 0;
    int d1i = 0;
    int d2i = 0;
    int pi = 0; //momentum

    // for holding the result of read operation ( how many units have been read)
    int result;

    // for holding the data read from file
    char * cbuff;


    // input file handle from where binary data is read
    FILE * fin = fopen(argv[1],"rb");

    // if the input file cannot be read for reading, close opened file handles, show an error message to the user, and exit
    if (fin==NULL)
    {
        fputs ("Error opening input file\n",stderr);
        exit (1);
    }

    FILE * fout = fopen(argv[2],"wt");

    // if the output file cannot be opened for writing, close opened file handles, show an error message to the user, and exit
    if (fout==NULL)
    {
        fclose(fin);
        fputs ("Error opening output file\n",stderr);
        exit (1);
    }


    // take input from the user
    // take input from the user
    printf("Enter the configuration number: ");
    scanf("%d",&userci);

    // allocate memory to contain the chunk of data for a time slice:
    cbuff = (char*)malloc(sizeof(dcomplex)*unitcountperconfig );

    // show error message and exit if memory allocation failed
    if(cbuff == NULL)
    {
        fputs("Buffer allocation failed.", stderr);
        exit(1);
    }

    // variable to hold a complex number read from the file
    dcomplex aComplexNumber;

    dcomplex sumpertimeslice[tcount];


    // loop on time slices
    for( ci = 0;  ci< configcount ; ci++){

        // index of the complex number being read
        unsigned int cNumberIdx = 0;

        // debugging message
        printf("reading data for configuration: %d\n",ci);

        // perform read operation to read the desired chunk of data
        result = fread(cbuff,  sizeof(char),  sizeof(dcomplex)*unitcountperconfig, fin );


        // if size of data successfully read is not equal to what we wanted to read, notify the user and exit
        if (result != sizeof(dcomplex)*unitcountperconfig) {
            fputs ("data reading error\n",stderr);
            exit (3);
        }

        double realP;
        double imagP;// variable to hold real and imaginary part of the complex number


        double realPSum;
        double imagPSum;// variable to hold sum of real and sum of imaginary part of the current sum per time slice


            for (mi =0; mi< mcount ; mi++){

                for (ti =0; ti< tcount ; ti++){



                // array to hold trace for each time slice
                sumpertimeslice[ti] = 0.0 + 0.0*_Complex_I;

                    for (d1i =0; d1i < d1count ; d1i++){
                        for (d2i =0; d2i < d2count ; d2i++){
                            for (pi =0; pi < pcount ; pi++){
                            aComplexNumber = constructfloat( &cbuff[cNumberIdx], 0 ) + constructfloat( &cbuff[cNumberIdx+ ((int)sizeof(dcomplex))/2 ], 0 )*_Complex_I;
                if (ci== userci)
                {
            cNumberIdx += (int)sizeof(dcomplex);
                            if (cimag(aComplexNumber)>0)
                            {fprintf( fout, "%d,%d,%d,%d,%d,%d,%e+%ei\n" ,ci+1, mi+1,ti+1, d1i+1,d2i+1,pi+1,creal( aComplexNumber ),cimag( aComplexNumber ) );}
                            else
                            {fprintf( fout, "%d,%d,%d,%d,%d,%d,%e%ei\n" ,ci+1, mi+1,ti+1, d1i+1,d2i+1,pi+1,creal(  aComplexNumber ),cimag( aComplexNumber ) );}
                }


                            }               
                        }           
                    }


                }           
            }

            }



    // free the allocated memory    
    free(cbuff);

    // close the opened file handles
    fclose(fin);
    fclose(fout);
    //fclose(complexNumberFileP);

}

Answer 1

使用seq实用程序生成0到0之间的数字列表。 100和send it as a string to stdin。

for CNUMBER in $(seq 0 100); do
    ./f.out proton-p000-1.bin out${CNUMBER}.txt <<< "${CNUMBER}"
done

或

for CNUMBER in $(seq 0 100); do
    echo $CNUMBER | ./f.out proton-p000-1.bin out${CNUMBER}.txt
done

循环使用C代码输入并立即使用终端将数据写入不同的文件

1 个答案: