libmesh解析: libMesh::EpetraVector< T

在文件 trilinos_epetra_vector.C 第 325 行定义.

 {
   _vec->Abs(*_vec);
 }

template<typename T >

void libMesh::EpetraVector< T >::add	(	const numeric_index_type	i,
		const T	value
	)

overridevirtual

将索引i处的值增加value。

参数

i	索引
value	要增加的值

在文件 trilinos_epetra_vector.C 第 214 行定义.

 {
   int i = static_cast<int> (i_in);
   T value = value_in;
 
   libmesh_assert_less (i_in, this->size());
 
   std::scoped_lock lock(this->_numeric_vector_mutex);
   SumIntoGlobalValues(1, &i, &value);
 
   this->_is_closed = false;
 }

template<typename T >

void libMesh::EpetraVector< T >::add ( const T s )

overridevirtual

将矢量中的所有元素增加s。

参数

s 要增加的值

在文件 trilinos_epetra_vector.C 第 271 行定义.

 {
   const unsigned int nl = _vec->MyLength();
 
   T * values = _vec->Values();
 
   for (unsigned int i=0; i<nl; i++)
     values[i]+=v_in;
 
   this->_is_closed = false;
 }

template<typename T >

void libMesh::EpetraVector< T >::add ( const NumericVector< T > & v )

overridevirtual

将矢量与另一个矢量相加并将结果存储在当前矢量中。

参数

v 另一个矢量

在文件 trilinos_epetra_vector.C 第 285 行定义.

 {
   this->add (1., v);
 }

template<typename T >

void libMesh::EpetraVector< T >::add	(	const T	a,
		const NumericVector< T > &	v
	)

overridevirtual

将矢量与另一个矢量的逐元素乘积与标量a相加并将结果存储在当前矢量中。

参数

a	标量
v	另一个矢量

在文件 trilinos_epetra_vector.C 第 292 行定义.

参考 libMesh::EpetraVector< T >::_vec , 以及 libMesh::EpetraVector< T >::size().

 {
   const EpetraVector<T> * v = cast_ptr<const EpetraVector<T> *>(&v_in);
 
   libmesh_assert_equal_to (this->size(), v->size());
 
   _vec->Update(a_in,*v->_vec, 1.);
 }

template<typename T >

void libMesh::EpetraVector< T >::add_vector	(	const T *	v,
		const std::vector< numeric_index_type > &	dof_indices
	)

overridevirtual

使用特定的自由度索引添加一个值数组到当前矢量中。

参数

v	值数组
dof_indices	自由度索引

在文件 trilinos_epetra_vector.C 第 230 行定义.

参考 libMesh::numeric_trilinos_cast().

 {
   libmesh_assert_equal_to (sizeof(numeric_index_type), sizeof(int));
 
   std::scoped_lock lock(this->_numeric_vector_mutex);
   SumIntoGlobalValues (cast_int<numeric_index_type>(dof_indices.size()),
                        numeric_trilinos_cast(dof_indices.data()),
                        const_cast<T *>(v));
 }

template<typename T >

void libMesh::EpetraVector< T >::add_vector	(	const NumericVector< T > &	v,
		const SparseMatrix< T > &	A
	)

overridevirtual

使用稀疏矩阵A的乘积将另一个矢量v添加到当前矢量中。

参数

v	另一个矢量
A	稀疏矩阵A

在文件 trilinos_epetra_vector.C 第 245 行定义.

参考 libMesh::EpetraVector< T >::_vec , 以及 libMesh::EpetraVector< T >::zero_clone().

 {
   const EpetraVector<T> * v = cast_ptr<const EpetraVector<T> *>(&v_in);
   const EpetraMatrix<T> * A = cast_ptr<const EpetraMatrix<T> *>(&A_in);
 
   // FIXME - does Trilinos let us do this *without* memory allocation?
   std::unique_ptr<NumericVector<T>> temp = v->zero_clone();
   EpetraVector<T> * temp_v = cast_ptr<EpetraVector<T> *>(temp.get());
   A->mat()->Multiply(false, *v->_vec, *temp_v->_vec);
   *this += *temp;
 }

template<typename T >

void libMesh::NumericVector< T >::add_vector	(	const std::vector< T > &	v,
		const std::vector< numeric_index_type > &	dof_indices
	)

inlineinherited

计算 $\vec{u} \leftarrow \vec{u} + \vec{v}$ ，其中 v 是一个 std::vector，每个 dof_indices[i] 指定了要添加的值 v[i] 的位置默认实现对每个索引调用 operator() ，但某些实现可能在此处提供更快的方法。请注意，此方法的库实现是线程安全的。

参数

v	包含要添加的值的 std::vector。
dof_indices	每个要添加值的位置的索引集合。

在文件 numeric_vector.h 第 1041 行定义.

 {
   libmesh_assert(v.size() == dof_indices.size());
   if (!v.empty())
     this->add_vector(v.data(), dof_indices);
 }

template<typename T >

void libMesh::NumericVector< T >::add_vector	(	const NumericVector< T > &	v,
		const std::vector< numeric_index_type > &	dof_indices
	)

inherited

计算 $\vec{u} \leftarrow \vec{u} + \vec{v}$ ，其中 v 是一个 NumericVector，每个 dof_indices[i] 指定了要添加的值 v(i) 的位置。此方法是线程安全的。

参数

v	另一个数值向量。
dof_indices	每个要添加值的位置的索引集合。

在文件 numeric_vector.C 第 391 行定义.

参考 libMesh::NumericVector< T >::readable() , 以及 libMesh::NumericVector< T >::size().

 {
   libmesh_assert(v.readable());
 
   const std::size_t n = dof_indices.size();
   libmesh_assert_equal_to(v.size(), n);
   for (numeric_index_type i=0; i != n; i++)
     this->add (dof_indices[i], v(i));
 }

template<typename T >

void libMesh::NumericVector< T >::add_vector	(	const DenseVector< T > &	v,
		const std::vector< numeric_index_type > &	dof_indices
	)

inlineinherited

计算 $\vec{u} \leftarrow \vec{u} + \vec{v}$ ，其中 v 是一个 DenseVector，每个 dof_indices[i] 指定了要添加的值 v(i) 的位置。此方法是线程安全的。

参数

v	另一个数值向量。
dof_indices	每个要添加值的位置的索引集合。

在文件 numeric_vector.h 第 1053 行定义.

参考 libMesh::DenseVector< T >::empty() , 以及 libMesh::DenseVector< T >::size().

 {
   libmesh_assert(v.size() == dof_indices.size());
   if (!v.empty())
     this->add_vector(&v(0), dof_indices);
 }

template<typename T >

void libMesh::NumericVector< T >::add_vector	(	const NumericVector< T > &	v,
		const ShellMatrix< T > &	A
	)

inherited

计算 $\vec{u} \leftarrow \vec{u} + A \vec{v}$ ，即将 ShellMatrix A 和 NumericVector v 的乘积添加到 this。

参数

v	另一个数值向量。
A	Shell 矩阵。

在文件 numeric_vector.C 第 405 行定义.

参考 libMesh::ShellMatrix< T >::vector_mult_add().

 {
   libmesh_assert(this->compatible(v));
 
   a.vector_mult_add(*this,v);
 }

template<typename T >

void libMesh::EpetraVector< T >::add_vector_transpose	(	const NumericVector< T > &	v,
		const SparseMatrix< T > &	A
	)

overridevirtual

使用稀疏矩阵A的转置乘积将另一个矢量v添加到当前矢量中。

参数

v	另一个矢量
A	稀疏矩阵A

参考自 libMesh::DiagonalMatrix< T >::DiagonalMatrix(), libMesh::DofMap::enforce_adjoint_constraints_exactly(), libMesh::DofMap::enforce_constraints_exactly(), libMesh::DofMap::enforce_constraints_on_residual() , 以及 libMesh::DofMap::process_mesh_constraint_rows().

在文件 trilinos_epetra_vector.C 第 262 行定义.

 {
   libmesh_not_implemented();
 }

template<typename T >

std::unique_ptr< NumericVector< T > > libMesh::NumericVector< T >::build	(	const Parallel::Communicator &	comm,
		const SolverPackage	solver_package = `libMesh::default_solver_package()`
	)

staticinherited

构建一个 NumericVector 对象。

参数

comm	通信器对象，表示通信方式。
solver_package	线性求解器包的类型。

返回: 指向构建的 NumericVector 对象的 unique_ptr。

在文件 numeric_vector.C 第 50 行定义.

参考 libMesh::EIGEN_SOLVERS, libMesh::LASPACK_SOLVERS , 以及 libMesh::TRILINOS_SOLVERS.

 {
   // Build the appropriate vector
   switch (solver_package)
     {
 
 #ifdef LIBMESH_HAVE_LASPACK
     case LASPACK_SOLVERS:
       return std::make_unique<LaspackVector<T>>(comm, AUTOMATIC);
 #endif
 
 #ifdef LIBMESH_HAVE_PETSC
     case PETSC_SOLVERS:
       return std::make_unique<PetscVector<T>>(comm, AUTOMATIC);
 #endif
 
 #ifdef LIBMESH_TRILINOS_HAVE_EPETRA
     case TRILINOS_SOLVERS:
       return std::make_unique<EpetraVector<T>>(comm, AUTOMATIC);
 #endif
 
 #ifdef LIBMESH_HAVE_EIGEN
     case EIGEN_SOLVERS:
       return std::make_unique<EigenSparseVector<T>>(comm, AUTOMATIC);
 #endif
 
     default:
       return std::make_unique<DistributedVector<T>>(comm, AUTOMATIC);
     }
 }

template<typename T >

void libMesh::EpetraVector< T >::clear ( )

inlineoverridevirtualnoexcept

清除矢量，析构函数中会调用此函数，因此不应抛出异常。

在文件 trilinos_epetra_vector.h 第 1065 行定义.

