Hello Stéphane,
you are using the hp-namespace. There you can have two different finite
elements associated to a face, one from each cell sharing it. Therefore
you have to tell the compiler, which one you want to use:
cellT->face(face_no)->get_dof_indices(Temperature_local_dof_indices_face,
cellT->active_fe_index ());
Best Regards,
Markus
Am 10.09.10 14:04, schrieb veys:
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
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
