MathLib::PETScVector Class Reference

Detailed Description

Wrapper class for PETSc vector.

Definition at line 32 of file PETScVector.h.

#include <PETScVector.h>

Public Types
using	IndexType = PetscInt
using	PETSc_Vec = Vec

Public Member Functions
	PETScVector ()
	PETScVector (const PetscInt vec_size, const bool is_global_size=true)
	Constructor. More...
	PETScVector (const PetscInt vec_size, const std::vector< PetscInt > &ghost_ids, const bool is_global_size=true)
	Constructor. More...
	PETScVector (const PETScVector &existing_vec, const bool deep_copy=true)
	Copy constructor. More...
	PETScVector (PETScVector &&other)
	~PETScVector ()
void	finalizeAssembly ()
	Perform MPI collection of assembled entries in buffer. More...
PetscInt	size () const
	Get the global size of the vector. More...
PetscInt	getLocalSize () const
	Get the number of entries in the same rank. More...
PetscInt	getGhostSize () const
	Get the number of ghost entries in the same rank. More...
PetscInt	getRangeBegin () const
	Get the start index of the local vector. More...
PetscInt	getRangeEnd () const
	Get the end index of the local vector. More...
void	set (const PetscInt i, const PetscScalar value)
void	add (const PetscInt i, const PetscScalar value)
template<class T_SUBVEC >
void	add (const std::vector< PetscInt > &e_idxs, const T_SUBVEC &sub_vec)
template<class T_SUBVEC >
void	set (const std::vector< PetscInt > &e_idxs, const T_SUBVEC &sub_vec)
void	setZero ()
PETScVector &	operator= (PETScVector &&)=delete
void	setLocalAccessibleVector () const
std::vector< PetscScalar >	get (std::vector< IndexType > const &indices) const
PetscScalar	operator[] (PetscInt idx) const
void	getGlobalVector (std::vector< PetscScalar > &u) const
PetscScalar	get (const PetscInt idx) const
PETSc_Vec &	getRawVector ()
	Exposes the underlying PETSc vector. More...
PETSc_Vec const &	getRawVector () const
void	copyValues (std::vector< PetscScalar > &u) const
void	viewer (const std::string &file_name, const PetscViewerFormat vw_format=PETSC_VIEWER_ASCII_MATLAB) const
void	shallowCopy (const PETScVector &v)

Private Member Functions
void	destroy ()
void	gatherLocalVectors (PetscScalar local_array[], PetscScalar global_array[]) const
	Collect local vectors. More...
PetscScalar *	getLocalVector () const
PetscInt	getLocalIndex (const PetscInt global_index) const
	Get local index by a global index. More...
void	restoreArray (PetscScalar *array) const
void	config ()
	A function called by constructors to configure members. More...

Private Attributes
PETSc_Vec	v_ = nullptr
PETSc_Vec	v_loc_ = nullptr
PetscInt	start_rank_
	Starting index in a rank. More...
PetscInt	end_rank_
	Ending index in a rank. More...
PetscInt	size_
	Size of the vector. More...
PetscInt	size_loc_
	Size of local entries. More...
PetscInt	size_ghosts_ = 0
	Size of local ghost entries. More...
bool	has_ghost_id_ = false
	Flag to indicate whether the vector is created with ghost entry indices. More...
std::vector< PetscScalar >	entry_array_
	Array containing the entries of the vector. If the vector is created without given ghost IDs, the array contains all entries of the global vector, v_. Otherwise it only contains the entries of the global vector owned by the current rank. More...
std::map< PetscInt, PetscInt >	global_ids2local_ids_ghost_
	Map global indices of ghost entries to local indices. More...

Friends
void	finalizeVectorAssembly (PETScVector &vec)
	Function to finalize the vector assembly. More...

Member Typedef Documentation

IndexType

using MathLib::PETScVector::IndexType = PetscInt

Definition at line 35 of file PETScVector.h.

PETSc_Vec

using MathLib::PETScVector::PETSc_Vec = Vec

Definition at line 38 of file PETScVector.h.

Constructor & Destructor Documentation

PETScVector() [1/5]