参考 libMesh::libMeshPrivateData::_is_initialized , 以及 libMesh::initialized().

 {
   if (this->initialized())
     {
       // We might just be an interface to a user-provided _vec
       if (this->_destroy_vec_on_exit)
         {
           delete _vec;
           _vec = nullptr;
         }
 
       // But we currently always own our own _map
       _map.reset();
     }
 
   this->_is_closed = this->_is_initialized = false;
 }

template<typename T >

std::unique_ptr< NumericVector< T > > libMesh::EpetraVector< T >::clone ( ) const

inlineoverridevirtual

创建矢量克隆。

返回: 矢量克隆的唯一指针

在文件 trilinos_epetra_vector.h 第 1110 行定义.

参考 libMesh::NumericVector< T >::init().

 {
   NumericVector<T> * cloned_vector = new EpetraVector<T>(this->comm(), AUTOMATIC);
   cloned_vector->init(*this, true);
   *cloned_vector = *this;
   return std::unique_ptr<NumericVector<T>>(cloned_vector);
 }

template<typename T >

void libMesh::EpetraVector< T >::close ( )

inlineoverridevirtual

关闭矢量，使其无法再次修改。

在文件 trilinos_epetra_vector.h 第 1041 行定义.

参考自 libMesh::DofMap::enforce_adjoint_constraints_exactly(), libMesh::DofMap::enforce_constraints_exactly(), libMesh::PetscVector< T >::localize(), libMesh::DofMap::max_constraint_error(), libMesh::EigenSparseVector< T >::operator=(), libMesh::LaspackVector< T >::operator=() , 以及 libMesh::PetscVector< T >::operator=().

参考自 libMesh::EpetraVector< T >::localize().

 {
   libmesh_assert (this->initialized());
 
   // Are we adding or inserting?
   unsigned char global_last_edit = last_edit;
   this->comm().max(global_last_edit);
   libmesh_assert(!last_edit || last_edit == global_last_edit);
 
   if (global_last_edit == 1)
     this->GlobalAssemble(Insert);
   else if (global_last_edit == 2)
     this->GlobalAssemble(Add);
   else
     libmesh_assert(!global_last_edit);
 
   this->_is_closed = true;
   this->last_edit = 0;
 }

template<typename T>

virtual bool libMesh::NumericVector< T >::closed ( ) const

inlinevirtualinherited

检查向量是否已经关闭并准备好进行计算。

返回: 如果向量已经关闭并准备好进行计算，则返回 true；否则返回 false。

在文件 numeric_vector.h 第 180 行定义.

参考 libMesh::NumericVector< T >::_is_closed.

180 { return _is_closed; }

libMesh::NumericVector::_is_closed

bool _is_closed

用于跟踪向量的值在在一些或全部处理器上进行插入或添加值操作后是否在所有处理器上保持一致的标志。

Definition: numeric_vector.h:913

template<typename T >

int libMesh::NumericVector< T >::compare	(	const NumericVector< T > &	other_vector,
		const Real	threshold = `TOLERANCE`
	)		const

virtualinherited

比较 this 与 other_vector 的等效性，（在给定的 threshold 内）如果等效则返回 -1 ，或者返回第一个索引，其中 abs(a[i]-b[i]) 超过阈值。

注解: "等效" 表示两个向量非常相似或几乎相同，其差异小于或等于给定的阈值。这个阈值通常用于考虑浮点数精度问题，因为在计算机中，浮点数运算可能会导致微小的数值差异。当两个向量被认为等效时，它们的差异小于或等于阈值。具体来说，这是通过计算两个向量的每个元素之间的差异并将其与阈值进行比较来确定的。如果两个向量的差异小于阈值，就可以说它们是等效的。否则，如果差异大于阈值，就可以说它们不等效。

参数

other_vector	另一个数值向量。
threshold	阈值。

返回: 如果向量等效则返回 -1 ，否则返回第一个不等效的索引。

在文件 numeric_vector.C 第 109 行定义.

参考 std::abs().

 {
   libmesh_assert(this->compatible(other_vector));
 
   int first_different_i = std::numeric_limits<int>::max();
   numeric_index_type i = first_local_index();
 
   while (first_different_i==std::numeric_limits<int>::max()
          && i<last_local_index())
   {
     if (std::abs((*this)(i) - other_vector(i)) > threshold)
       first_different_i = i;
     else
       i++;
   }
 
   // Find the correct first differing index in parallel
   this->comm().min(first_different_i);
 
   if (first_different_i == std::numeric_limits<int>::max())
     return -1;
 
   return first_different_i;
 }

template<typename T >

bool libMesh::NumericVector< T >::compatible ( const NumericVector< T > & v ) const

inherited

检查该向量和向量 v 是否能够一起用于全局操作。

参数

v	要比较的另一个向量。

返回: 如果这两个向量可以一起用于全局操作，则返回 true，否则返回 false。

在文件 numeric_vector.C 第 423 行定义.

参考 libMesh::NumericVector< T >::first_local_index(), libMesh::NumericVector< T >::last_local_index(), libMesh::NumericVector< T >::local_size(), libMesh::NumericVector< T >::readable() , 以及 libMesh::NumericVector< T >::size().

 {
   return this->readable() && v.readable() &&
          this->size() == v.size() &&
          this->local_size() == v.local_size() &&
          this->first_local_index() == v.first_local_index() &&
          this->last_local_index() == v.last_local_index();
 }

template<typename T >

void libMesh::EpetraVector< T >::conjugate ( )

overridevirtual

计算当前矢量的逐元素共轭。

在文件 trilinos_epetra_vector.C 第 206 行定义.

 {
   // EPetra is real, rendering this a no-op.
 }

template<typename T >

void libMesh::EpetraVector< T >::create_subvector	(	NumericVector< T > &	subvector,
		const std::vector< numeric_index_type > &	rows
	)		const

overridevirtual

创建当前矢量的子矢量，包含指定行的值。

参数

subvector	子矢量
rows	指定行的索引列表

在文件 trilinos_epetra_vector.C 第 581 行定义.

 {
   libmesh_not_implemented();
 }

template<typename T >

void libMesh::EpetraVector< T >::destroyNonlocalData ( )

private

销毁非本地数据。

在文件 trilinos_epetra_vector.C 第 915 行定义.

 {
   if (allocatedNonlocalLength_ > 0) {
     delete [] nonlocalIDs_;
     delete [] nonlocalElementSize_;
     nonlocalIDs_ = nullptr;
     nonlocalElementSize_ = nullptr;
     for (int i=0; i<numNonlocalIDs_; ++i) {
       delete [] nonlocalCoefs_[i];
     }
     delete [] nonlocalCoefs_;
     nonlocalCoefs_ = nullptr;
     numNonlocalIDs_ = 0;
     allocatedNonlocalLength_ = 0;
   }
   return;
 }

void libMesh::ReferenceCounter::disable_print_counter_info ( )

staticinherited

在文件 reference_counter.C 第 100 行定义.

参考 libMesh::ReferenceCounter::_enable_print_counter.

 {
   _enable_print_counter = false;
   return;
 }

template<typename T >

T libMesh::EpetraVector< T >::dot ( const NumericVector< T > & v ) const

overridevirtual

计算当前矢量与另一个矢量的点积。

参数

v 另一个矢量

返回: 点积结果

在文件 trilinos_epetra_vector.C 第 332 行定义.

参考 libMesh::EpetraVector< T >::_vec, libMesh::EpetraVector< T >::allocatedNonlocalLength_, libMesh::EpetraVector< T >::nonlocalCoefs_, libMesh::EpetraVector< T >::nonlocalElementSize_, libMesh::EpetraVector< T >::nonlocalIDs_ , 以及 libMesh::EpetraVector< T >::numNonlocalIDs_.

 {
   const EpetraVector<T> * v = cast_ptr<const EpetraVector<T> *>(&v_in);
 
   T result=0.0;
 
   _vec->Dot(*v->_vec, &result);
 
   return result;
 }

template<typename T>

virtual T libMesh::NumericVector< T >::el ( const numeric_index_type i ) const

inlinevirtualinherited

获取向量的第 i 个条目。

参数

i	要获取的条目的索引。

返回: 第 i 个条目的值。

在文件 numeric_vector.h 第 416 行定义.

416 { return (*this)(i); }

void libMesh::ReferenceCounter::enable_print_counter_info ( )

staticinherited

Methods to enable/disable the reference counter output from print_info()

在文件 reference_counter.C 第 94 行定义.

参考 libMesh::ReferenceCounter::_enable_print_counter.

 {
   _enable_print_counter = true;
   return;
 }

template<typename T >

void libMesh::EpetraVector< T >::FEoperatorequals ( const EpetraVector< T > & source )

private

从另一个 EpetraVector 对象中复制操作符。

参数

source 源 EpetraVector 对象。

在文件 trilinos_epetra_vector.C 第 887 行定义.

 {
   (*_vec) = *(source._vec);
 
   destroyNonlocalData();
 
   if (source.allocatedNonlocalLength_ > 0) {
     allocatedNonlocalLength_ = source.allocatedNonlocalLength_;
     numNonlocalIDs_ = source.numNonlocalIDs_;
     nonlocalIDs_ = new int[allocatedNonlocalLength_];
     nonlocalElementSize_ = new int[allocatedNonlocalLength_];
     nonlocalCoefs_ = new double *[allocatedNonlocalLength_];
     for (int i=0; i<numNonlocalIDs_; ++i) {
       int elemSize = source.nonlocalElementSize_[i];
       nonlocalCoefs_[i] = new double[elemSize];
       nonlocalIDs_[i] = source.nonlocalIDs_[i];
       nonlocalElementSize_[i] = elemSize;
       for (int j=0; j<elemSize; ++j) {
         nonlocalCoefs_[i][j] = source.nonlocalCoefs_[i][j];
       }
     }
   }
 }

template<typename T >

numeric_index_type libMesh::EpetraVector< T >::first_local_index ( ) const

inlineoverridevirtual

获取矢量的第一个本地索引。

