Dear Daniel,
I think that I can address my problem with the source now after looking
through the way we create dofs in FE_NedelecSZ.
I tested in my codes using
template<int dim>
CSEM<dim>::CSEM()
://mapping(1),
dof_handler(triangulation),
fe(FE_NedelecSZ<dim>(order), 1, //(order), multiplicity
FE_NedelecSZ<dim>(order), 1) {}
When I changed the parameter order, I saw only the number of dofs was
only of lowest order dofs no matter how big the order was
I think that there should be a way to get higher order dofs, but I
cannot find it now.
Thank you very much.
Vào Th 4, 20 thg 2, 2019 vào lúc 18:54 Daniel Arndt <
daniel.ar...@iwr.uni-heidelberg.de> đã viết:
> 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.
>
--
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.
