[petsc-users] Petsc creates a random vector
Xujun Zhao
xzhao99 at gmail.com
Mon Aug 17 11:02:24 CDT 2015
No. It gives the following error msg:
mpirun -np 2 ex43
[proxy:0:0 at mcswl164.mcs.anl.gov] [proxy:0:0 at mcswl164.mcs.anl.gov]
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 <knepley at gmail.com> wrote:
> On Mon, Aug 17, 2015 at 10:49 AM, Xujun Zhao <xzhao99 at gmail.com> 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150817/0aa68bd9/attachment-0001.html>
More information about the petsc-users
mailing list