返回: 第一个本地索引

在文件 trilinos_epetra_vector.h 第 1142 行定义.

 {
   libmesh_assert (this->initialized());
 
   return _vec->Map().MinMyGID();
 }

template<typename T >

void libMesh::NumericVector< T >::get	(	const std::vector< numeric_index_type > &	index,
		T *	values
	)		const

inlinevirtualinherited

一次访问多个组件。 values 将 *不会* 重新分配空间；它应该已经具有足够的空间。默认实现对每个索引调用 operator() ，但某些实现可能在此处提供更快的方法。

参数

index	要获取的索引集合。
values	存储获取值的数组。

被 libMesh::PetscVector< T > 重载.

在文件 numeric_vector.h 第 1012 行定义.

参考自 libMesh::DofMap::enforce_adjoint_constraints_exactly(), libMesh::DofMap::enforce_constraints_exactly(), libMesh::DofMap::enforce_constraints_on_residual() , 以及 libMesh::DofMap::max_constraint_error().

 {
   const std::size_t num = index.size();
   for (std::size_t i=0; i<num; i++)
     {
       values[i] = (*this)(index[i]);
     }
 }

template<typename T >

void libMesh::NumericVector< T >::get	(	const std::vector< numeric_index_type > &	index,
		std::vector< T > &	values
	)		const

inlineinherited

一次访问多个组件。 values 将被调整大小，如果需要，将被填充。默认实现对每个索引调用 operator() ，但某些实现可能在此处提供更快的方法。

参数

index	要获取的索引集合。
values	存储获取值的向量。

在文件 numeric_vector.h 第 1026 行定义.

 {
   const std::size_t num = index.size();
   values.resize(num);
   if (!num)
     return;
 
   this->get(index, values.data());
 }

std::string libMesh::ReferenceCounter::get_info ( )

staticinherited

Gets a string containing the reference information.

在文件 reference_counter.C 第 47 行定义.

参考 libMesh::ReferenceCounter::_counts.

参考自 libMesh::ReferenceCounter::print_info().

 {
 #if defined(LIBMESH_ENABLE_REFERENCE_COUNTING) && defined(DEBUG)
 
   std::ostringstream oss;
 
   oss << '\n'
       << " ---------------------------------------------------------------------------- \n"
       << "| Reference count information                                                |\n"
       << " ---------------------------------------------------------------------------- \n";
 
   for (const auto & [name, cd] : _counts)
     oss << "| " << name << " reference count information:\n"
         << "|  Creations:    " << cd.first    << '\n'
         << "|  Destructions: " << cd.second << '\n';
 
   oss << " ---------------------------------------------------------------------------- \n";
 
   return oss.str();
 
 #else
 
   return "";
 
 #endif
 }

template<typename T >

int libMesh::EpetraVector< T >::GlobalAssemble ( Epetra_CombineMode mode = Add )

private

将所有重叠/共享数据收集到由 Map 在构造函数中传递给该向量定义的非重叠分区中。从其他处理器导入的数据以“sumInto”或累积操作存储在拥有处理器上。这是一种集体方法，每个处理器在任何处理器完成它之前都必须进入它。

参数

mode	聚合模式（默认为 Add）。

返回: 操作的结果代码。

在文件 trilinos_epetra_vector.C 第 841 行定义.

 {
   //In this method we need to gather all the non-local (overlapping) data
   //that's been input on each processor, into the (probably) non-overlapping
   //distribution defined by the map that 'this' vector was constructed with.
 
   //We don't need to do anything if there's only one processor or if
   //ignoreNonLocalEntries_ is true.
   if (_vec->Map().Comm().NumProc() < 2 || ignoreNonLocalEntries_) {
     return(0);
   }
 
 
 
   //First build a map that describes the data in nonlocalIDs_/nonlocalCoefs_.
   //We'll use the arbitrary distribution constructor of Map.
 
   Epetra_BlockMap sourceMap(-1, numNonlocalIDs_,
                             nonlocalIDs_, nonlocalElementSize_,
                             _vec->Map().IndexBase(), _vec->Map().Comm());
 
   //Now build a vector to hold our nonlocalCoefs_, and to act as the source-
   //vector for our import operation.
   Epetra_MultiVector nonlocalVector(sourceMap, 1);
 
   int i,j;
   for (i=0; i<numNonlocalIDs_; ++i) {
     for (j=0; j<nonlocalElementSize_[i]; ++j) {
       nonlocalVector.ReplaceGlobalValue(nonlocalIDs_[i], j, 0,
                                         nonlocalCoefs_[i][j]);
     }
   }
 
   Epetra_Export exporter(sourceMap, _vec->Map());
 
   EPETRA_CHK_ERR( _vec->Export(nonlocalVector, exporter, mode) );
 
   destroyNonlocalData();
 
   return(0);
 }

void libMesh::ReferenceCounter::increment_constructor_count ( const std::string & name )

inlineprotectednoexceptinherited

Increments the construction counter.

Should be called in the constructor of any derived class that will be reference counted.

在文件 reference_counter.h 第 183 行定义.

参考 libMesh::err.

参考自 libMesh::ReferenceCountedObject< SparseMatrix< T > >::ReferenceCountedObject().

 {
   libmesh_try
   {
     Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
     std::pair<unsigned int, unsigned int> & p = _counts[name];
     p.first++;
   }
   libmesh_catch (...)
   {
     auto stream = libMesh::err.get();
     stream->exceptions(stream->goodbit); // stream must not throw
     libMesh::err << "Encountered unrecoverable error while calling "
                  << "ReferenceCounter::increment_constructor_count() "
                  << "for a(n) " << name << " object." << std::endl;
     std::terminate();
   }
 }

void libMesh::ReferenceCounter::increment_destructor_count ( const std::string & name )

inlineprotectednoexceptinherited

Increments the destruction counter.

Should be called in the destructor of any derived class that will be reference counted.

在文件 reference_counter.h 第 207 行定义.

参考 libMesh::err.

参考自 libMesh::ReferenceCountedObject< SparseMatrix< T > >::~ReferenceCountedObject().

 {
   libmesh_try
   {
     Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
     std::pair<unsigned int, unsigned int> & p = _counts[name];
     p.second++;
   }
   libmesh_catch (...)
   {
     auto stream = libMesh::err.get();
     stream->exceptions(stream->goodbit); // stream must not throw
     libMesh::err << "Encountered unrecoverable error while calling "
                  << "ReferenceCounter::increment_destructor_count() "
                  << "for a(n) " << name << " object." << std::endl;
     std::terminate();
   }
 }

template<typename T >

void libMesh::EpetraVector< T >::init	(	const numeric_index_type	N,
		const numeric_index_type	n_local,
		const bool	fast = `false`,
		const ParallelType	type = `AUTOMATIC`
	)

inlineoverridevirtual

初始化矢量。

参数

N	全局维度
n_local	本地维度
fast	是否启用快速初始化
type	并行类型（默认为AUTOMATIC）

在文件 trilinos_epetra_vector.h 第 960 行定义.

参考 libMesh::libMeshPrivateData::_is_initialized , 以及 libMesh::zero.

参考自 libMesh::EpetraVector< T >::EpetraVector() , 以及 libMesh::EpetraVector< T >::localize().

 {
   // We default to allocating n_local local storage
   numeric_index_type my_n_local = n_local;
 
   if (type == AUTOMATIC)
     {
       if (n == n_local)
         this->_type = SERIAL;
       else
         this->_type = PARALLEL;
     }
   else if (type == GHOSTED)
     {
       // We don't yet support GHOSTED Epetra vectors, so to get the
       // same functionality we need a SERIAL vector with local
       // storage allocated for every entry.
       this->_type = SERIAL;
       my_n_local = n;
     }
   else
     this->_type = type;
 
   libmesh_assert ((this->_type==SERIAL && n==my_n_local) ||
                   this->_type==PARALLEL);
 
   _map = std::make_unique<Epetra_Map>
     (static_cast<int>(n),
      my_n_local,
      0,
      Epetra_MpiComm (this->comm().get()));
 
   _vec = new Epetra_Vector(*_map);
 
   myFirstID_ = _vec->Map().MinMyGID();
   myNumIDs_ = _vec->Map().NumMyElements();
 
   // Currently we impose the restriction that NumVectors==1, so we won't
   // need the LDA argument when calling ExtractView. Hence the "dummy" arg.
   int dummy;
   _vec->ExtractView(&myCoefs_, &dummy);
 
   this->_is_initialized = true;
   this->_is_closed = true;
   this->last_edit = 0;
 
   if (fast == false)
     this->zero ();
 }

template<typename T >

void libMesh::EpetraVector< T >::init	(	const numeric_index_type	N,
		const bool	fast = `false`,
		const ParallelType	type = `AUTOMATIC`
	)

inlineoverridevirtual

初始化矢量。

参数

N	全局维度
fast	是否启用快速初始化
type	并行类型（默认为AUTOMATIC）

在文件 trilinos_epetra_vector.h 第 1030 行定义.

 {
   this->init(n,n,fast,type);
 }

template<typename T >

void libMesh::EpetraVector< T >::init	(	const numeric_index_type	N,
		const numeric_index_type	n_local,
		const std::vector< numeric_index_type > &	ghost,
		const bool	fast = `false`,
		const ParallelType	type = `AUTOMATIC`
	)

inlineoverridevirtual

初始化矢量。

参数

N	全局维度
n_local	本地维度
ghost	需要保留内存的索引
fast	是否启用快速初始化
type	并行类型（默认为AUTOMATIC）

在文件 trilinos_epetra_vector.h 第 1016 行定义.

 {
   // TODO: we shouldn't ignore the ghost sparsity pattern
   this->init(n, n_local, fast, type);
 }

template<class T >

void libMesh::EpetraVector< T >::init	(	const NumericVector< T > &	other,
		const bool	fast = `false`
	)

overridevirtual

初始化矢量，使用其他矢量进行初始化。

参数