MathLib::PETScVector::PETScVector ( )

inline

Definition at line 41 of file PETScVector.h.

{}

PETScVector() [2/5]

MathLib::PETScVector::PETScVector	(	const PetscInt	vec_size,
		const bool	is_global_size = true
	)

Constructor.

Parameters

vec_size	The size of the vector, either global or local
is_global_size	The flag of the type of vec_size, i.e. whether it is a global size or local size. The default is true. If is_global_size is true, the vector is created by the global size, the local size of the vector is determined by PETSc, and vice versa is the same.

Definition at line 30 of file PETScVector.cpp.

{
    if (is_global_size)
    {
        VecCreate(PETSC_COMM_WORLD, &v_);
        VecSetSizes(v_, PETSC_DECIDE, vec_size);
    }
    else
    {
        // Fix size partitioning
        // the size can be associated to specific memory allocation of a matrix
        VecCreateMPI(PETSC_COMM_WORLD, vec_size, PETSC_DECIDE, &v_);
    }
 
    config();
}

References config(), and v_.

PETScVector() [3/5]

MathLib::PETScVector::PETScVector	(	const PetscInt	vec_size,
		const std::vector< PetscInt > &	ghost_ids,
		const bool	is_global_size = true
	)

Constructor.

Parameters

vec_size	The size of the vector, either global or local
ghost_ids	The global indices of ghost entries
is_global_size	The flag of the type of vec_size, i.e. whether it is a global size or local size. The default is true.

Definition at line 47 of file PETScVector.cpp.

    : size_ghosts_{static_cast<PetscInt>(ghost_ids.size())}, has_ghost_id_{true}
{
    PetscInt nghosts = static_cast<PetscInt>(ghost_ids.size());
    if (is_global_size)
    {
        VecCreateGhost(PETSC_COMM_WORLD, PETSC_DECIDE, vec_size, nghosts,
                       ghost_ids.data(), &v_);
    }
    else
    {
        VecCreate(PETSC_COMM_WORLD, &v_);
        VecSetType(v_, VECMPI);
        VecSetSizes(v_, vec_size, PETSC_DECIDE);
        VecMPISetGhost(v_, nghosts, ghost_ids.data());
    }
 
    config();
 
    for (PetscInt i = 0; i < nghosts; i++)
        global_ids2local_ids_ghost_.emplace(ghost_ids[i], size_loc_ + i);
}

References config(), global_ids2local_ids_ghost_, size_loc_, and v_.

PETScVector() [4/5]

MathLib::PETScVector::PETScVector ( const PETScVector &  existing_vec, const bool  deep_copy = true  )

explicit

Copy constructor.

Parameters

existing_vec	The vector to be copied
deep_copy	The flag for a deep copy, which means to copy the values as well, the default is true

Definition at line 72 of file PETScVector.cpp.

{
    shallowCopy(existing_vec);
 
    // Copy values
    if (deep_copy)
    {
        VecCopy(existing_vec.v_, v_);
    }
}

References shallowCopy(), and v_.

PETScVector() [5/5]

MathLib::PETScVector::PETScVector

(

PETScVector &&

other

)

Definition at line 83 of file PETScVector.cpp.

    : v_{std::move(other.v_)},
      v_loc_{std::move(other.v_loc_)},
      start_rank_{other.start_rank_},
      end_rank_{other.end_rank_},
      size_{other.size_},
      size_loc_{other.size_loc_},
      size_ghosts_{other.size_ghosts_},
      has_ghost_id_{other.has_ghost_id_},
      global_ids2local_ids_ghost_{other.global_ids2local_ids_ghost_}
{
}

~PETScVector()

MathLib::PETScVector::~PETScVector ( )

inline

Definition at line 77 of file PETScVector.h.

{ destroy(); }

References destroy().

Member Function Documentation

add() [1/2]

void MathLib::PETScVector::add ( const PetscInt  i, const PetscScalar  value  )

inline

Add a value to an entry.

Parameters

i	Number of the entry
value	Value.

Definition at line 107 of file PETScVector.h.

    {
        VecSetValue(v_, i, value, ADD_VALUES);
    }

