Actual source code: veccreate.c


  2: #include <petsc/private/vecimpl.h>

  4: /*@
  5:   VecCreate - Creates an empty vector object. The type can then be set with VecSetType(),
  6:   or VecSetFromOptions().

  8:    If you never  call VecSetType() or VecSetFromOptions() it will generate an
  9:    error when you try to use the vector.

 11:   Collective

 13:   Input Parameter:
 14: . comm - The communicator for the vector object

 16:   Output Parameter:
 17: . vec  - The vector object

 19:   Level: beginner

 21: .seealso: `VecSetType()`, `VecSetSizes()`, `VecCreateMPIWithArray()`, `VecCreateMPI()`, `VecDuplicate()`,
 22:           `VecDuplicateVecs()`, `VecCreateGhost()`, `VecCreateSeq()`, `VecPlaceArray()`
 23: @*/
 24: PetscErrorCode VecCreate(MPI_Comm comm, Vec *vec)
 25: {
 26:   Vec v;

 29:   *vec = NULL;
 30:   VecInitializePackage();

 32:   PetscHeaderCreate(v, VEC_CLASSID, "Vec", "Vector", "Vec", comm, VecDestroy, VecView);

 34:   PetscLayoutCreate(comm, &v->map);
 35:   v->array_gotten = PETSC_FALSE;
 36:   v->petscnative  = PETSC_FALSE;
 37:   v->offloadmask  = PETSC_OFFLOAD_UNALLOCATED;
 38: #if defined(PETSC_HAVE_VIENNACL) || defined(PETSC_HAVE_CUDA) || defined(PETSC_HAVE_HIP)
 39:   v->minimum_bytes_pinned_memory = 0;
 40:   v->pinned_memory               = PETSC_FALSE;
 41: #endif
 42: #if defined(PETSC_HAVE_DEVICE)
 43:   v->boundtocpu = PETSC_TRUE;
 44: #endif
 45:   PetscStrallocpy(PETSCRANDER48, &v->defaultrandtype);
 46:   *vec = v;
 47:   return 0;
 48: }