other	源矢量
fast	是否启用快速初始化

libMesh::NumericVector::_is_initialized

在文件 trilinos_epetra_vector.h 第 941 行定义.

参考 libMesh::NumericVector< T >::local_size(), libMesh::NumericVector< T >::size() , 以及 libMesh::NumericVector< T >::type().

 {
   this->init(other.size(),other.local_size(),fast,other.type());
 }

template<typename T>

virtual bool libMesh::NumericVector< T >::initialized ( ) const

inlinevirtualinherited

检查向量是否已经初始化。

返回: 如果向量已经初始化，则返回 true；否则返回 false。

在文件 numeric_vector.h 第 159 行定义.

参考 libMesh::NumericVector< T >::_is_initialized.

参考自 libMesh::PetscVector< T >::create_subvector() , 以及 libMesh::PetscVector< T >::init().

159 { return _is_initialized; }

bool _is_initialized

在调用 init() 后设置为 true。

Definition: numeric_vector.h:918

template<typename T >

int libMesh::EpetraVector< T >::inputNonlocalValue	(	int	GID,
		double	value,
		bool	accumulate
	)

private

输入非本地值到向量中，覆盖或累积到指定 GID 的已存在的任何值。

参数

GID	全局 ID。
value	输入的值。
accumulate	如果为 true，则将值累积到已存在的值中。

返回: 操作的结果代码。

在文件 trilinos_epetra_vector.C 第 739 行定义.

 {
   int insertPoint = -1;
 
   //find offset of GID in nonlocalIDs_
   int offset = Epetra_Util_binary_search(GID, nonlocalIDs_, numNonlocalIDs_,
                                          insertPoint);
   if (offset >= 0) {
     //if offset >= 0
     //  put value in nonlocalCoefs_[offset][0]
 
     if (accumulate) {
       nonlocalCoefs_[offset][0] += value;
     }
     else {
       nonlocalCoefs_[offset][0] = value;
     }
   }
   else {
     //else
     //  insert GID in nonlocalIDs_
     //  insert 1   in nonlocalElementSize_
     //  insert value in nonlocalCoefs_
 
     int tmp1 = numNonlocalIDs_;
     int tmp2 = allocatedNonlocalLength_;
     int tmp3 = allocatedNonlocalLength_;
     EPETRA_CHK_ERR( Epetra_Util_insert(GID, insertPoint, nonlocalIDs_,
                                        tmp1, tmp2) );
     --tmp1;
     EPETRA_CHK_ERR( Epetra_Util_insert(1, insertPoint, nonlocalElementSize_,
                                        tmp1, tmp3) );
     double * values = new double[1];
     values[0] = value;
     EPETRA_CHK_ERR( Epetra_Util_insert(values, insertPoint, nonlocalCoefs_,
                                        numNonlocalIDs_, allocatedNonlocalLength_) );
   }
 
   return(0);
 }

template<typename T >

int libMesh::EpetraVector< T >::inputNonlocalValues	(	int	GID,
		int	numValues,
		const double *	values,
		bool	accumulate
	)

private

输入非本地值到向量中，覆盖或累积到指定 GID 的已存在的任何值。

参数

GID	全局 ID。
numValues	值的数量。
values	包含系数值的数组。
accumulate	如果为 true，则将值累积到已存在的值中。

返回: 操作的结果代码。

在文件 trilinos_epetra_vector.C 第 782 行定义.

参考 libMesh::err.

 {
   int insertPoint = -1;
 
   //find offset of GID in nonlocalIDs_
   int offset = Epetra_Util_binary_search(GID, nonlocalIDs_, numNonlocalIDs_,
                                          insertPoint);
   if (offset >= 0) {
     //if offset >= 0
     //  put value in nonlocalCoefs_[offset][0]
 
     if (numValues != nonlocalElementSize_[offset]) {
       libMesh::err << "Epetra_FEVector ERROR: block-size for GID " << GID << " is "
                    << numValues<<" which doesn't match previously set block-size of "
                    << nonlocalElementSize_[offset] << std::endl;
       return(-1);
     }
 
     if (accumulate) {
       for (int j=0; j<numValues; ++j) {
         nonlocalCoefs_[offset][j] += values[j];
       }
     }
     else {
       for (int j=0; j<numValues; ++j) {
         nonlocalCoefs_[offset][j] = values[j];
       }
     }
   }
   else {
     //else
     //  insert GID in nonlocalIDs_
     //  insert numValues   in nonlocalElementSize_
     //  insert values in nonlocalCoefs_
 
     int tmp1 = numNonlocalIDs_;
     int tmp2 = allocatedNonlocalLength_;
     int tmp3 = allocatedNonlocalLength_;
     EPETRA_CHK_ERR( Epetra_Util_insert(GID, insertPoint, nonlocalIDs_,
                                        tmp1, tmp2) );
     --tmp1;
     EPETRA_CHK_ERR( Epetra_Util_insert(numValues, insertPoint, nonlocalElementSize_,
                                        tmp1, tmp3) );
     double * newvalues = new double[numValues];
     for (int j=0; j<numValues; ++j) {
       newvalues[j] = values[j];
     }
     EPETRA_CHK_ERR( Epetra_Util_insert(newvalues, insertPoint, nonlocalCoefs_,
                                        numNonlocalIDs_, allocatedNonlocalLength_) );
   }
 
   return(0);
 }

template<typename T >

int libMesh::EpetraVector< T >::inputValues	(	int	numIDs,
		const int *	GIDs,
		const double *	values,
		bool	accumulate
	)

private

输入值到向量中，覆盖或累积到指定索引的已存在的任何值。

参数

numIDs	全局 ID 的数量。
GIDs	包含全局 ID 的数组。
values	包含系数值的数组。
accumulate	如果为 true，则将值累积到已存在的值中。

返回: 操作的结果代码。

在文件 trilinos_epetra_vector.C 第 659 行定义.

 {
   if (accumulate) {
     libmesh_assert(last_edit == 0 || last_edit == 2);
     last_edit = 2;
   } else {
     libmesh_assert(last_edit == 0 || last_edit == 1);
     last_edit = 1;
   }
 
   //Important note!! This method assumes that there is only 1 point
   //associated with each element.
 
   for (int i=0; i<numIDs; ++i) {
     if (_vec->Map().MyGID(GIDs[i])) {
       if (accumulate) {
         _vec->SumIntoGlobalValue(GIDs[i], 0, 0, values[i]);
       }
       else {
         _vec->ReplaceGlobalValue(GIDs[i], 0, 0, values[i]);
       }
     }
     else {
       if (!ignoreNonLocalEntries_) {
         EPETRA_CHK_ERR( inputNonlocalValue(GIDs[i], values[i], accumulate) );
       }
     }
   }
 
   return(0);
 }

template<typename T >

int libMesh::EpetraVector< T >::inputValues	(	int	numIDs,
		const int *	GIDs,
		const int *	numValuesPerID,
		const double *	values,
		bool	accumulate
	)

private

输入值到向量中，替换或累积到指定全局 ID 的已存在的任何值。

参数

numIDs	全局 ID 的数量。
GIDs	包含全局 ID 的数组。
numValuesPerID	包含每个全局 ID 对应值的数量的数组。
values	包含系数值的数组。
accumulate	如果为 true，则将值累积到已存在的值中。

返回: 操作的结果代码。

在文件 trilinos_epetra_vector.C 第 696 行定义.

 {
   if (accumulate) {
     libmesh_assert(last_edit == 0 || last_edit == 2);
     last_edit = 2;
   } else {
     libmesh_assert(last_edit == 0 || last_edit == 1);
     last_edit = 1;
   }
 
   int offset=0;
   for (int i=0; i<numIDs; ++i) {
     int numValues = numValuesPerID[i];
     if (_vec->Map().MyGID(GIDs[i])) {
       if (accumulate) {
         for (int j=0; j<numValues; ++j) {
           _vec->SumIntoGlobalValue(GIDs[i], j, 0, values[offset+j]);
         }
       }
       else {
         for (int j=0; j<numValues; ++j) {
           _vec->ReplaceGlobalValue(GIDs[i], j, 0, values[offset+j]);
         }
       }
     }
     else {
       if (!ignoreNonLocalEntries_) {
         EPETRA_CHK_ERR( inputNonlocalValues(GIDs[i], numValues,
                                             &(values[offset]), accumulate) );
       }
     }
     offset += numValues;
   }
 
   return(0);
 }

template<typename T >

void libMesh::EpetraVector< T >::insert	(	const T *	v,
		const std::vector< numeric_index_type > &	dof_indices
	)

overridevirtual

使用特定的自由度索引将一个值数组插入到当前矢量中。

参数

v	值数组
dof_indices	自由度索引

在文件 trilinos_epetra_vector.C 第 304 行定义.

参考 libMesh::numeric_trilinos_cast().

 {
   libmesh_assert_equal_to (sizeof(numeric_index_type), sizeof(int));
 
   std::scoped_lock lock(this->_numeric_vector_mutex);
   ReplaceGlobalValues (cast_int<numeric_index_type>(dof_indices.size()),
                        numeric_trilinos_cast(dof_indices.data()),
                        const_cast<T *>(v));
   this->_is_closed = false;
 }

template<typename T >

void libMesh::NumericVector< T >::insert	(	const NumericVector< T > &	v,
		const std::vector< numeric_index_type > &	dof_indices
	)

inherited

将 v 的条目插入到 *this 中，位置由 v 指定。只要使用 NumericVector 的库实现，此方法是线程安全的。

参数

v	一个数值向量。
dof_indices	每个要插入值的位置的索引集合。

在文件 numeric_vector.C 第 96 行定义.

参考 libMesh::NumericVector< T >::readable() , 以及 libMesh::NumericVector< T >::size().

 {
   libmesh_assert_equal_to (V.size(), dof_indices.size());
   libmesh_assert (V.readable());
 
   for (auto i : index_range(dof_indices))
     this->set (dof_indices[i], V(i));
 }

