Hi,
The attachment is my small program that I created.
My fe system have 24 dofs for both real and imaginary of 12 edges of a cell.
No matter how I change the parameter "order" the number of dofs is
unchanged.
I think I have problem when constructing the object.
I would like to thank you very much.
On 2/20/19 7:01 PM, Phạm Ngọc Kiên wrote:
> I tested in my codes using
>
> template
> CSEM::CSEM()
> ://mapping(1),
> dof_handler(triangulation),
> fe(FE_NedelecSZ(order),1,//(order), multiplicity
> FE_NedelecSZ(order),1) {}
>
> When I changed the parameter order, I saw only the
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
CSEM::CSEM()
://mapping(1),
dof_handler(triangulation),
fe(FE_NedelecSZ(order), 1, //(order), multiplic
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
Yes.
Let me go from my basic idea when setting up my source term.
First, I create a computational domain such that a source can be placed
exactly on an edge. It is always possible to do with unstructured meshes.
Then I go through each element and check if its edge overlap the source. If
it does the
On 2/19/19 12:03 AM, Phạm Ngọc Kiên wrote:
> I currently try to solve the electromagnetics problem, followed by the
> research of Ross Kynch.
> My goal is to compute the right hand side for the problem:
> curl curl E + i*omega*mu*sigma*E = - Js
> I only want to set the Js term here to be none zero
Dear Prof. Wolfgang Bangerth,
I would like to thank you very much.
I am reading through the FE_NedelecSZ codes in Deal.II library.
I currently try to solve the electromagnetics problem, followed by the
research of Ross Kynch.
My goal is to compute the right hand side for the problem:
curl curl E +
On 2/18/19 10:47 PM, Phạm Ngọc Kiên wrote:
>
> 2. I have tried some thing with the third method. And here-below is my code:
>
> Point p{0.5,0.5,0.5};//position in reference cell
> Quadrature q(p);
> FEValues fe_values_q(fe, q,update_quadrature_points);
> fe_values_q.reinit(cell);
>
> //position
I am sorry as my question is not clear enough.
What I mean is when we loop over all degree of freedom in a cell, for
example,
for (unsigned int i = 0; i < dofs_per_cell; ++i) { .. }
Do we have any method to get the position of the i_th degree of
freedom inside this loop?
I would like to
Hi,
I have tried with the methods described in the FAQ and figured out that:
1. As I used edge-based element (FE_NedelecSZ), I queried if my element had
support points by *fe.has_support_points()*. As a consequence, it returned
zero.
The *FiniteElement::get_unit_support_points() *also returned an
Yes, thank you very much.
I will dig deeper into this.
Best,
Pham Ngoc Kien
Vào Th 3, 19 thg 2, 2019 vào lúc 11:23 Bruno Turcksin <
bruno.turck...@gmail.com> đã viết:
> Pham Ngoc Kien,
>
> Le lun. 18 févr. 2019 à 19:33, Phạm Ngọc Kiên
> a écrit :
>
>> With the DoFTools::map_dofs_to_support_poi
Pham Ngoc Kien,
Le lun. 18 févr. 2019 à 19:33, Phạm Ngọc Kiên a
écrit :
> With the DoFTools::map_dofs_to_support_points() , I saw in the library
> manual that it can be used with no edge elements or the like.
> As my finite element is edge-based one, I wonder if this function work?
>
No the func
Hi Bruno,
With the DoFTools::map_dofs_to_support_points() , I saw in the library
manual that it can be used with no edge elements or the like.
As my finite element is edge-based one, I wonder if this function work?
I have check the function fe.has_support_points() in my code and it
returned 0.
That
13 matches
Mail list logo