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