template<typename T >

void libMesh::NumericVector< T >::insert	(	const DenseVector< T > &	v,
		const std::vector< numeric_index_type > &	dof_indices
	)

inlineinherited

将 v 的条目插入到 *this 中，位置由 v 指定。只要使用 NumericVector 的库实现，此方法是线程安全的。

参数

v	一个 DenseVector。
dof_indices	每个要插入值的位置的索引集合。

在文件 numeric_vector.h 第 1077 行定义.

参考 libMesh::DenseVector< T >::empty() , 以及 libMesh::DenseVector< T >::size().

 {
   libmesh_assert(v.size() == dof_indices.size());
   if (!v.empty())
     this->insert(&v(0), dof_indices);
 }

template<typename T >

void libMesh::NumericVector< T >::insert	(	const DenseSubVector< T > &	v,
		const std::vector< numeric_index_type > &	dof_indices
	)

inlineinherited

将 v 的条目插入到 *this 中，位置由 v 指定。只要使用 NumericVector 的库实现，此方法是线程安全的。

参数

v	一个 DenseSubVector。
dof_indices	每个要插入值的位置的索引集合。

在文件 numeric_vector.h 第 1089 行定义.

参考 libMesh::DenseSubVector< T >::empty() , 以及 libMesh::DenseSubVector< T >::size().

 {
   libmesh_assert(v.size() == dof_indices.size());
   if (!v.empty())
     this->insert(&v(0), dof_indices);
 }

template<typename T >

Real libMesh::EpetraVector< T >::l1_norm ( ) const

overridevirtual

计算矢量的L1范数。

返回: L1范数

在文件 trilinos_epetra_vector.C 第 69 行定义.

参考 libMesh::closed() , 以及 libMesh::Real.

 {
   libmesh_assert(this->closed());
 
   Real value;
 
   _vec->Norm1(&value);
 
   return value;
 }

template<typename T >

Real libMesh::EpetraVector< T >::l2_norm ( ) const

overridevirtual

计算矢量的L2范数。

返回: L2范数

在文件 trilinos_epetra_vector.C 第 81 行定义.

参考 libMesh::closed() , 以及 libMesh::Real.

 {
   libmesh_assert(this->closed());
 
   Real value;
 
   _vec->Norm2(&value);
 
   return value;
 }

template<class T >

Real libMesh::NumericVector< T >::l2_norm_diff ( const NumericVector< T > & other_vec ) const

inherited

获取 $\ell_2$ -范数的向量差值 $\vec{u} - \vec{v}$ ，其中 $\vec{u}$ 是 this。

参数

other_vec 另一个数值向量。

返回: $\ell_2$ -范数的向量差值。

在文件 numeric_vector.C 第 363 行定义.

参考 std::norm(), libMesh::TensorTools::norm_sq(), libMesh::Real , 以及 std::sqrt().

 {
   libmesh_assert(this->compatible(v));
 
   Real norm = 0;
   for (const auto i : make_range(this->first_local_index(), this->last_local_index()))
     norm += TensorTools::norm_sq((*this)(i) - v(i));
 
   this->comm().sum(norm);
 
   return std::sqrt(norm);
 }

template<typename T >

numeric_index_type libMesh::EpetraVector< T >::last_local_index ( ) const

inlineoverridevirtual

获取矢量的最后一个本地索引。

返回: 最后一个本地索引

在文件 trilinos_epetra_vector.h 第 1153 行定义.

 {
   libmesh_assert (this->initialized());
 
   return _vec->Map().MaxMyGID()+1;
 }

template<typename T >

Real libMesh::EpetraVector< T >::linfty_norm ( ) const

overridevirtual

计算矢量的L∞范数。

返回: L∞范数

在文件 trilinos_epetra_vector.C 第 93 行定义.

参考 libMesh::closed() , 以及 libMesh::Real.

 {
   libmesh_assert(this->closed());
 
   Real value;
 
   _vec->NormInf(&value);
 
   return value;
 }

template<typename T >

int libMesh::NumericVector< T >::local_relative_compare	(	const NumericVector< T > &	other_vector,
		const Real	threshold = `TOLERANCE`
	)		const

virtualinherited

比较该向量与另一个向量的局部相对差异。

如果两个向量在给定的局部相对阈值 abs(a[i]-b[i])/max(a[i],b[i]) 下等效，则返回-1，否则返回首个不等效的索引。

参数

other_vector	要比较的另一个向量。
threshold	允许的局部相对差异阈值（可选，默认值为 TOLERANCE）。

返回: -1 表示等效，否则返回不等效的索引。

在文件 numeric_vector.C 第 137 行定义.

参考 std::abs().

 {
   libmesh_assert(this->compatible(other_vector));
 
   int first_different_i = std::numeric_limits<int>::max();
   numeric_index_type i = first_local_index();
 
   do
     {
       if (std::abs((*this)(i) - other_vector(i)) > threshold *
           std::max(std::abs((*this)(i)), std::abs(other_vector(i))))
         first_different_i = i;
       else
         i++;
     }
   while (first_different_i==std::numeric_limits<int>::max()
          && i<last_local_index());
 
   // Find the correct first differing index in parallel
   this->comm().min(first_different_i);
 
   if (first_different_i == std::numeric_limits<int>::max())
     return -1;
 
   return first_different_i;
 }

template<typename T >

numeric_index_type libMesh::EpetraVector< T >::local_size ( ) const

inlineoverridevirtual

获取矢量的本地大小。

返回: 矢量的本地大小

在文件 trilinos_epetra_vector.h 第 1133 行定义.

 {
   libmesh_assert (this->initialized());
 
   return _vec->MyLength();
 }

template<typename T >

void libMesh::EpetraVector< T >::localize ( std::vector< T > & v_local ) const

overridevirtual

将当前矢量的值本地化到一个本地矢量中。

参数

v_local 存储本地化值的矢量

在文件 trilinos_epetra_vector.C 第 532 行定义.

参考自 libMesh::EpetraVector< T >::localize().

 {
   // This function must be run on all processors at once
   parallel_object_only();
 
   const unsigned int n  = this->size();
   const unsigned int nl = this->local_size();
 
   libmesh_assert(this->_vec);
 
   v_local.clear();
   v_local.reserve(n);
 
   // build up my local part
   for (unsigned int i=0; i<nl; i++)
     v_local.push_back((*this->_vec)[i]);
 
   this->comm().allgather (v_local);
 }

template<typename T >

void libMesh::EpetraVector< T >::localize ( NumericVector< T > & v_local ) const

overridevirtual

将当前矢量的值本地化到另一个矢量的相应位置。

参数

v_local 存储本地化值的矢量

在文件 trilinos_epetra_vector.C 第 420 行定义.

 {
   EpetraVector<T> * v_local = cast_ptr<EpetraVector<T> *>(&v_local_in);
 
   Epetra_Map rootMap = Epetra_Util::Create_Root_Map( *_map, -1);
   v_local->_vec->ReplaceMap(rootMap);
 
   Epetra_Import importer(v_local->_vec->Map(), *_map);
   v_local->_vec->Import(*_vec, importer, Insert);
 }

template<typename T >

void libMesh::EpetraVector< T >::localize	(	NumericVector< T > &	v_local,
		const std::vector< numeric_index_type > &	send_list
	)		const

overridevirtual

将当前矢量的部分值本地化到另一个矢量的相应位置。

参数

v_local	存储本地化值的矢量
send_list	本地化的部分值的索引列表

在文件 trilinos_epetra_vector.C 第 434 行定义.

 {
   // TODO: optimize to sync only the send list values
   this->localize(v_local_in);
 
   //   EpetraVector<T> * v_local =
   //   cast_ptr<EpetraVector<T> *>(&v_local_in);
 
   //   libmesh_assert(this->_map.get());
   //   libmesh_assert(v_local->_map.get());
   //   libmesh_assert_equal_to (v_local->local_size(), this->size());
   //   libmesh_assert_less_equal (send_list.size(), v_local->size());
 
   //   Epetra_Import importer (*v_local->_map, *this->_map);
 
   //   v_local->_vec->Import (*this->_vec, importer, Insert);
 }

template<typename T >

void libMesh::EpetraVector< T >::localize	(	std::vector< T > &	v_local,
		const std::vector< numeric_index_type > &	indices
	)		const

overridevirtual

将当前矢量的部分值本地化到本地矢量的相应位置。

参数

v_local	存储本地化值的矢量
indices	本地化的部分值的索引列表

在文件 trilinos_epetra_vector.C 第 456 行定义.

 {
   // Create a "replicated" map for importing values.  This is
   // equivalent to creating a general Epetra_Map with
   // NumGlobalElements == NumMyElements.
   Epetra_LocalMap import_map(static_cast<int>(indices.size()),
                              /*IndexBase=*/0,
                              _map->Comm());
 
   // Get a pointer to the list of global elements for the map, and set
   // all the values from indices.
   int * import_map_global_elements = import_map.MyGlobalElements();
   for (auto i : index_range(indices))
     import_map_global_elements[i] = indices[i];
 
   // Create a new EpetraVector to import values into.
   Epetra_Vector import_vector(import_map);
 
   // Set up an "Import" object which associates the two maps.
   Epetra_Import import_object(import_map, *_map);
 
   // Import the values
   import_vector.Import(*_vec, import_object, Insert);
 
   // Get a pointer to the imported values array and the length of the
   // array.
   T * values = import_vector.Values();
   int import_vector_length = import_vector.MyLength();
 
   // Copy the imported values into v_local
   v_local.resize(import_vector_length);
   for (int i=0; i<import_vector_length; ++i)
     v_local[i] = values[i];
 }

template<typename T >

void libMesh::EpetraVector< T >::localize	(	const numeric_index_type	first_local_idx,
		const numeric_index_type	last_local_idx,
		const std::vector< numeric_index_type > &	send_list
	)

overridevirtual

将指定范围内的部分值本地化到另一个矢量的相应位置。

