Re: [OMPI users] problems on parallel writing

2010-02-25 Thread w k
Hi Jody, I tried your suggestion but it still failed. Attached is the modified code. If your machine has fortran compiler as well, you can try it. BTW, how many processors did you use for testing your C code? Thanks, Kan On Thu, Feb 25, 2010 at 3:35 AM, jody wrote: >

Re: [OMPI users] problems on parallel writing

2010-02-25 Thread jody
Hi Just wanted to let you know: I translated your program to C ran it, and it crashed at MPI_FILE_SET_VIEW in a similar way than yours did. then i added an if-clause to prevent the call of MPI_FILE_WRITE with the undefined value. if (myid == 0) { MPI_File_write(fh, temp, count,

Re: [OMPI users] problems on parallel writing

2010-02-24 Thread Terry Frankcombe
On Wed, 2010-02-24 at 13:40 -0500, w k wrote: > Hi Jordy, > > I don't think this part caused the problem. For fortran, it doesn't > matter if the pointer is NULL as long as the count requested from the > processor is 0. Actually I tested the code and it passed this part > without problem. I

Re: [OMPI users] problems on parallel writing

2010-02-24 Thread jody
Hi I can't answer your question about the array q offhand, but i will try to translate your program to C and see if it fails the same way. Jody On Wed, Feb 24, 2010 at 7:40 PM, w k wrote: > Hi Jordy, > > I don't think this part caused the problem. For fortran, it doesn't

Re: [OMPI users] problems on parallel writing

2010-02-24 Thread w k
Hi Jordy, I don't think this part caused the problem. For fortran, it doesn't matter if the pointer is NULL as long as the count requested from the processor is 0. Actually I tested the code and it passed this part without problem. I believe it aborted at MPI_FILE_SET_VIEW part. Just curious,

Re: [OMPI users] problems on parallel writing

2010-02-24 Thread jody
Hi I know nearly nothing about fortran but it looks to me as the pointer 'temp' in > call MPI_FILE_WRITE(FH, temp, COUNT, MPI_REAL8, STATUS, IERR) is not defined (or perhaps NULL?) for all processors except processor 0 : > if ( myid == 0 ) then > count = 1 > else > count = 0 > end if

[OMPI users] problems on parallel writing

2010-02-23 Thread w k
Hello everyone, I'm trying to implement some functions in my code using parallel writing. Each processor has an array, say q, whose length is single_no(could be zero on some processors). I want to write q down to a common file, but the elements of q would be scattered to their locations in this