References v_.

add() [2/2]

template<class T_SUBVEC >

void MathLib::PETScVector::add ( const std::vector< PetscInt > &  e_idxs, const T_SUBVEC &  sub_vec  )

inline

Add values to several entries.

Parameters

e_idxs	Indices of entries to be added Note: std::size_t cannot be the type of e_idxs template argument.
sub_vec	Entries to be added.

Definition at line 120 of file PETScVector.h.

    {
        VecSetValues(v_, e_idxs.size(), &e_idxs[0], &sub_vec[0], ADD_VALUES);
    }

References v_.

config()

void MathLib::PETScVector::config ( )

private

A function called by constructors to configure members.

Definition at line 96 of file PETScVector.cpp.

{
    VecSetFromOptions(v_);
    VecGetOwnershipRange(v_, &start_rank_, &end_rank_);
 
    VecGetLocalSize(v_, &size_loc_);
    VecGetSize(v_, &size_);
 
    VecSetOption(v_, VEC_IGNORE_NEGATIVE_INDICES, PETSC_TRUE);
}

References end_rank_, size_, size_loc_, start_rank_, and v_.

Referenced by PETScVector().

copyValues()

void MathLib::PETScVector::copyValues

(

std::vector< PetscScalar > &

u

)

const

Copy local entries including ghost ones to an array

Parameters

u	Preallocated vector for the values of local entries.

Definition at line 197 of file PETScVector.cpp.

{
    assert(u.size() == (std::size_t)(getLocalSize() + getGhostSize()));
 
    PetscScalar* loc_x = getLocalVector();
    std::copy_n(loc_x, getLocalSize() + getGhostSize(), u.begin());
    restoreArray(loc_x);
}

References getGhostSize(), getLocalSize(), getLocalVector(), and restoreArray().

destroy()

void MathLib::PETScVector::destroy ( )

inlineprivate

Definition at line 221 of file PETScVector.h.

    {
        if (v_)
            VecDestroy(&v_);
        v_ = nullptr;
    }

References v_.

Referenced by ~PETScVector(), and shallowCopy().

finalizeAssembly()

void MathLib::PETScVector::finalizeAssembly

(

)

Perform MPI collection of assembled entries in buffer.

Definition at line 107 of file PETScVector.cpp.

{
    VecAssemblyBegin(v_);
    VecAssemblyEnd(v_);
}

References v_.

Referenced by MathLib::applyKnownSolution(), and MathLib::LinAlg::finalizeAssembly().

gatherLocalVectors()

void MathLib::PETScVector::gatherLocalVectors ( PetscScalar  local_array[], PetscScalar  global_array[]  ) const

private

Collect local vectors.

Parameters

local_array	Local array
global_array	Global array

Definition at line 113 of file PETScVector.cpp.

{
    // Collect vectors from processors.
    int size_rank;
    MPI_Comm_size(PETSC_COMM_WORLD, &size_rank);
 
    // number of elements to be sent for each rank
    std::vector<PetscInt> i_cnt(size_rank);
 
    MPI_Allgather(&size_loc_, 1, MPI_INT, &i_cnt[0], 1, MPI_INT,
                  PETSC_COMM_WORLD);
 
    // collect local array
    PetscInt offset = 0;
    // offset in the receive vector of the data from each rank
    std::vector<PetscInt> i_disp(size_rank);
    for (PetscInt i = 0; i < size_rank; i++)
    {
        i_disp[i] = offset;
        offset += i_cnt[i];
    }
 
    MPI_Allgatherv(local_array, size_loc_, MPI_DOUBLE, global_array, &i_cnt[0],
                   &i_disp[0], MPI_DOUBLE, PETSC_COMM_WORLD);
}

References size_loc_.

Referenced by getGlobalVector().

get() [1/2]

PetscScalar MathLib::PETScVector::get

(

const PetscInt

idx

)

const

Definition at line 206 of file PETScVector.cpp.

{
    if (!global_ids2local_ids_ghost_.empty())
    {
        return entry_array_[getLocalIndex(idx)];
    }
 
    // Ghost entries, and its original index is 0.
    const PetscInt id_p = (idx == -size_) ? 0 : std::abs(idx);
    return entry_array_[id_p];
}