参数

first_local_idx	第一个本地索引
last_local_idx	最后一个本地索引
send_list	本地化的部分值的索引列表

在文件 trilinos_epetra_vector.C 第 495 行定义.

参考 libMesh::EpetraVector< T >::close(), libMesh::EpetraVector< T >::init(), libMesh::EpetraVector< T >::localize() , 以及 libMesh::EpetraVector< T >::set().

 {
   // Only good for serial vectors.
   libmesh_assert_equal_to (this->size(), this->local_size());
   libmesh_assert_greater (last_local_idx, first_local_idx);
   libmesh_assert_less_equal (send_list.size(), this->size());
   libmesh_assert_less (last_local_idx, this->size());
 
   const unsigned int my_size       = this->size();
   const unsigned int my_local_size = (last_local_idx - first_local_idx + 1);
 
   // Don't bother for serial cases
   if ((first_local_idx == 0) &&
       (my_local_size == my_size))
     return;
 
   // Build a parallel vector, initialize it with the local
   // parts of (*this)
   EpetraVector<T> parallel_vec(this->comm(), PARALLEL);
 
   parallel_vec.init (my_size, my_local_size, true, PARALLEL);
 
   // Copy part of *this into the parallel_vec
   for (numeric_index_type i=first_local_idx; i<=last_local_idx; i++)
     parallel_vec.set(i,this->el(i));
 
   // localize like normal
   parallel_vec.close();
   parallel_vec.localize (*this, send_list);
   this->close();
 }

template<typename T >

void libMesh::EpetraVector< T >::localize_to_one	(	std::vector< T > &	v_local,
		const processor_id_type	proc_id = `0`
	)		const

overridevirtual

将所有值本地化到一个本地矢量中，仅保留一个处理器上的值。

参数

v_local	存储本地化值的矢量
proc_id	用于本地化的处理器ID

在文件 trilinos_epetra_vector.C 第 555 行定义.

 {
   // This function must be run on all processors at once
   parallel_object_only();
 
   const unsigned int n  = this->size();
   const unsigned int nl = this->local_size();
 
   libmesh_assert_less (pid, this->n_processors());
   libmesh_assert(this->_vec);
 
   v_local.clear();
   v_local.reserve(n);
 
 
   // build up my local part
   for (unsigned int i=0; i<nl; i++)
     v_local.push_back((*this->_vec)[i]);
 
   this->comm().gather (pid, v_local);
 }

template<typename T >

Real libMesh::EpetraVector< T >::max ( ) const

inlineoverridevirtual

获取矢量中的最大值。

返回: 矢量中的最大值

在文件 trilinos_epetra_vector.h 第 1191 行定义.

 {
   libmesh_assert (this->initialized());
 
   T value;
 
   _vec->MaxValue(&value);
 
   return value;
 }

template<typename T >

std::size_t libMesh::EpetraVector< T >::max_allowed_id ( ) const

inlineoverridevirtual

获取允许的最大ID大小。

返回: 最大允许的ID大小

在文件 trilinos_epetra_vector.h 第 1231 行定义.

 {
   // Epetra_Vector seems to use hard-coded ints in its various indexing routines.
   return std::numeric_limits<int>::max();
 }

template<typename T >

Real libMesh::EpetraVector< T >::min ( ) const

inlineoverridevirtual

获取矢量中的最小值。

返回: 矢量中的最小值

在文件 trilinos_epetra_vector.h 第 1176 行定义.

libMesh::ReferenceCounter::_n_objects

 {
   libmesh_assert (this->initialized());
 
   T value;
 
   _vec->MinValue(&value);
 
   return value;
 }

static unsigned int libMesh::ReferenceCounter::n_objects ( )

inlinestaticinherited

Prints the number of outstanding (created, but not yet destroyed) objects.

在文件 reference_counter.h 第 85 行定义.

参考 libMesh::ReferenceCounter::_n_objects.

参考自 libMesh::LibMeshInit::~LibMeshInit().

86 { return _n_objects; }

static Threads::atomic< unsigned int > _n_objects

The number of objects.

Definition: reference_counter.h:132

template<typename T >

T libMesh::EpetraVector< T >::operator() ( const numeric_index_type i ) const

inlineoverridevirtual

获取矢量中索引i处的值。

参数

i 要获取的索引

返回: 索引i处的值

在文件 trilinos_epetra_vector.h 第 1163 行定义.

 {
   libmesh_assert (this->initialized());
   libmesh_assert ( ((i >= this->first_local_index()) &&
                     (i <  this->last_local_index())) );
 
   return (*_vec)[i-this->first_local_index()];
 }

template<typename T >

NumericVector< T > & libMesh::EpetraVector< T >::operator*= ( const NumericVector< T > & v )

overridevirtual

将矢量与另一个矢量逐元素相乘并将结果存储在当前矢量中。

参数

v 另一个矢量

返回: 引用到当前矢量

libMesh::NumericVector::scale

在文件 trilinos_epetra_vector.C 第 131 行定义.

参考 libMesh::closed() , 以及 libMesh::NumericVector< T >::size().

 {
   libmesh_assert(this->closed());
   libmesh_assert_equal_to(size(), v.size());
 
   const EpetraVector<T> & v_vec = cast_ref<const EpetraVector<T> &>(v);
 
   _vec->Multiply(1.0, *v_vec._vec, *_vec, 0.0);
 
   return *this;
 }

template<typename T>

NumericVector<T>& libMesh::NumericVector< T >::operator*= ( const T a )

inlineinherited

将向量缩放为 a ， $\vec{u} \leftarrow a\vec{u}$ 。等价于 u.scale(a)。

参数

a 缩放因子。

返回: 对 *this 的引用。

在文件 numeric_vector.h 第 463 行定义.

参考 libMesh::NumericVector< T >::scale().

463 { this->scale(a); return *this; }

virtual void scale(const T factor)=0

缩放向量的每个元素。

template<typename T >

NumericVector< T > & libMesh::EpetraVector< T >::operator+= ( const NumericVector< T > & v )

overridevirtual

将矢量与另一个矢量相加并将结果存储在当前矢量中。

参数

v 另一个矢量

返回: 引用到当前矢量

在文件 trilinos_epetra_vector.C 第 106 行定义.

参考 libMesh::closed().

 {
   libmesh_assert(this->closed());
 
   this->add(1., v);
 
   return *this;
 }

template<typename T >

NumericVector< T > & libMesh::EpetraVector< T >::operator-= ( const NumericVector< T > & v )

overridevirtual

将矢量与另一个矢量相减并将结果存储在当前矢量中。

参数

v 另一个矢量

返回: 引用到当前矢量

在文件 trilinos_epetra_vector.C 第 119 行定义.

参考 libMesh::closed().

 {
   libmesh_assert(this->closed());
 
   this->add(-1., v);
 
   return *this;
 }

template<typename T >

NumericVector< T > & libMesh::EpetraVector< T >::operator/= ( const NumericVector< T > & v )

overridevirtual

将矢量与另一个矢量逐元素相除并将结果存储在当前矢量中。

参数

v 另一个矢量

返回: 引用到当前矢量

libMesh::NumericVector::scale

在文件 trilinos_epetra_vector.C 第 146 行定义.

参考 libMesh::closed() , 以及 libMesh::NumericVector< T >::size().

 {
   libmesh_assert(this->closed());
   libmesh_assert_equal_to(size(), v.size());
 
   const EpetraVector<T> & v_vec = cast_ref<const EpetraVector<T> &>(v);
 
   _vec->ReciprocalMultiply(1.0, *v_vec._vec, *_vec, 0.0);
 
   return *this;
 }

template<typename T>

NumericVector<T>& libMesh::NumericVector< T >::operator/= ( const T a )

inlineinherited

将向量缩放为 1/a ， $\vec{u} \leftarrow \frac{1}{a}\vec{u}$ 。等价于 u.scale(1.

/a)。

参数

a 缩放因子。

返回: 对 *this 的引用。

在文件 numeric_vector.h 第 481 行定义.

参考 libMesh::NumericVector< T >::scale().

481 { this->scale(1./a); return *this; }

virtual void scale(const T factor)=0

缩放向量的每个元素。

template<typename T>

EpetraVector& libMesh::EpetraVector< T >::operator= ( const EpetraVector< T > & )

delete

template<typename T>

EpetraVector& libMesh::EpetraVector< T >::operator= ( EpetraVector< T > && )

delete

template<typename T >

NumericVector< T > & libMesh::EpetraVector< T >::operator= ( const T s )

overridevirtual

将矢量中的所有元素设置为标量值s。

参数

s 标量值

返回: 引用到当前矢量

在文件 trilinos_epetra_vector.C 第 357 行定义.

 {
   _vec->PutScalar(s_in);
 
   return *this;
 }

template<typename T >

NumericVector< T > & libMesh::EpetraVector< T >::operator= ( const NumericVector< T > & v )

overridevirtual

从其他矢量复制元素。

参数

v 源矢量

返回: 引用到当前矢量

在文件 trilinos_epetra_vector.C 第 368 行定义.

 {
   // This function could be implemented in terms of the copy
   // assignment operator (see other NumericVector subclasses) but that
   // function is currently deleted so calling this function is an error.
   // const EpetraVector<T> * v = cast_ptr<const EpetraVector<T> *>(&v_in);
   // *this = *v;
   libmesh_not_implemented();
   return *this;
 }

template<typename T >

NumericVector< T > & libMesh::EpetraVector< T >::operator= ( const std::vector< T > & v )

overridevirtual

从std::vector复制元素。

参数

v	源std::vector

返回: 引用到当前矢量

Case 1: The vector is the same size of The global vector. Only add the local components.

Case 2: The vector is the same size as our local piece. Insert directly to the local piece.

