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 >