References entry_array_, getLocalIndex(), global_ids2local_ids_ghost_, and size_.

get() [2/2]

std::vector< PetscScalar > MathLib::PETScVector::get

(

std::vector< IndexType > const &

indices

)

const

Get several entries. setLocalAccessibleVector() must be called beforehand.

Definition at line 218 of file PETScVector.cpp.

{
    std::vector<PetscScalar> local_x(indices.size());
    // If VecGetValues can get values from different processors,
    // use VecGetValues(v_, indices.size(), indices.data(),
    //                    local_x.data());
 
    if (!global_ids2local_ids_ghost_.empty())
    {
        for (std::size_t i = 0; i < indices.size(); i++)
        {
            local_x[i] = entry_array_[getLocalIndex(indices[i])];
        }
    }
    else
    {
        for (std::size_t i = 0; i < indices.size(); i++)
        {
            // Ghost entries, and its original index is 0.
            const IndexType id_p =
                (indices[i] == -size_) ? 0 : std::abs(indices[i]);
            local_x[i] = entry_array_[id_p];
        }
    }
 
    return local_x;
}

References entry_array_, getLocalIndex(), global_ids2local_ids_ghost_, and size_.

Referenced by operator[]().

getGhostSize()

PetscInt MathLib::PETScVector::getGhostSize ( ) const

inline

Get the number of ghost entries in the same rank.

Definition at line 86 of file PETScVector.h.

{ return size_ghosts_; }

References size_ghosts_.

Referenced by copyValues().

getGlobalVector()

void MathLib::PETScVector::getGlobalVector

(

std::vector< PetscScalar > &

u

)

const

Get global vector

Parameters

u	Array to store the global vector. Memory allocation is needed in advance

Definition at line 140 of file PETScVector.cpp.

{
#ifdef TEST_MEM_PETSC
    PetscLogDouble mem1, mem2;
    PetscMemoryGetCurrentUsage(&mem1);
#endif
 
    assert(static_cast<PetscInt>(u.size()) == size_);
 
    PetscInt state;
    VecLockGet(v_, &state);
    if (state != 0)
    {
        OGS_FATAL("PETSc vector is already locked for {:s} access.",
                  state > 0 ? "read" : "write");
    }
    PetscScalar* xp = nullptr;
    VecGetArray(v_, &xp);
 
    gatherLocalVectors(xp, u.data());
 
    // This following line may be needed late on
    //  for a communication load balance:
    // MPI_Barrier(PETSC_COMM_WORLD);
 
    VecRestoreArray(v_, &xp);
 
// TEST
#ifdef TEST_MEM_PETSC
    PetscMemoryGetCurrentUsage(&mem2);
    PetscPrintf(
        PETSC_COMM_WORLD,
        "### Memory usage by Updating. Before :%f After:%f Increase:%d\n", mem1,
        mem2, (int)(mem2 - mem1));
#endif
}

References gatherLocalVectors(), OGS_FATAL, size_, and v_.

Referenced by setLocalAccessibleVector().

getLocalIndex()

PetscInt MathLib::PETScVector::getLocalIndex ( const PetscInt  global_index ) const

private

Get local index by a global index.

Definition at line 262 of file PETScVector.cpp.

{
    if (global_index >= 0)  // non-ghost entry.
        return global_index - start_rank_;
 
    // A special case for a ghost location with global index equal to
    // the size of the local vector:
    PetscInt real_global_index = (-global_index == size_) ? 0 : -global_index;
 
    return global_ids2local_ids_ghost_.at(real_global_index);
}

References global_ids2local_ids_ghost_, size_, and start_rank_.

Referenced by get().

getLocalSize()

PetscInt MathLib::PETScVector::getLocalSize ( ) const

inline

Get the number of entries in the same rank.

Definition at line 84 of file PETScVector.h.

{ return size_loc_; }

References size_loc_.

Referenced by copyValues().

getLocalVector()

PetscScalar * MathLib::PETScVector::getLocalVector ( ) const

