16 #ifndef dealii_trilinos_epetra_vector_h 17 #define dealii_trilinos_epetra_vector_h 20 #include <deal.II/base/config.h> 22 #if defined(DEAL_II_WITH_TRILINOS) && defined(DEAL_II_WITH_MPI) 24 # include <deal.II/base/index_set.h> 25 # include <deal.II/base/subscriptor.h> 27 # include <deal.II/lac/trilinos_epetra_communication_pattern.h> 28 # include <deal.II/lac/vector_operation.h> 29 # include <deal.II/lac/vector_space_vector.h> 30 # include <deal.II/lac/vector_type_traits.h> 32 # include <Epetra_FEVector.h> 37 DEAL_II_NAMESPACE_OPEN
42 template <
typename Number>
43 class ReadWriteVector;
45 namespace EpetraWrappers
78 const MPI_Comm &communicator);
88 const MPI_Comm &communicator,
89 const bool omit_zeroing_entries =
false);
97 const bool omit_zeroing_entries =
false)
override;
127 std::shared_ptr<const CommunicationPatternBase> communication_pattern =
128 std::shared_ptr<const CommunicationPatternBase>())
override;
165 add(
const double a)
override;
279 size()
const override;
305 const Epetra_FEVector &
319 print(std::ostream & out,
320 const unsigned int precision = 3,
321 const bool scientific =
true,
322 const bool across =
true)
const override;
350 <<
"An error with error number " << arg1
351 <<
" occurred while calling a Trilinos function");
361 const MPI_Comm &mpi_comm);
397 DEAL_II_NAMESPACE_CLOSE
virtual double l1_norm() const override
static::ExceptionBase & ExcDifferentParallelPartitioning()
virtual bool all_zero() const override
const Epetra_FEVector & trilinos_vector() const
virtual ::IndexSet locally_owned_elements() const override
virtual double l2_norm() const override
virtual size_type size() const override
virtual void scale(const VectorSpaceVector< double > &scaling_factors) override
virtual std::size_t memory_consumption() const override
std::unique_ptr< Epetra_FEVector > vector
virtual void equ(const double a, const VectorSpaceVector< double > &V) override
#define DeclException1(Exception1, type1, outsequence)
void create_epetra_comm_pattern(const IndexSet &source_index_set, const MPI_Comm &mpi_comm)
virtual Vector & operator*=(const double factor) override
bool has_ghost_elements() const
#define DeclException0(Exception0)
Vector & operator=(const Vector &V)
static::ExceptionBase & ExcTrilinosError(int arg1)
static::ExceptionBase & ExcVectorTypeNotCompatible()
void reinit(const IndexSet ¶llel_partitioner, const MPI_Comm &communicator, const bool omit_zeroing_entries=false)
std::shared_ptr< const CommunicationPattern > epetra_comm_pattern
virtual void sadd(const double s, const double a, const VectorSpaceVector< double > &V) override
virtual Vector & operator+=(const VectorSpaceVector< double > &V) override
virtual Vector & operator/=(const double factor) override
virtual void print(std::ostream &out, const unsigned int precision=3, const bool scientific=true, const bool across=true) const override
virtual double mean_value() const override
MPI_Comm get_mpi_communicator() const
virtual Vector & operator-=(const VectorSpaceVector< double > &V) override
virtual double operator*(const VectorSpaceVector< double > &V) const override
virtual void add(const double a) override
virtual double linfty_norm() const override
virtual double add_and_dot(const double a, const VectorSpaceVector< double > &V, const VectorSpaceVector< double > &W) override
::IndexSet source_stored_elements