198 const std::function<
void(DOFINFO &,
typename INFOBOX::CellInfo &)>
200 const std::function<
void(DOFINFO &,
typename INFOBOX::CellInfo &)>
202 const std::function<
void(DOFINFO &,
204 typename INFOBOX::CellInfo &,
205 typename INFOBOX::CellInfo &)> &face_worker,
208 const bool ignore_subdomain =
209 (cell->get_triangulation().locally_owned_subdomain() ==
213 cell->level_subdomain_id() :
214 cell->subdomain_id();
216 const bool own_cell =
218 (csid == cell->get_triangulation().locally_owned_subdomain());
225 dof_info.
cell.reinit(cell);
228 const bool integrate_cell = (cell_worker !=
nullptr);
229 const bool integrate_boundary = (boundary_worker !=
nullptr);
230 const bool integrate_interior_face = (face_worker !=
nullptr);
233 info.cell.reinit(dof_info.
cell);
240 cell_worker(dof_info.
cell, info.cell);
246 info.post_cell(dof_info);
248 if (integrate_interior_face || integrate_boundary)
249 for (
const unsigned int face_no : cell->face_indices())
251 typename IteratorType::AccessorType::Container::face_iterator face =
253 if (cell->at_boundary(face_no) &&
254 !cell->has_periodic_neighbor(face_no))
257 if (integrate_boundary && own_cell)
260 dof_info.
interior[face_no].reinit(cell, face, face_no);
261 info.boundary.reinit(dof_info.
interior[face_no]);
262 boundary_worker(dof_info.
interior[face_no], info.boundary);
265 else if (integrate_interior_face)
269 cell->neighbor_or_periodic_neighbor(face_no);
272 if (neighbor->is_level_cell())
273 neighbid = neighbor->level_subdomain_id();
275 else if (neighbor->is_active())
276 neighbid = neighbor->subdomain_id();
278 const bool own_neighbor =
281 cell->get_triangulation().locally_owned_subdomain());
284 if (!own_cell && !own_neighbor)
288 if (own_cell && own_neighbor &&
293 if (own_cell != own_neighbor &&
300 const bool periodic_neighbor =
301 cell->has_periodic_neighbor(face_no);
303 if ((!periodic_neighbor && cell->neighbor_is_coarser(face_no)) ||
304 (periodic_neighbor &&
305 cell->periodic_neighbor_is_coarser(face_no)))
316 const std::pair<unsigned int, unsigned int> neighbor_face_no =
318 cell->periodic_neighbor_of_coarser_periodic_neighbor(
320 cell->neighbor_of_coarser_neighbor(face_no);
321 const typename IteratorType::AccessorType::Container::
322 face_iterator nface =
323 neighbor->face(neighbor_face_no.first);
327 dof_info.
interior[face_no].reinit(cell, face, face_no);
328 info.face.reinit(dof_info.
interior[face_no]);
329 dof_info.
exterior[face_no].reinit(neighbor,
331 neighbor_face_no.first,
332 neighbor_face_no.second);
333 info.subface.reinit(dof_info.
exterior[face_no]);
335 face_worker(dof_info.
interior[face_no],
345 neighbor->has_children())
350 "Assembling from both sides for own_faces is not "
351 "supported with hanging nodes!"));
356 Assert(cell->level() == neighbor->level(),
362 if (own_cell && own_neighbor &&
376 if (own_cell && !own_neighbor &&
381 const unsigned int neighbor_face_no =
383 cell->periodic_neighbor_face_no(face_no) :
384 cell->neighbor_face_no(face_no);
385 Assert(periodic_neighbor ||
386 neighbor->face(neighbor_face_no) == face,
391 dof_info.
interior[face_no].reinit(cell, face, face_no);
392 info.face.reinit(dof_info.
interior[face_no]);
393 dof_info.
exterior[face_no].reinit(neighbor,
397 info.neighbor.reinit(dof_info.
exterior[face_no]);
399 face_worker(dof_info.
interior[face_no],
410 info.post_faces(dof_info);
417 cell_worker(dof_info.
cell, info.cell);
446 const std::function<
void(DOFINFO &,
typename INFOBOX::CellInfo &)>
448 const std::function<
void(DOFINFO &,
typename INFOBOX::CellInfo &)>
450 const std::function<
void(DOFINFO &,
452 typename INFOBOX::CellInfo &,
453 typename INFOBOX::CellInfo &)> &face_worker,
454 AssemblerType &assembler,
459 assembler.initialize_info(dof_info.
cell,
false);
462 assembler.initialize_info(dof_info.
interior[i],
true);
463 assembler.initialize_info(dof_info.
exterior[i],
true);
470 [&cell_worker, &boundary_worker, &face_worker, &lctrl](
472 cell_action<INFOBOX, DOFINFO, dim, spacedim, IteratorType>(
482 ::internal::assemble<dim, DOFINFO, AssemblerType>(dinfo,
virtual void face(DoFInfo< dim, spacedim, number > &dinfo1, DoFInfo< dim, spacedim, number > &dinfo2, IntegrationInfo< dim, spacedim > &info1, IntegrationInfo< dim, spacedim > &info2) const
void cell_action(IteratorType cell, DoFInfoBox< dim, DOFINFO > &dof_info, INFOBOX &info, const std::function< void(DOFINFO &, typename INFOBOX::CellInfo &)> &cell_worker, const std::function< void(DOFINFO &, typename INFOBOX::CellInfo &)> &boundary_worker, const std::function< void(DOFINFO &, DOFINFO &, typename INFOBOX::CellInfo &, typename INFOBOX::CellInfo &)> &face_worker, const LoopControl &loop_control)
void integration_loop(IteratorType begin, std_cxx20::type_identity_t< IteratorType > end, DoFInfo< dim, spacedim > &dof_info, IntegrationInfoBox< dim, spacedim > &box, const LocalIntegrator< dim, spacedim > &integrator, AssemblerType &assembler, const LoopControl &lctrl=LoopControl())
void loop(IteratorType begin, std_cxx20::type_identity_t< IteratorType > end, DOFINFO &dinfo, INFOBOX &info, const std::function< void(DOFINFO &, typename INFOBOX::CellInfo &)> &cell_worker, const std::function< void(DOFINFO &, typename INFOBOX::CellInfo &)> &boundary_worker, const std::function< void(DOFINFO &, DOFINFO &, typename INFOBOX::CellInfo &, typename INFOBOX::CellInfo &)> &face_worker, AssemblerType &assembler, const LoopControl &lctrl=LoopControl())
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)