private

Get local vector, i.e. entries in the same rank

Definition at line 247 of file PETScVector.cpp.

{
    PetscScalar* loc_array;
    if (has_ghost_id_)
    {
        VecGhostUpdateBegin(v_, INSERT_VALUES, SCATTER_FORWARD);
        VecGhostUpdateEnd(v_, INSERT_VALUES, SCATTER_FORWARD);
        VecGhostGetLocalForm(v_, &v_loc_);
        VecGetArray(v_loc_, &loc_array);
    }
    else
        VecGetArray(v_, &loc_array);
    return loc_array;
}

References has_ghost_id_, v_, and v_loc_.

Referenced by copyValues(), and setLocalAccessibleVector().

getRangeBegin()

PetscInt MathLib::PETScVector::getRangeBegin ( ) const

inline

Get the start index of the local vector.

Definition at line 88 of file PETScVector.h.

{ return start_rank_; }

References start_rank_.

getRangeEnd()

PetscInt MathLib::PETScVector::getRangeEnd ( ) const

inline

Get the end index of the local vector.

Definition at line 90 of file PETScVector.h.

{ return end_rank_; }

References end_rank_.

getRawVector() [1/2]

PETSc_Vec& MathLib::PETScVector::getRawVector ( )

inline

Exposes the underlying PETSc vector.

Definition at line 171 of file PETScVector.h.

{ return v_; }

References v_.

Referenced by MathLib::LinAlg::axpby(), MathLib::LinAlg::axpy(), MathLib::LinAlg::aypx(), MathLib::LinAlg::componentwiseDivide(), MathLib::LinAlg::copy(), MathLib::LinAlg::matMult(), MathLib::LinAlg::matMultAdd(), MathLib::LinAlg::norm1(), MathLib::LinAlg::norm2(), MathLib::LinAlg::normMax(), MathLib::LinAlg::scale(), MathLib::LinAlg::set(), shallowCopy(), and MathLib::PETScLinearSolver::solve().

getRawVector() [2/2]

PETSc_Vec const& MathLib::PETScVector::getRawVector ( ) const

inline

Exposes the underlying PETSc vector.

Warning

This method is dangerous insofar as you can do arbitrary things also with a const PETSc vector!

Definition at line 178 of file PETScVector.h.

{ return v_; }

References v_.

operator=()

PETScVector& MathLib::PETScVector::operator= ( PETScVector &&  )

delete

Disallow moving.

Todo:

This operator should be implemented properly when doing a general cleanup of all matrix and vector classes.

operator[]()

PetscScalar MathLib::PETScVector::operator[] ( PetscInt  idx ) const

inline

Get the value of an entry by [] operator. setLocalAccessibleVector() must be called beforehand.

Definition at line 155 of file PETScVector.h.

{ return get(idx); }

References get().

restoreArray()

void MathLib::PETScVector::restoreArray ( PetscScalar *  array ) const

inlineprivate

Restore array after finish access local array

Parameters

array

Pointer to the local array fetched by VecGetArray

Definition at line 274 of file PETScVector.cpp.

{
    if (has_ghost_id_)
    {
        VecRestoreArray(v_loc_, &array);
        VecGhostRestoreLocalForm(v_, &v_loc_);
    }
    else
        VecRestoreArray(v_, &array);
}

References has_ghost_id_, v_, and v_loc_.

Referenced by copyValues(), and setLocalAccessibleVector().

set() [1/2]

void MathLib::PETScVector::set ( const PetscInt  i, const PetscScalar  value  )

inline

Insert a single entry with value.

Parameters

i	Entry index
value	Entry value

Definition at line 97 of file PETScVector.h.

    {
        VecSetValue(v_, i, value, INSERT_VALUES);
    }

References v_.

Referenced by MathLib::applyKnownSolution(), and MathLib::setVector().

set() [2/2]

template<class T_SUBVEC >

void MathLib::PETScVector::set ( const std::vector< PetscInt > &  e_idxs, const T_SUBVEC &  sub_vec  )

inline

Add values to several entries

Parameters