在文件 trilinos_epetra_vector.C 第 383 行定义.

 {
   T * values = _vec->Values();
 
   if (this->size() == v.size())
     {
       const unsigned int nl=this->local_size();
       const unsigned int fli=this->first_local_index();
 
       for (unsigned int i=0;i<nl;i++)
         values[i]=v[fli+i];
     }
 
   else
     {
       libmesh_assert_equal_to (v.size(), this->local_size());
 
       const unsigned int nl=this->local_size();
 
       for (unsigned int i=0;i<nl;i++)
         values[i]=v[i];
     }
 
   return *this;
 }

template<typename T>

virtual void libMesh::NumericVector< T >::pointwise_divide	(	const NumericVector< T > &	vec1,
		const NumericVector< T > &	vec2
	)

pure virtualinherited

计算该向量与另一个向量的逐点除法。

执行逐点（分量间）除法操作 $u_i \leftarrow \frac{v_{1,i}}{v_{2,i}}$ ，并将结果存储在该向量中。

参数

vec1	被除数的输入向量。
vec2	除数的输入向量。

在 libMesh::PetscVector< T >, libMesh::LaspackVector< T >, libMesh::EigenSparseVector< T > , 以及 libMesh::DistributedVector< T > 内被实现.

template<typename T >

void libMesh::EpetraVector< T >::pointwise_mult	(	const NumericVector< T > &	vec1,
		const NumericVector< T > &	vec2
	)

overridevirtual

计算当前矢量与另一个矢量的逐元素乘积。

参数

vec1	另一个矢量
vec2	另一个矢量

在文件 trilinos_epetra_vector.C 第 345 行定义.

 {
   const EpetraVector<T> * v1 = cast_ptr<const EpetraVector<T> *>(&vec1);
   const EpetraVector<T> * v2 = cast_ptr<const EpetraVector<T> *>(&vec2);
 
   _vec->Multiply(1.0, *v1->_vec, *v2->_vec, 0.0);
 }

template<typename T >

void libMesh::NumericVector< T >::print ( std::ostream & os = libMesh::out ) const

inlinevirtualinherited

打印本地向量的内容，默认输出到 libMesh::out 流。

参数

os	输出流，用于指定输出位置（可选，默认为 libMesh::out）。

在文件 numeric_vector.h 第 1121 行定义.

 {
   libmesh_assert (this->initialized());
   os << "Size\tglobal =  " << this->size()
      << "\t\tlocal =  " << this->local_size() << std::endl;
 
   os << "#\tValue" << std::endl;
   for (auto i : index_range(*this))
     os << i << "\t" << (*this)(i) << std::endl;
 }

template<>

void libMesh::NumericVector< Complex >::print ( std::ostream & os ) const

inlineinherited

在文件 numeric_vector.h 第 1103 行定义.

参考 std::imag(), libMesh::initialized() , 以及 std::real().

 {
   libmesh_assert (this->initialized());
   os << "Size\tglobal =  " << this->size()
      << "\t\tlocal =  " << this->local_size() << std::endl;
 
   // std::complex<>::operator<<() is defined, but use this form
   os << "#\tReal part\t\tImaginary part" << std::endl;
   for (auto i : index_range(*this))
     os << i << "\t"
        << (*this)(i).real() << "\t\t"
        << (*this)(i).imag() << std::endl;
 }

template<typename T >

void libMesh::NumericVector< T >::print_global ( std::ostream & os = libMesh::out ) const

inlinevirtualinherited

打印全局向量的内容，默认输出到 libMesh::out 流。

参数

os	输出流，用于指定输出位置（可选，默认为 libMesh::out）。

在文件 numeric_vector.h 第 1158 行定义.

 {
   libmesh_assert (this->initialized());
 
   std::vector<T> v(this->size());
   this->localize(v);
 
   // 现在我们只需要输出的一个副本
   if (this->processor_id())
     return;
 
   os << "Size\tglobal =  " << this->size() << std::endl;
   os << "#\tValue" << std::endl;
   for (auto i : make_range(v.size()))
     os << i << "\t" << v[i] << std::endl;
 }

template<>

void libMesh::NumericVector< Complex >::print_global ( std::ostream & os ) const

inlineinherited

在文件 numeric_vector.h 第 1136 行定义.

 {
   libmesh_assert (this->initialized());
 
   std::vector<Complex> v(this->size());
   this->localize(v);
 
   // 现在我们只需要输出的一个副本
   if (this->processor_id())
     return;
 
   os << "Size\tglobal =  " << this->size() << std::endl;
   os << "#\tReal part\t\tImaginary part" << std::endl;
   for (auto i : make_range(v.size()))
     os << i << "\t"
        << v[i].real() << "\t\t"
        << v[i].imag() << std::endl;
 }

void libMesh::ReferenceCounter::print_info ( std::ostream & out_stream = libMesh::out )

staticinherited

Prints the reference information, by default to libMesh::out.

在文件 reference_counter.C 第 81 行定义.

参考 libMesh::ReferenceCounter::_enable_print_counter , 以及 libMesh::ReferenceCounter::get_info().

参考自 libMesh::LibMeshInit::~LibMeshInit().

 {
   if (_enable_print_counter)
     out_stream << ReferenceCounter::get_info();
 }

template<typename T>

virtual void libMesh::NumericVector< T >::print_matlab ( const std::string & = "" ) const

inlinevirtualinherited

以Matlab的稀疏矩阵格式打印向量内容。可选择将向量打印到名为 name 的文件中。如果未指定 name，则内容将被打印到屏幕上。

参数

name	用于指定输出文件名的可选参数（默认为空）。

被 libMesh::PetscVector< T > 重载.

在文件 numeric_vector.h 第 859 行定义.

   {
     libmesh_not_implemented();
   }

template<typename T >

bool libMesh::NumericVector< T >::readable ( ) const

inherited

检查该向量是否能够用于全局操作。

返回: 如果向量可以用于全局操作，则返回 true，否则返回 false。

在文件 numeric_vector.C 第 416 行定义.

参考 libMesh::closed() , 以及 libMesh::initialized().

参考自 libMesh::NumericVector< T >::add_vector(), libMesh::NumericVector< T >::compatible() , 以及 libMesh::NumericVector< T >::insert().

 {
   return this->initialized() && this->closed();
 }

template<typename T >

void libMesh::EpetraVector< T >::reciprocal ( )

overridevirtual

计算当前矢量的逐元素倒数。

在文件 trilinos_epetra_vector.C 第 178 行定义.

参考 std::abs().

 {
   // The Epetra::reciprocal() function takes a constant reference to *another* vector,
   // and fills _vec with its reciprocal.  Does that mean we can't pass *_vec as the
   // argument?
   // _vec->reciprocal( *_vec );
 
   // Alternatively, compute the reciprocal by hand... see also the add(T) member that does this...
   const unsigned int nl = _vec->MyLength();
 
   T * values = _vec->Values();
 
   for (unsigned int i=0; i<nl; i++)
     {
       // Don't divide by zero (maybe only check this in debug mode?)
       libmesh_error_msg_if(std::abs(values[i]) < std::numeric_limits<T>::min(),
                            "Error, divide by zero in DistributedVector<T>::reciprocal()!");
 
       values[i] = 1. / values[i];
     }
 
   // Leave the vector in a closed state...
   this->close();
 }

template<typename T >

int libMesh::EpetraVector< T >::ReplaceGlobalValues	(	int	numIDs,
		const int *	GIDs,
		const double *	values
	)

private

将值复制到向量中，覆盖指定索引已存在的任何值。

参数

numIDs	全局 ID 的数量。
GIDs	包含全局 ID 的数组。
values	包含系数值的数组。

返回: 操作的结果代码。

在文件 trilinos_epetra_vector.C 第 628 行定义.

 {
   return( inputValues( numIDs, GIDs, values, false) );
 }

template<typename T >

int libMesh::EpetraVector< T >::ReplaceGlobalValues	(	const Epetra_IntSerialDenseVector &	GIDs,
		const Epetra_SerialDenseVector &	values
	)

private

将值复制到向量中，替换指定 GID 已存在的任何值。

参数

GIDs	包含全局 ID 的数组，必须与伴随的 values 数组的长度相同。
values	包含系数值的数组，必须与伴随的 GIDs 数组的长度相同。

返回: 操作的结果代码。

在文件 trilinos_epetra_vector.C 第 637 行定义.

 {
   if (GIDs.Length() != values.Length()) {
     return(-1);
   }
 
   return( inputValues( GIDs.Length(), GIDs.Values(), values.Values(), false) );
 }

template<typename T >

int libMesh::EpetraVector< T >::ReplaceGlobalValues	(	int	numIDs,
		const int *	GIDs,
		const int *	numValuesPerID,
		const double *	values
	)

private

将值复制到向量中，替换指定索引已存在的任何值。

参数

numIDs	全局 ID 的数量。
GIDs	包含全局 ID 的数组。
numValuesPerID	包含每个全局 ID 对应值的数量的数组。
values	包含系数值的数组。

返回: 操作的结果代码。

在文件 trilinos_epetra_vector.C 第 649 行定义.

 {
   return( inputValues( numIDs, GIDs, numValuesPerID, values, false) );
 }

template<typename T >

void libMesh::EpetraVector< T >::scale ( const T factor )

overridevirtual

缩放矢量的所有元素。

参数

factor 缩放因子

在文件 trilinos_epetra_vector.C 第 319 行定义.

 {
   _vec->Scale(factor_in);
 }

template<typename T >

void libMesh::EpetraVector< T >::set	(	const numeric_index_type	i,
		const T	value
	)

overridevirtual

设置索引i处的值为value。

参数

i	索引
value	要设置的值

在文件 trilinos_epetra_vector.C 第 162 行定义.

参考自 libMesh::EpetraVector< T >::localize().

 {
   int i = static_cast<int> (i_in);
   T value = value_in;
 
   libmesh_assert_less (i_in, this->size());
 
   std::scoped_lock lock(this->_numeric_vector_mutex);
   ReplaceGlobalValues(1, &i, &value);
 
   this->_is_closed = false;
 }

template<typename T>

