16 #ifndef dealii_vector_element_access_h 17 #define dealii_vector_element_access_h 20 #include <deal.II/lac/trilinos_epetra_vector.h> 22 DEAL_II_NAMESPACE_OPEN
26 template <
typename VectorType>
31 add(
const typename VectorType::value_type value,
36 set(
typename VectorType::value_type value,
40 static typename VectorType::value_type
46 template <
typename VectorType>
48 ElementAccess<VectorType>::add(
const typename VectorType::value_type value,
57 template <
typename VectorType>
59 ElementAccess<VectorType>::set(
const typename VectorType::value_type value,
68 template <
typename VectorType>
69 inline typename VectorType::value_type
70 ElementAccess<VectorType>::get(
const VectorType & V,
78 #if defined(DEAL_II_WITH_TRILINOS) && defined(DEAL_II_WITH_MPI) 81 ElementAccess<LinearAlgebra::EpetraWrappers::Vector>::add(
88 TrilinosWrappers::types::int_type trilinos_i =
89 vector.Map().LID(static_cast<TrilinosWrappers::types::int_type>(i));
91 vector[0][trilinos_i] += value;
98 ElementAccess<LinearAlgebra::EpetraWrappers::Vector>::set(
105 TrilinosWrappers::types::int_type trilinos_i =
106 vector.Map().LID(static_cast<TrilinosWrappers::types::int_type>(i));
108 vector[0][trilinos_i] = value;
114 ElementAccess<LinearAlgebra::EpetraWrappers::Vector>::get(
120 TrilinosWrappers::types::int_type trilinos_i =
121 vector.Map().LID(static_cast<TrilinosWrappers::types::int_type>(i));
123 return vector[0][trilinos_i];
128 DEAL_II_NAMESPACE_CLOSE
const Epetra_FEVector & trilinos_vector() const
unsigned long long int global_dof_index