e_idxs	Indices of entries to be added. Note: std::size_t cannot be the type of e_idxs template argument
sub_vec	Entries to be added

Definition at line 133 of file PETScVector.h.

    {
        VecSetValues(v_, e_idxs.size(), &e_idxs[0], &sub_vec[0], INSERT_VALUES);
    }

References v_.

setLocalAccessibleVector()

void MathLib::PETScVector::setLocalAccessibleVector

(

)

const

Set local accessible vector in order to get entries. Call this before call operator[] or get(...).

Definition at line 177 of file PETScVector.cpp.

{
    if (entry_array_.empty())
    {
        const PetscInt array_size = global_ids2local_ids_ghost_.empty()
                                        ? size_
                                        : size_loc_ + size_ghosts_;
        entry_array_.resize(array_size);
    }
 
    if (!global_ids2local_ids_ghost_.empty())
    {
        PetscScalar* loc_x = getLocalVector();
        std::copy_n(loc_x, size_loc_ + size_ghosts_, entry_array_.begin());
        restoreArray(loc_x);
    }
    else
        getGlobalVector(entry_array_);
}

References entry_array_, getGlobalVector(), getLocalVector(), global_ids2local_ids_ghost_, restoreArray(), size_, size_ghosts_, and size_loc_.

Referenced by MathLib::LinAlg::setLocalAccessibleVector().

setZero()

void MathLib::PETScVector::setZero ( )

inline

Definition at line 139 of file PETScVector.h.

{ VecSet(v_, 0.0); }

References v_.

shallowCopy()

void MathLib::PETScVector::shallowCopy

(

const PETScVector &

v

)

Definition at line 303 of file PETScVector.cpp.

{
    destroy();
 
    VecDuplicate(v.getRawVector(), &v_);
 
    start_rank_ = v.start_rank_;
    end_rank_ = v.end_rank_;
    size_ = v.size_;
    size_loc_ = v.size_loc_;
    size_ghosts_ = v.size_ghosts_;
    has_ghost_id_ = v.has_ghost_id_;
    global_ids2local_ids_ghost_ = v.global_ids2local_ids_ghost_;
 
    VecSetOption(v_, VEC_IGNORE_NEGATIVE_INDICES, PETSC_TRUE);
}

References destroy(), end_rank_, getRawVector(), global_ids2local_ids_ghost_, has_ghost_id_, size_, size_ghosts_, size_loc_, start_rank_, and v_.

Referenced by PETScVector(), MathLib::LinAlg::copy(), MathLib::LinAlg::matMult(), and MathLib::LinAlg::matMultAdd().

size()

PetscInt MathLib::PETScVector::size ( ) const

inline

Get the global size of the vector.

Definition at line 82 of file PETScVector.h.

{ return size_; }

References size_.

viewer()

void MathLib::PETScVector::viewer	(	const std::string &	file_name,
		const PetscViewerFormat	vw_format = PETSC_VIEWER_ASCII_MATLAB
	)		const

View the global vector for test purpose. Do not use it for output a big vector.

Parameters

file_name

File name for output

vw_format

File format listed as: PETSC_VIEWER_DEFAULT Default format PETSC_VIEWER_ASCII_MATLAB MATLAB format PETSC_VIEWER_ASCII_DENSE Print matrix as dense PETSC_VIEWER_ASCII_IMPL Implementation-specific format (which is in many cases the same as the default) PETSC_VIEWER_ASCII_INFO Basic information about object PETSC_VIEWER_ASCII_INFO_DETAIL More detailed info about object PETSC_VIEWER_ASCII_COMMON Identical output format for all objects of a particular type PETSC_VIEWER_ASCII_INDEX (for vectors) Prints the vector element number next to each vector entry PETSC_VIEWER_ASCII_SYMMODU Print parallel vectors without indicating the processor ranges PETSC_VIEWER_ASCII_VTK Outputs the object to a VTK file PETSC_VIEWER_NATIVE Store the object to the binary file in its native format (for example, dense matrices are stored as dense), DMDA vectors are dumped directly to the file instead of being first put in the natural ordering PETSC_VIEWER_DRAW_BASIC Views the vector with a simple 1d plot PETSC_VIEWER_DRAW_LG Views the vector with a line graph PETSC_VIEWER_DRAW_CONTOUR Views the vector with a contour plot

