On 1/8/21 12:15 PM, Lixing Zhu wrote:
This way I can freely use bubble functions and only construct local K and F of
FE_1. However, one issue is that, apparently, triangulation.cell does not
support cell->get_dof_indices(local_dof_indices). Is there any way to call the
DOFhandler by the cell in triangulation and then access the local_dof_indices?
You want a statement like this (from step-43):
auto cell = darcy_dof_handler.begin_active();
const auto endc = darcy_dof_handler.end();
auto saturation_cell = saturation_dof_handler.begin_active();
for (; cell != endc; ++cell, ++saturation_cell)
Here, we're walking iterators over two DoFHandler objects in parallel. In your
case, one will be the DoFHandler iterator and the other one the Triangulation
iterator. Or, you could do
fe_values_b.reinit
(static_cast<typename Triangulation<dim>::active_cell_iterator>
(cell));
to cast the DoFHandler iterator to a Triangulation iterator.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/70d80b6d-d7f5-df0b-1262-9758a7dd4dd9%40colostate.edu.
