For the record, the lines of code in my last post were wrong. It should be
for the boundary worker :
// new cell iterator for other dof_handler
typename DoFHandler<dim>::active_cell_iterator cell_2 (&cell->
get_triangulation(), cell->level(), cell->index(), &dof_handler_HE);
// Initialise FEValues with FACE quadrature formula and reinitialise
with current cell
FEFaceValues<dim> fe_values_HE_bound(mappingQ, fe_HE,
quadrature_face_formula, update_values);
fe_values_HE_bound.reinit(cell_2, face_no);
std::vector<double> T_old_data_bound(q_points.size());
fe_values_HE_bound.get_function_values(old_T,T_old_data_bound);
Best,
Sylvain M.
El viernes, 8 de octubre de 2021 a la(s) 14:30:54 UTC+2, Wolfgang Bangerth
escribió:
> On 10/8/21 5:13 AM, Sylvain Mathonnière wrote:
> > How can one suspect there is such problem with the code in those cases ?
> Just
> > with experience I guess and doing multiple simulations ?
>
> Making every mistake myself at least once in the last 20 years has
> certainly
> helped me have an eye for things :-)
>
> 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/1fcdc32c-3943-4c49-b8bc-f5a218f0e2f0n%40googlegroups.com.
