Dear Bastian and Daniel,
> That's what I was assuming. So yesterday I did it with this code :
> // This block add to the periodicity constraint the little
> compression we want
> {
> IndexSet selected_dofs_left;
> std::set< types::boundary_id > boundary_ids_left=
>
Dear Daniel,
When I am at the face I need the solution at the interior quadrature
points as well. However from dinfo1.indices we could get the the dof
indices for the cell1. But is it possible to get the shape value at the
interior quadrature points ?
Thanks a lot
On Friday, Nove
Sudarshan,
template
> void Step12::integrate_face_term (DoFInfo& dinfo1, DoFInfo& dinfo2,
>CellInfo& info1, CellInfo&
> info2)
> {
>const FEValuesBase& fe_v = info1.fe_values();
>const FEValuesBase& fe_v_neighbor = info2.fe_values();
>
