Hi Wolfgang, Many thanks! The issue is resolved and the code is working.
Regards, Lixing On Saturday, January 9, 2021 at 4:07:27 AM UTC+8 Wolfgang Bangerth wrote: > 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/d85fad4e-baa7-4984-bde2-b1f1f5093bd0n%40googlegroups.com.
