petsc_utils_01_createVec.cc

Go to the documentation of this file.

#include "petsc_utils.h"
 
#include "chi_log.h"
 
//###################################################################
/**Creates a general vector.
 *
This is a macro for:
\code
Vec x;
VecCreate(PETSC_COMM_WORLD,&x);
VecSetType(x,VECMPI);
VecSetSizes(x, local_size, global_size);
VecSetOption(x,VEC_IGNORE_NEGATIVE_INDICES,PETSC_TRUE);
 
return x;
\endcode*/
Vec chi_math::PETScUtils::
CreateVector(int64_t local_size, int64_t global_size)
{
  Vec x;
  VecCreate(PETSC_COMM_WORLD,&x);
  VecSetType(x,VECMPI);
  VecSetSizes(x, local_size, global_size);
  VecSetOption(x,VEC_IGNORE_NEGATIVE_INDICES,PETSC_TRUE);
 
  return x;
}
 
//###################################################################
/**Creates a general vector.
 *
This is a macro for:
\code
VecCreate(PETSC_COMM_WORLD,&x);
VecSetType(x,VECMPI);
VecSetSizes(x, local_size, global_size);
VecSetOption(x,VEC_IGNORE_NEGATIVE_INDICES,PETSC_TRUE);
\endcode*/
void chi_math::PETScUtils::
CreateVector(Vec& x, int64_t local_size, int64_t global_size)
{
  VecCreate(PETSC_COMM_WORLD,&x);
  VecSetType(x,VECMPI);
  VecSetSizes(x, local_size, global_size);
  VecSetOption(x,VEC_IGNORE_NEGATIVE_INDICES,PETSC_TRUE);
}
 
//###################################################################
/**Creates a general vector with ghost value support.
 *
This is a macro for:
\code
Vec x;
VecCreateGhost(PETSC_COMM_WORLD,
               local_size,
               global_size,
               nghosts,
               ghost_indices.data(),
               &x);
 
VecSetOption(x,VEC_IGNORE_NEGATIVE_INDICES,PETSC_TRUE);
 
return x;
\endcode*/
Vec chi_math::PETScUtils::
CreateVectorWithGhosts(int64_t local_size, int64_t global_size,
                       int64_t nghosts,
                       const std::vector<int64_t>& ghost_indices)
{
  Vec x;
  VecCreateGhost(PETSC_COMM_WORLD,
                 local_size,
                 global_size,
                 nghosts,
                 (ghost_indices.empty())? NULL : ghost_indices.data(),
                 &x);
 
  VecSetOption(x,VEC_IGNORE_NEGATIVE_INDICES,PETSC_TRUE);
 
  return x;
}