42template <
typename Number>
46 <<
"Called set_dof_values_by_interpolation(), but"
47 <<
" the element to be set, value " << std::setprecision(16)
48 << arg1 <<
", does not match with the non-zero value "
49 << std::setprecision(16) << arg2 <<
" already set before.");
59 template <
typename Number>
60 std::enable_if_t<!std::is_unsigned_v<Number>,
67 template <
typename Number>
68 std::enable_if_t<std::is_unsigned_v<Number>, Number>
81 decltype(std::declval<const T>().set_ghost_state(std::declval<bool>()));
85 is_supported_operation<set_ghost_state_t, T>;
89 std::enable_if_t<has_set_ghost_state<VectorType>, VectorType> * =
nullptr>
93 vector.set_ghost_state(ghosted);
98 std::enable_if_t<!has_set_ghost_state<VectorType>, VectorType> * =
nullptr>
119 OutputVector &values,
120 const bool perform_check)
126 using VectorNumber =
typename OutputVector::value_type;
127 constexpr bool is_dealii_vector =
128 std::is_same_v<OutputVector, Vector<VectorNumber>> ||
129 std::is_same_v<OutputVector, BlockVector<VectorNumber>> ||
130 std::is_same_v<OutputVector,
132 std::is_same_v<OutputVector,
134 if (perform_check && is_dealii_vector)
136 const bool old_ghost_state = values.has_ghost_elements();
142 for (
unsigned int i = 0; i < cell.
get_fe().n_dofs_per_cell(); ++i)
147 local_values_old[i] == number() ||
148 get_abs(local_values_old[i] - local_values[i]) <=
149 get_abs(local_values_old[i] + local_values[i]) * 100000. *
151 number>::real_type>::epsilon(),
152 ExcNonMatchingElementsSetDofValuesByInterpolation<number>(
153 local_values[i], local_values_old[i]));
173 OutputVector &values,
177 OutputVector &values)> &processor)
185 if (cell.is_active() && !cell.is_artificial())
194 processor(cell, local_values, values);
199 ExcMessage(
"Incorrect size of local_values vector."));
202 cell.
get_fe().n_dofs_per_cell(),
205 cell.
get_fe().get_interpolation_matrix(
212 if ((tmp.
size() > 0) && (local_values.
size() > 0))
213 interpolation.
vmult(tmp, local_values);
216 processor(cell, tmp, values);
225 "You cannot call this function on non-active cells "
226 "of DoFHandler objects unless you provide an explicit "
227 "finite element index because they do not have naturally "
228 "associated finite element spaces associated: degrees "
229 "of freedom are only distributed on active cells for which "
230 "the active FE index has been set."));
238 ExcVectorDoesNotMatch()));
241 ExcVectorDoesNotMatch()));
245 for (
unsigned int child = 0; child < cell.n_children(); ++child)
249 .
vmult(tmp, local_values);
251 *cell.
child(child), tmp, values, fe_index, processor);
260template <
int dim,
int spacedim,
bool lda>
261template <
class OutputVector,
typename number>
265 OutputVector &values,
267 const bool perform_check)
const
269 internal::process_by_interpolation<dim, spacedim, lda, OutputVector, number>(
276 OutputVector &values) {
282template <
int dim,
int spacedim,
bool lda>
283template <
class OutputVector,
typename number>
288 OutputVector &values,
291 internal::process_by_interpolation<dim, spacedim, lda, OutputVector, number>(
298 OutputVector &values) {
299 std::vector<types::global_dof_index> dof_indices(
300 cell.
get_fe().n_dofs_per_cell());
311#include "dofs/dof_accessor_set.inst"
void distribute_local_to_global(const InVector &local_vector, const std::vector< size_type > &local_dof_indices, OutVector &global_vector) const
const DoFHandler< dim, spacedim > & get_dof_handler() const
void get_dof_values(const InputVector &values, Vector< number > &local_values) const
const FiniteElement< dimension_, space_dimension_ > & get_fe() const
TriaIterator< DoFCellAccessor< dimension_, space_dimension_, level_dof_access > > child(const unsigned int i) const
void distribute_local_to_global_by_interpolation(const Vector< number > &local_values, OutputVector &values, const types::fe_index fe_index=numbers::invalid_fe_index) const
void set_dof_values(const Vector< number > &local_values, OutputVector &values) const
void set_dof_values_by_interpolation(const Vector< number > &local_values, OutputVector &values, const types::fe_index fe_index=numbers::invalid_fe_index, const bool perform_check=false) const
void get_dof_indices(std::vector< types::global_dof_index > &dof_indices) const
types::fe_index active_fe_index() const
const FiniteElement< dim, spacedim > & get_fe(const types::fe_index index=0) const
bool has_hp_capabilities() const
types::global_dof_index n_dofs() const
unsigned int n_dofs_per_cell() const
virtual const FullMatrix< double > & get_prolongation_matrix(const unsigned int child, const RefinementCase< dim > &refinement_case=RefinementCase< dim >::isotropic_refinement) const
void vmult(Vector< number2 > &w, const Vector< number2 > &v, const bool adding=false) const
virtual size_type size() const override
#define DEAL_II_NAMESPACE_OPEN
constexpr bool running_in_debug_mode()
#define DEAL_II_NAMESPACE_CLOSE
#define Assert(cond, exc)
#define DeclException2(Exception2, type1, type2, outsequence)
static ::ExceptionBase & ExcNonMatchingElementsSetDofValuesByInterpolation(Number arg1, Number arg2)
static ::ExceptionBase & ExcMessage(std::string arg1)
decltype(std::declval< const T >().set_ghost_state(std::declval< bool >())) set_ghost_state_t
constexpr bool has_set_ghost_state
std::enable_if_t<!std::is_unsigned_v< Number >, typename numbers::NumberTraits< Number >::real_type > get_abs(const Number a)
void set_ghost_state(VectorType &vector, const bool ghosted)
void set_dof_values(const DoFCellAccessor< dim, spacedim, lda > &cell, const Vector< number > &local_values, OutputVector &values, const bool perform_check)
void process_by_interpolation(const DoFCellAccessor< dim, spacedim, lda > &cell, const Vector< number > &local_values, OutputVector &values, const types::fe_index fe_index_, const std::function< void(const DoFCellAccessor< dim, spacedim, lda > &cell, const Vector< number > &local_values, OutputVector &values)> &processor)
constexpr types::fe_index invalid_fe_index
::VectorizedArray< Number, width > abs(const ::VectorizedArray< Number, width > &)
unsigned short int fe_index