Definition at line 285 of file PETScVector.cpp.

{
    PetscViewer viewer;
    PetscViewerASCIIOpen(PETSC_COMM_WORLD, file_name.c_str(), &viewer);
    PetscViewerPushFormat(viewer, vw_format);
 
    PetscObjectSetName((PetscObject)v_, file_name.c_str());
    VecView(v_, viewer);
 
#define nEXIT_TEST
#ifdef EXIT_TEST
    VecDestroy(v_);
    PetscFinalize();
    exit(0);
#endif
}

References v_.

Friends And Related Function Documentation

finalizeVectorAssembly

void finalizeVectorAssembly ( PETScVector &  vec )

friend

Function to finalize the vector assembly.

Definition at line 320 of file PETScVector.cpp.

{
    vec.finalizeAssembly();
}

Member Data Documentation

end_rank_

PetscInt MathLib::PETScVector::end_rank_

private

Ending index in a rank.

Definition at line 236 of file PETScVector.h.

Referenced by config(), getRangeEnd(), and shallowCopy().

entry_array_

std::vector<PetscScalar> MathLib::PETScVector::entry_array_

mutableprivate

Array containing the entries of the vector. If the vector is created without given ghost IDs, the array contains all entries of the global vector, v_. Otherwise it only contains the entries of the global vector owned by the current rank.

Definition at line 255 of file PETScVector.h.

Referenced by get(), and setLocalAccessibleVector().

global_ids2local_ids_ghost_

std::map<PetscInt, PetscInt> MathLib::PETScVector::global_ids2local_ids_ghost_

mutableprivate

Map global indices of ghost entries to local indices.

Definition at line 258 of file PETScVector.h.

Referenced by PETScVector(), get(), getLocalIndex(), setLocalAccessibleVector(), and shallowCopy().

has_ghost_id_

bool MathLib::PETScVector::has_ghost_id_ = false

private

Flag to indicate whether the vector is created with ghost entry indices.

Definition at line 246 of file PETScVector.h.

Referenced by getLocalVector(), restoreArray(), and shallowCopy().

size_

PetscInt MathLib::PETScVector::size_

private

Size of the vector.

Definition at line 239 of file PETScVector.h.

Referenced by config(), get(), getGlobalVector(), getLocalIndex(), setLocalAccessibleVector(), shallowCopy(), and size().

size_ghosts_

PetscInt MathLib::PETScVector::size_ghosts_ = 0

private

Size of local ghost entries.

Definition at line 243 of file PETScVector.h.

Referenced by getGhostSize(), setLocalAccessibleVector(), and shallowCopy().

size_loc_

PetscInt MathLib::PETScVector::size_loc_

private

Size of local entries.

Definition at line 241 of file PETScVector.h.

Referenced by PETScVector(), config(), gatherLocalVectors(), getLocalSize(), setLocalAccessibleVector(), and shallowCopy().

start_rank_

PetscInt MathLib::PETScVector::start_rank_

private

Starting index in a rank.

Definition at line 234 of file PETScVector.h.

Referenced by config(), getLocalIndex(), getRangeBegin(), and shallowCopy().

v_

PETSc_Vec MathLib::PETScVector::v_ = nullptr

private

Definition at line 228 of file PETScVector.h.

Referenced by PETScVector(), add(), config(), destroy(), finalizeAssembly(), getGlobalVector(), getLocalVector(), getRawVector(), restoreArray(), set(), setZero(), shallowCopy(), and viewer().

v_loc_

PETSc_Vec MathLib::PETScVector::v_loc_ = nullptr

mutableprivate

Local vector, which is only for the case that v_ is created with ghost entries.

Definition at line 231 of file PETScVector.h.

Referenced by getLocalVector(), and restoreArray().

---

The documentation for this class was generated from the following files:

• MathLib/LinAlg/PETSc/PETScVector.h
• MathLib/LinAlg/PETSc/PETScVector.cpp