void libMesh::EpetraVector< T >::setIgnoreNonLocalEntries ( bool flag )

inlineprivate

设置是否应忽略非本地数据值。

参数

flag	如果为 true，则忽略非本地数据值。

在文件 trilinos_epetra_vector.h 第 668 行定义.

参考 libMesh::EpetraVector< T >::ignoreNonLocalEntries_.

   {
     ignoreNonLocalEntries_ = flag;
   }

template<typename T >

numeric_index_type libMesh::EpetraVector< T >::size ( ) const

inlineoverridevirtual

获取矢量的全局大小。

返回: 矢量的全局大小

在文件 trilinos_epetra_vector.h 第 1122 行定义.

参考自 libMesh::EpetraVector< T >::add().

 {
   libmesh_assert (this->initialized());
 
   return _vec->GlobalLength();
 }

template<class T >

Real libMesh::NumericVector< T >::subset_l1_norm ( const std::set< numeric_index_type > & indices ) const

virtualinherited

获取指定条目的向量的 $\ell_1$ -范数，即指定条目的绝对值之和。

参数

indices 指定的索引集合。

返回: 指定条目的向量的 $\ell_1$ -范数。

注解: 索引必须位于此处理器上。

在文件 numeric_vector.C 第 306 行定义.

参考 std::abs(), std::norm() , 以及 libMesh::Real.

 {
   libmesh_assert (this->readable());
 
   const NumericVector<T> & v = *this;
 
   Real norm = 0;
 
   for (const auto & index : indices)
     norm += std::abs(v(index));
 
   this->comm().sum(norm);
 
   return norm;
 }

template<class T >

Real libMesh::NumericVector< T >::subset_l2_norm ( const std::set< numeric_index_type > & indices ) const

virtualinherited

获取指定条目的向量的 $\ell_2$ -范数，即指定条目平方和的平方根。

参数

indices 指定的索引集合。

返回: 指定条目的向量的 $\ell_2$ -范数。

注解: 索引必须位于此处理器上。

在文件 numeric_vector.C 第 323 行定义.

参考 std::norm(), libMesh::TensorTools::norm_sq(), libMesh::Real , 以及 std::sqrt().

 {
   libmesh_assert (this->readable());
 
   const NumericVector<T> & v = *this;
 
   Real norm = 0;
 
   for (const auto & index : indices)
     norm += TensorTools::norm_sq(v(index));
 
   this->comm().sum(norm);
 
   return std::sqrt(norm);
 }

template<class T >

Real libMesh::NumericVector< T >::subset_linfty_norm ( const std::set< numeric_index_type > & indices ) const

virtualinherited

获取指定条目的向量的最大绝对值，即指定条目的 $\ell_{\infty}$ -范数。

参数

indices 指定的索引集合。

返回: 指定条目的向量的最大绝对值。

注解: 索引必须位于此处理器上。

在文件 numeric_vector.C 第 340 行定义.

参考 std::abs(), std::norm() , 以及 libMesh::Real.

 {
   libmesh_assert (this->readable());
 
   const NumericVector<T> & v = *this;
 
   Real norm = 0;
 
   for (const auto & index : indices)
     {
       Real value = std::abs(v(index));
       if (value > norm)
         norm = value;
     }
 
   this->comm().max(norm);
 
   return norm;
 }

template<typename T >

T libMesh::EpetraVector< T >::sum ( ) const

overridevirtual

计算矢量中所有元素的总和。

返回: 所有元素的总和

在文件 trilinos_epetra_vector.C 第 50 行定义.

参考 libMesh::closed().

 {
   libmesh_assert(this->closed());
 
   const unsigned int nl = _vec->MyLength();
 
   T sum=0.0;
 
   T * values = _vec->Values();
 
   for (unsigned int i=0; i<nl; i++)
     sum += values[i];
 
   this->comm().sum(sum);
 
   return sum;
 }

template<typename T >

int libMesh::EpetraVector< T >::SumIntoGlobalValues	(	int	numIDs,
		const int *	GIDs,
		const double *	values
	)

private

将值累积到向量中，将它们添加到指定索引已存在的任何值中。

参数

numIDs	全局 ID 的数量。
GIDs	包含全局 ID 的数组。
values	包含系数值的数组。

返回: 操作的结果代码。

在文件 trilinos_epetra_vector.C 第 597 行定义.

 {
   return( inputValues( numIDs, GIDs, values, true) );
 }

template<typename T >

int libMesh::EpetraVector< T >::SumIntoGlobalValues	(	const Epetra_IntSerialDenseVector &	GIDs,
		const Epetra_SerialDenseVector &	values
	)

private

将值累积到向量中，将它们添加到指定 GID 已存在的任何值中。

参数

GIDs	包含全局 ID 的数组，必须与伴随的 values 数组的长度相同。
values	包含系数值的数组，必须与伴随的 GIDs 数组的长度相同。

返回: 操作的结果代码。

在文件 trilinos_epetra_vector.C 第 606 行定义.

 {
   if (GIDs.Length() != values.Length()) {
     return(-1);
   }
 
   return( inputValues( GIDs.Length(), GIDs.Values(), values.Values(), true) );
 }

template<typename T >

int libMesh::EpetraVector< T >::SumIntoGlobalValues	(	int	numIDs,
		const int *	GIDs,
		const int *	numValuesPerID,
		const double *	values
	)

private

将值累积到向量中，将它们添加到指定索引已存在的任何值中。

参数

numIDs	全局 ID 的数量。
GIDs	包含全局 ID 的数组。
numValuesPerID	包含每个全局 ID 对应值的数量的数组。
values	包含系数值的数组。

返回: 操作的结果代码。

在文件 trilinos_epetra_vector.C 第 618 行定义.

 {
   return( inputValues( numIDs, GIDs, numValuesPerID, values, true) );
 }

template<typename T >

void libMesh::EpetraVector< T >::swap ( NumericVector< T > & v )

inlineoverridevirtual

交换当前矢量与另一个矢量的值。

参数

v 另一个矢量

参考自 libMesh::DofMap::constrain_element_residual(), libMesh::DofMap::enforce_adjoint_constraints_exactly(), libMesh::DofMap::enforce_constraints_exactly(), libMesh::DofMap::enforce_constraints_on_residual(), libMesh::EpetraVector< T >::EpetraVector(), libMesh::DofMap::heterogeneously_constrain_element_jacobian_and_residual(), libMesh::DofMap::heterogeneously_constrain_element_residual(), libMesh::DistributedVector< T >::init(), libMesh::EigenSparseVector< T >::init(), libMesh::EpetraVector< T >::init(), libMesh::LaspackVector< T >::init() , 以及 libMesh::PetscVector< T >::localize().

在文件 trilinos_epetra_vector.h 第 1206 行定义.

 {
   NumericVector<T>::swap(other);
 
   EpetraVector<T> & v = cast_ref<EpetraVector<T> &>(other);
 
   std::swap(_vec, v._vec);
   _map.swap(v._map);
   std::swap(_destroy_vec_on_exit, v._destroy_vec_on_exit);
   std::swap(myFirstID_, v.myFirstID_);
   std::swap(myNumIDs_, v.myNumIDs_);
   std::swap(myCoefs_, v.myCoefs_);
   std::swap(nonlocalIDs_, v.nonlocalIDs_);
   std::swap(nonlocalElementSize_, v.nonlocalElementSize_);
   std::swap(numNonlocalIDs_, v.numNonlocalIDs_);
   std::swap(allocatedNonlocalLength_, v.allocatedNonlocalLength_);
   std::swap(nonlocalCoefs_, v.nonlocalCoefs_);
   std::swap(last_edit, v.last_edit);
   std::swap(ignoreNonLocalEntries_, v.ignoreNonLocalEntries_);
 }

template<typename T>

ParallelType libMesh::NumericVector< T >::type ( ) const

inlineinherited

获取向量的类型。

返回: 向量的类型（SERIAL、PARALLEL、GHOSTED）。

在文件 numeric_vector.h 第 166 行定义.

参考 libMesh::NumericVector< T >::_type.

166 { return _type; }

libMesh::NumericVector::_type

ParallelType _type

向量的类型。

Definition: numeric_vector.h:923

template<typename T>

ParallelType& libMesh::NumericVector< T >::type ( )

inlineinherited

获取向量的类型。

返回: 向量的类型（SERIAL、PARALLEL、GHOSTED）。

在文件 numeric_vector.h 第 173 行定义.

参考 libMesh::NumericVector< T >::_type.

173 { return _type; }

libMesh::NumericVector::_type

ParallelType _type

向量的类型。

Definition: numeric_vector.h:923

template<typename T>

Epetra_Vector* libMesh::EpetraVector< T >::vec ( )

inline

返回原始的 Epetra_Vector 指针。

注解: 通常情况下，用户级别的代码不需要使用这个函数。; 请不要做一些疯狂的事情，比如删除指针，否则可能会发生非常糟糕的事情！

返回: 指向 Epetra_Vector 对象的指针。

在文件 trilinos_epetra_vector.h 第 553 行定义.

libMesh::EpetraVector::_vec

参考自 libMesh::EpetraMatrix< T >::get_diagonal().

553 { libmesh_assert(_vec); return _vec; }

Epetra_Vector * _vec

用于保存向量条目的实际 Epetra 向量数据类型。

Definition: trilinos_epetra_vector.h:560

template<typename T >

void libMesh::EpetraVector< T >::zero ( )

inlineoverridevirtual

将所有元素置零。

在文件 trilinos_epetra_vector.h 第 1087 行定义.

参考 libMesh::closed() , 以及 libMesh::initialized().

 {
   libmesh_assert (this->initialized());
   libmesh_assert (this->closed());
 
   _vec->PutScalar(0.0);
 }

template<typename T >

std::unique_ptr< NumericVector< T > > libMesh::EpetraVector< T >::zero_clone ( ) const

inlineoverridevirtual

创建零克隆矢量。

返回: 零克隆矢量的唯一指针