108 patch.reference_cell = cell->face(face_number)->reference_cell();
114 for (
const unsigned int vertex : cell->face(face_number)->vertex_indices())
118 cell->reference_cell().face_to_cell_vertices(
119 face_number, vertex, cell->combined_face_orientation(face_number)));
122 data.mapping_collection[0].transform_unit_to_real_cell(
129 if (
data.n_datasets > 0)
131 data.reinit_all_fe_values(this->dof_data, cell, face_number);
138 const std::vector<Point<dim>> &q_points =
142 patch.data.reinit(
data.n_datasets + dim, q_points.size());
145 patch.points_are_available =
true;
147 for (
unsigned int i = 0; i < dim; ++i)
148 for (
unsigned int q = 0;
q < n_q_points; ++
q)
149 patch.data(
patch.data.size(0) - dim + i,
q) = q_points[
q][i];
152 unsigned int offset = 0;
159 const unsigned int n_components =
162 this->dof_data[
dataset]->postprocessor;
163 if (postprocessor !=
nullptr)
168 postprocessor->get_needed_update_flags();
170 if (n_components == 1)
175 this->dof_data[
dataset]->get_function_values(
179 data.patch_values_scalar.solution_values);
181 this->dof_data[
dataset]->get_function_gradients(
185 data.patch_values_scalar.solution_gradients);
187 this->dof_data[
dataset]->get_function_hessians(
191 data.patch_values_scalar.solution_hessians);
194 data.patch_values_scalar.evaluation_points =
198 data.patch_values_scalar.normals =
202 dh_cell(&cell->get_triangulation(),
209 postprocessor->evaluate_scalar_field(
210 data.patch_values_scalar,
217 data.resize_system_vectors(n_components);
219 this->dof_data[
dataset]->get_function_values(
223 data.patch_values_system.solution_values);
225 this->dof_data[
dataset]->get_function_gradients(
229 data.patch_values_system.solution_gradients);
231 this->dof_data[
dataset]->get_function_hessians(
235 data.patch_values_system.solution_hessians);
238 data.patch_values_system.evaluation_points =
242 data.patch_values_system.normals =
246 dh_cell(&cell->get_triangulation(),
253 postprocessor->evaluate_vector_field(
254 data.patch_values_system,
258 for (
unsigned int q = 0;
q < n_q_points; ++
q)
259 for (
unsigned int component = 0;
260 component < this->dof_data[
dataset]->n_output_variables;
262 patch.data(offset + component,
q) =
269 if (n_components == 1)
271 this->dof_data[
dataset]->get_function_values(
275 data.patch_values_scalar.solution_values);
276 for (
unsigned int q = 0;
q < n_q_points; ++
q)
278 data.patch_values_scalar.solution_values[
q];
282 data.resize_system_vectors(n_components);
283 this->dof_data[
dataset]->get_function_values(
287 data.patch_values_system.solution_values);
288 for (
unsigned int component = 0; component < n_components;
290 for (
unsigned int q = 0;
q < n_q_points; ++
q)
291 patch.data(offset + component,
q) =
292 data.patch_values_system.solution_values[
q](component);
295 offset += this->dof_data[
dataset]->n_output_variables;
308 "The current function is trying to generate cell-data output "
309 "for a face that does not belong to an active cell. This is "
315 const double value = this->cell_data[
dataset]->get_cell_data_value(
318 for (
unsigned int q = 0;
q < n_q_points; ++
q)
349 const unsigned int n_subdivisions =
352 Assert(n_subdivisions >= 1,
359 this->validate_dataset_names();
361 unsigned int n_datasets = this->cell_data.size();
362 for (
unsigned int i = 0; i < this->dof_data.size(); ++i)
363 n_datasets += this->dof_data[i]->n_output_variables;
373 face = next_face(face))
377 this->patches.
clear();
384 if (this->dof_data[
dataset]->postprocessor)
386 this->dof_data[
dataset]->n_output_variables;
391 for (
unsigned int i = 0; i < this->dof_data.size(); ++i)
392 if (this->dof_data[i]->postprocessor)
394 this->dof_data[i]->postprocessor->get_needed_update_flags();
415 this->build_one_patch(cell_and_face, data, patch);
418 internal::DataOutFacesImplementation::append_patch_to_list<dim, spacedim>(
419 patch, this->patches);
void build_one_patch(const FaceDescriptor *cell_and_face, internal::DataOutFacesImplementation::ParallelData< dim, spacedim > &data, DataOutBase::Patch< patch_dim, patch_spacedim > &patch)
void run(const std::vector< std::vector< Iterator > > &colored_iterators, Worker worker, Copier copier, const ScratchData &sample_scratch_data, const CopyData &sample_copy_data, const unsigned int queue_length=2 *MultithreadInfo::n_threads(), const unsigned int chunk_size=8)
void append_patch_to_list(const DataOutBase::Patch< DataOutFaces< dim, spacedim >::patch_dim, DataOutFaces< dim, spacedim >::patch_spacedim > &patch, std::vector< DataOutBase::Patch< DataOutFaces< dim, spacedim >::patch_dim, DataOutFaces< dim, spacedim >::patch_spacedim > > &patches)
ParallelData(const unsigned int n_datasets, const unsigned int n_subdivisions, const std::vector< unsigned int > &n_postprocessor_outputs, const Mapping< dim, spacedim > &mapping, const std::vector< std::shared_ptr<::hp::FECollection< dim, spacedim > > > &finite_elements, const UpdateFlags update_flags)