0

多くのプロセスをデータファイルから読み取らせるC/MPIプログラムを書いています。stdioの標準関数(fopen、fread、fseek)を使用すると、すべてがうまくいきます。4Goオフセットを超えることができないという問題。そのため、MPI-IO関数を使用して大きなファイルを読み取りましたが、現時点ではメモリが十分に解放されていません。

実際、私はバッファを読み取り、それを処理してから、割り当てられたメモリを解放します。プロセスごとのメモリ使用量は完璧ですが、グローバルメモリ使用量の増加は止まりません。mpi_file_read atをfreadに置き換えるだけで、この問題は発生しません。

私のコードがあります:

double CPUtime(){ return ((double) clock())/CLOCKS_PER_SEC;}int main(int argc, char* argv []){

if(argc != 5) {
    printf("\t[Dictionary file] [Dictionary] [Input file] [Buffer size]\n");
    exit(0);
}

char* sInput    = malloc (sizeof(char)*maxLength);
char* sOutput   = malloc (sizeof(char)*maxLength);
char* compl     = malloc (sizeof(char)*maxLength);

char* sDictionaryFileName   =   argv[1];
char* sDictionaryName       =   argv[2];
char* filename              =   argv[3];
int Mbuffer                 =   atoi(argv[4]);



int maxBuffer = Mbuffer*1024*1024;
int over      = 10000;

int rank,numprocess;
long int offset;


char* buffer;
char* opbuffer;

double tstart=CPUtime();

MPI_Init( &argc, &argv );
MPI_Comm_rank( MPI_COMM_WORLD, &rank );

/* mpi version */
/* open the file*/
MPI_File fh;
int err;
err = MPI_File_open(MPI_COMM_WORLD, filename, MPI_MODE_RDONLY, MPI_INFO_NULL, &fh);
if (err != MPI_SUCCESS) {
    char errstr[MPI_MAX_ERROR_STRING];
    int  errlen;
    MPI_Error_string(err, errstr, &errlen);
    printf("Error at opening file %s (%s)\n",filename,errstr);
    MPI_Finalize();
    exit(1);
}

// get offsets and buffer size
MPI_Offset sfile;
MPI_File_get_size(fh,&sfile);
MPI_Status status;


/* C version */
/*FILE* fh;
long int sfile;
fh =fopen( filename,"rb");
if (fh==NULL) {
    printf("Error at opening file %s\n",filename);
    exit(1);
}
// get offsets and buffer size
fseek(fh, 0L, SEEK_END);
sfile = ftell(fh);
fseek(fh, 0L, SEEK_SET);*/


MPI_Comm_size( MPI_COMM_WORLD, &numprocess );

/* number of iterations */
long int data_size = (long int)(sfile/(numprocess));
int nbIter = data_size/maxBuffer;
if(nbIter<=1){
    nbIter = 1;
    maxBuffer = data_size;
}

/* offsets */
offset = data_size*(rank);
long int cursor = offset;
char* header;
if(rank==0){
    FILE* fh;
    fh =fopen( filename,"rb");
    if (fh==NULL) {
        printf("Error at opening file %s\n",filename);
        exit(1);
    }
    /* read the header and broadcast it */
    header = malloc(sizeof(char)*1000);
    fgets(header,1000,fh);
    fclose(fh);

    //broadcast header
    int sndHeader = strlen(header);
    //cursor+=sndHeader;
    int process_counter;
    for(process_counter=1;process_counter<numprocess;process_counter++){
        int ierr = MPI_Send(&sndHeader, 1, MPI_INT, process_counter, 42,MPI_COMM_WORLD);
        if (ierr != MPI_SUCCESS) {
            int errclass,resultlen;
            char err_buffer[MPI_MAX_ERROR_STRING];
            MPI_Error_class(ierr,&errclass);
            if (errclass== MPI_ERR_RANK) {
                fprintf(stderr,"Invalid rank used in MPI send call\n");
                MPI_Error_string(ierr,err_buffer,&resultlen);
                fprintf(stderr,err_buffer);
                MPI_Finalize();
            }
        }
        MPI_Send(header, sndHeader, MPI_CHAR, process_counter, 43, MPI_COMM_WORLD);
    }
}
else{
    /* receive the header */
    int sizeofHeader;
    MPI_Status s ;
    MPI_Recv(&sizeofHeader,1,MPI_INT,0,42,MPI_COMM_WORLD,&s);
    header = malloc (sizeof(char)*sizeofHeader+1);
    MPI_Recv(header,sizeofHeader,MPI_CHAR,0,43,MPI_COMM_WORLD,&s);
}


/* Synchronization barrier */
MPI_Barrier(MPI_COMM_WORLD);


int count;

opbuffer = malloc(sizeof(char)*maxBuffer);

/* C version */
//fseek(fh,cursor,SEEK_SET);

for(count=0;count<nbIter;count++){

    if(count==0 && rank==numprocess-1){ //init ring
        //send the token to p0
        int token=1;
        MPI_Send(&token,sizeof(int),MPI_INT,0,55,MPI_COMM_WORLD);
    }

    //recv
    int token;
    int sender;
    if(rank==0)
        sender = numprocess-1;
    else
        sender=rank-1;

    MPI_Status s;
    MPI_Recv(&token,sizeof(int),MPI_INT,sender,55,MPI_COMM_WORLD,&s);
    fflush(stdout);printf("P%d got the token at %G\n",rank,CPUtime());
    //read
    double start=CPUtime();
    /*double readtime;
    double sread=CPUtime();//read time*/

    //read
    if(token==1){
        /* MPI version */
        int err=MPI_File_read_at(fh, cursor,opbuffer,  sizeof(char)*maxBuffer, MPI_CHAR, &status);
        if(err!=MPI_SUCCESS){
            /*char errstr[MPI_MAX_ERROR_STRING];
            int  errlen;
            MPI_Error_string(err, errstr, &errlen);
            printf("Error reading file %s (%s)\n",filename,errstr);*/
            MPI_Finalize();
            exit(0);
        }

        /* C version of read */
        /*int k=fread(opbuffer,sizeof(char),maxBuffer,fh);
        if(k==0)
            perror("fread");*/

        cursor+=maxBuffer;
        buffer=opbuffer;

    }
    else{
        printf("Error token!\n");
        token=1;
    }
    //printf("P%d readtime=%G\n",rank,CPUtime()-sread);
    //Isend
    int next = (rank+1)%numprocess;
    MPI_Send(&token,sizeof(int),MPI_INT,next,55,MPI_COMM_WORLD);

    /* start processing*/ 
    /* end processing */


}
free(opbuffer);
int er=MPI_File_close(&fh);
if(er!=MPI_SUCCESS){
    printf("Error closing file\n");
    MPI_Finalize();
    exit(1);
}
MPI_Finalize();

printf("Global time : %G\n",CPUtime()-tstart);
return 0;
}

誰かがそれが何であるかについて何か考えを持っているならば、私はそれを認めるでしょう。ありがとうございました。

4

1 に答える 1

1

おそらく、あなたは一度も電話していないのでしょうMPI_File_close。これにより、ファイルに対する中間操作がリークします。if(err!=MPI_SUCCESS)本当にきれいなコードを書きたい場合は、エラー状態でも閉じる必要があることに注意してください。

于 2012-07-11T16:52:05.030 に答える