38template <
int dim,
int spacedim>
48template <
int dim,
int spacedim>
60 for (
typename std::vector<DataVector>::const_iterator i =
dof_data.begin();
64 for (
typename std::vector<DataVector>::const_iterator i =
cell_data.begin();
71template <
int dim,
int spacedim>
80template <
int dim,
int spacedim>
85 std::vector<std::string> names;
86 names.push_back(name);
91template <
int dim,
int spacedim>
94 const std::vector<std::string> &names,
107 for (
const auto &name : names)
134template <
int dim,
int spacedim>
135template <
typename number>
138 const std::string &name)
140 const unsigned int n_components =
dof_handler->get_fe(0).n_components();
142 std::vector<std::string> names;
146 if ((n_components == 1) ||
147 (
vec.size() ==
dof_handler->get_triangulation().n_active_cells()))
149 names.resize(1, name);
155 names.resize(n_components);
156 for (
unsigned int i = 0; i < n_components; ++i)
160 names[i] = name +
namebuf.str();
168template <
int dim,
int spacedim>
169template <
typename number>
173 const std::vector<std::string> &names)
180 (names.size() == 1)) ||
182 (names.size() ==
dof_handler->get_fe(0).n_components())),
184 names.size(),
dof_handler->get_fe(0).n_components()));
185 for (
const auto &name : names)
188 Assert(name.find_first_not_of(
"abcdefghijklmnopqrstuvwxyz"
189 "ABCDEFGHIJKLMNOPQRSTUVWXYZ"
190 "0123456789_<>()") == std::string::npos,
193 name.find_first_not_of(
"abcdefghijklmnopqrstuvwxyz"
194 "ABCDEFGHIJKLMNOPQRSTUVWXYZ"
195 "0123456789_<>()")));
200 typename std::vector<DataVector>::iterator data_vector =
dof_data.begin();
201 for (; data_vector !=
dof_data.end(); ++data_vector)
202 if (data_vector->names == names)
204 data_vector->data.reinit(
vec.size());
205 std::copy(
vec.begin(),
vec.end(), data_vector->data.begin());
222 dof_handler->get_triangulation().n_active_cells()))
224 typename std::vector<DataVector>::iterator data_vector =
226 for (; data_vector !=
cell_data.end(); ++data_vector)
227 if (data_vector->names == names)
229 data_vector->data.reinit(
vec.size());
230 std::copy(
vec.begin(),
vec.end(), data_vector->data.begin());
243template <
int dim,
int spacedim>
249 unsigned int n_subdivisions =
252 Assert(n_subdivisions >= 1,
260 const unsigned int n_components =
dof_handler->get_fe(0).n_components();
261 const unsigned int n_datasets =
300 std::vector<Vector<double>> patch_values_system(n_q_points,
310 default_patch.reference_cell = ReferenceCells::get_hypercube<dim + 1>();
311 default_patch.data.reinit(n_datasets, n_q_points * (n_subdivisions + 1));
316 typename std::vector<
392 if (n_components == 1)
396 for (
unsigned int i = 0; i < n_subdivisions + 1; ++i)
397 for (
unsigned int q = 0;
q < n_q_points; ++
q)
405 patch_values_system);
406 for (
unsigned int component = 0; component < n_components;
408 for (
unsigned int i = 0; i < n_subdivisions + 1; ++i)
409 for (
unsigned int q = 0;
q < n_q_points; ++
q)
411 q + n_q_points * i) =
412 patch_values_system[
q](component);
420 for (
unsigned int q = 0;
q < n_q_points; ++
q)
421 for (
unsigned int i = 0; i < n_subdivisions + 1; ++i)
423 q * (n_subdivisions + 1) + i) = value;
430template <
int dim,
int spacedim>
436 for (
typename std::vector<DataVector>::iterator i =
dof_data.begin();
441 for (
typename std::vector<DataVector>::iterator i =
cell_data.begin();
449template <
int dim,
int spacedim>
464template <
int dim,
int spacedim>
475template <
int dim,
int spacedim>
476std::vector<std::string>
479 std::vector<std::string> names;
480 for (
typename std::vector<DataVector>::const_iterator
dataset =
484 names.insert(names.end(),
dataset->names.begin(),
dataset->names.end());
485 for (
typename std::vector<DataVector>::const_iterator
dataset =
489 names.insert(names.end(),
dataset->names.begin(),
dataset->names.end());
497#include "numerics/data_out_stack.inst"
unsigned int default_subdivisions
void validate_dataset_names() const
std::size_t memory_consumption() const
ObserverPointer< const DoFHandler< dim, spacedim >, DataOutStack< dim, spacedim > > dof_handler
void declare_data_vector(const std::string &name, const VectorType vector_type)
virtual std::vector< std::string > get_dataset_names() const override
std::vector< DataVector > cell_data
void new_parameter_value(const double parameter_value, const double parameter_step)
std::vector< DataVector > dof_data
void finish_parameter_value()
void attach_dof_handler(const DoFHandler< dim, spacedim > &dof_handler)
void build_patches(const unsigned int n_subdivisions=0)
std::vector<::DataOutBase::Patch< patch_dim, patch_spacedim > > patches
virtual const std::vector< ::DataOutBase::Patch< DataOutStack< dim, spacedim >::patch_dim, DataOutStack< dim, spacedim >::patch_spacedim > > & get_patches() const override
void add_data_vector(const Vector< number > &vec, const std::string &name)
#define DEAL_II_NAMESPACE_OPEN
constexpr bool running_in_debug_mode()
#define DEAL_II_NAMESPACE_CLOSE
#define DEAL_II_ASSERT_UNREACHABLE()
#define DEAL_II_NOT_IMPLEMENTED()
static ::ExceptionBase & ExcInvalidNumberOfSubdivisions(int arg1)
static ::ExceptionBase & ExcDataNotCleared()
static ::ExceptionBase & ExcDataAlreadyAdded()
static ::ExceptionBase & ExcVectorNotDeclared(std::string arg1)
static ::ExceptionBase & ExcNotImplemented()
static ::ExceptionBase & ExcInvalidCharacter(std::string arg1, std::size_t arg2)
static ::ExceptionBase & ExcInvalidNumberOfNames(int arg1, int arg2)
#define Assert(cond, exc)
static ::ExceptionBase & ExcNoDoFHandlerSelected()
static ::ExceptionBase & ExcNameAlreadyUsed(std::string arg1)
static ::ExceptionBase & ExcInternalError()
typename ActiveSelector::active_cell_iterator active_cell_iterator
@ update_values
Shape function values.
std::vector< index_type > data
std::enable_if_t< std::is_fundamental_v< T >, std::size_t > memory_consumption(const T &t)
std::size_t memory_consumption() const
std::vector< std::string > names