Xujun, Regarding your original question: please send the complete error message. Dmitry.
On Mon, Aug 17, 2015 at 11:15 AM Xujun Zhao <[email protected]> wrote: > Ahhhh, I should drink some coffee in the morning. > Now it passed the test! > > On Mon, Aug 17, 2015 at 11:13 AM, Matthew Knepley <[email protected]> > wrote: > >> On Mon, Aug 17, 2015 at 11:02 AM, Xujun Zhao <[email protected]> wrote: >> >>> No. It gives the following error msg: >>> >> >> Did you build the executable? >> >> cd src/vec/vec/examples/tutorials >> make ex43 >> >> Matt >> >> >>> mpirun -np 2 ex43 >>> >>> [proxy:0:[email protected]] [proxy:0:[email protected]] >>> HYDU_create_process (utils/launch/launch.c:75): HYDU_create_process >>> (utils/launch/launch.c:75): execvp error on file ex43 (No such file or >>> directory) >>> >>> execvp error on file ex43 (No such file or directory) >>> >>> On Mon, Aug 17, 2015 at 10:57 AM, Matthew Knepley <[email protected]> >>> wrote: >>> >>>> On Mon, Aug 17, 2015 at 10:49 AM, Xujun Zhao <[email protected]> wrote: >>>> >>>>> Hi all, >>>>> >>>>> I want PETSc to generate random vector using VecSetRandom() following >>>>> given examples, but failed and showed some "out of memory" error. The >>>>> following is the code, which goes well until it reaches VecSetRandom(). >>>>> Can >>>>> anyone help me figure out the reason? Thanks a lot. >>>>> >>>> >>>> Does src/vec/vec/examples/tests/ex43.c run for you? >>>> >>>> Thanks, >>>> >>>> Matt >>>> >>>> >>>>> XZ >>>>> >>>>> >>>>> >>>>> -------------------------------------------------------------------------------------------- >>>>> Vec u; >>>>> PetscRandom rand_ctx; /* random number generator context */ >>>>> PetscMPIInt size, rank; >>>>> PetscInt n, dn; >>>>> >>>>> >>>>> MPI_Comm_rank(PETSC_COMM_WORLD,&rank);//CHKERRQ(ierr); >>>>> MPI_Comm_size(PETSC_COMM_WORLD,&size);//CHKERRQ(ierr); >>>>> n = N/size + 1; >>>>> dn = n*size - N; >>>>> if ( dn>0 && rank<dn ) n -= 1; >>>>> printf("--->test in petsc_random_vector(): rank = %d, n = >>>>> %d\n",rank,n); >>>>> >>>>> >>>>> VecCreate(PETSC_COMM_WORLD,&u); >>>>> VecSetSizes(u,n,N); >>>>> PetscRandomCreate(PETSC_COMM_WORLD, &rand_ctx); >>>>> #if defined(PETSC_HAVE_DRAND48) >>>>> PetscRandomSetType(rand_ctx,PETSCRAND48); >>>>> #elif defined(PETSC_HAVE_RAND) >>>>> PetscRandomSetType(rand_ctx,PETSCRAND); >>>>> #endif >>>>> PetscRandomSetFromOptions(rand_ctx); >>>>> >>>>> >>>>> VecSetRandom(u,rand_ctx); >>>>> PetscRandomDestroy(&rand_ctx); >>>>> >>>> >>>> >>>> >>>> -- >>>> What most experimenters take for granted before they begin their >>>> experiments is infinitely more interesting than any results to which their >>>> experiments lead. >>>> -- Norbert Wiener >>>> >>> >>> >> >> >> -- >> What most experimenters take for granted before they begin their >> experiments is infinitely more interesting than any results to which their >> experiments lead. >> -- Norbert Wiener >> > >
