Kien, I am not quite sure I understand how the source term you want to use looks like.
If you want to describe a point source, you should gave a look at https://www.dealii.org/developer/doxygen/deal.II/namespaceVectorTools.html#ac4ffc41e72e846029d04e7034de60d09. Be aware of the last paragraphs in its documentation, though. On the other hand, if you have a usual right hand side Js(x), you need to follow the approach Wolfgang described and that is used in almost all tutorial progams, e.g. https://www.dealii.org/developer/doxygen/deal.II/step_4.html#Step4assemble_system. How does Js(x) look like? I also tried to use FE_NedelecSZ<dim>(order) with different order. > And I saw that no matter how I change the order value, the functions > fe.n_dofs_per_quad() and fe.n_dofs_per_hex() returned 0. The number of dofs > was unchanged and was equal to fe.n_dofs_per_line(). > Does this mean that there is no face and cell dof? > According to the implementation there should be - 0 dofs per vertex - degree+1 dofs per line - 2*degree*(degree+1) dofs per quad - 3*degree*degree*(degree+1) dofs per hex > Do I need to use hp::DoFHandler to get the higher order dofs for > calculation? > No, you would use this if you want to use a different ansatz space, e.g. the same type of element with a different polynomial degree, for each cell. Best, Daniel -- 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
