Ok, thanks Markus :)
Now I have written the following code :
const unsigned int Temperature_dofs_per_face =
fft.get_fe().dofs_per_face;
Temperature_local_dof_indices_face.resize(Temperature_dofs_per_face);
cellT->face(face_no)->get_dof_indices(Temperature_local_dof_indices_face);
but there is a faillure at execution on the last line (get_dof_indices)
Has somebody an idea how to avoid this problem ? (Certainly I got
something wrong!)
Best,
Stéphane
On Fri, 2010-09-10 at 13:52 +0200, Markus Bürg wrote:
> Hello Stéphane,
>
> I think it does not, but you can use cell->get_fe ().dofs_per_face or
> fft.get_fe ().dofs_per_face.
>
> Best Regards,
> Markus
>
>
>
> Am 10.09.10 13:43, schrieb veys:
> > Hi all
> >
> > My question is simple I think :
> > How to evaluate a field (like temperature) at quadrature points of a
> > face "face_no" ?
> >
> > To do this for quadrature points on a cell I use following code :
> >
> > std::vector<unsigned int> Temperature_local_dof_indices;
> > const unsigned int Temperature_dofs_per_cell =
> > cell->get_fe().dofs_per_cell;
> > Temperature_local_dof_indices.resize
> > cell->get_dof_indices(Temperature_local_dof_indices);
> > for(int i=0;i<Temperature_dofs_per_cell;i++){
> > int dof = Temperature_local_dof_indices[i];
> > T += Temperature_old(dof) * Temperature_fe_values.shape_value(i,q);
> > }
> >
> >
> > I saw that it's possible to write :
> > const FEFaceValues<dim> &fft= hp_fft.get_present_fe_values();
> > const unsigned int Temperature_dofs_per_cell = fft.dofs_per_cell;
> > but does it exist "fft.dofs_per_face" ? or something equivalent for
> > face ?
> >
> > Best,
> > Stéphane
> >
> > _______________________________________________
> > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
> >
> _______________________________________________
> dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
