Thanks Ronan, that is very interesting. I checked the source file
getfem_fem.cc,
around lines 1300-1320 where the computation of the Nedelec basis functions for
the reference element is done, but could find no such division by an edge
length.
At that point in the source code I did not find any code that describes how
this
is converted to the actual tetrahedron element (not the reference), but I'm
sure
I'll find it somewhere :-)
I'll also do some tests by multiplying the GetFEM edge basis mass matrix
elements by the corresponding edge lengths to see if I get the same results as
I
get when manually calculating mass matrix elements from the definitions of the
Whitney 1-form basis that I mentioned in the first email.
Thanks again for that very useful info!
Torquil
On 23/09/10 20:23, Ronan Perrussel wrote:
> Hi Torquil,
>
> I have an element which is not an answer to your issue but it can help.
> I never used the RT0 elements but I used the Nedelec elements and the
> shape function associated to edge {ij} is :
> 1/length({ij})*(lambda_i grad lambda_j - lambda_j grad lambda_i).
>
> Maybe something similar is done for RT0 elements (for instance a dof
> associated to a face is not the flux through this face but the mean
> value of the normal component of the vector field on this face).
>
> Best regards,
> Ronan
>
> Le 23/09/2010 20:01, Torquil Macdonald Sørensen a écrit :
>> Hi!
>>
>> I'm having trouble with calculation of mass matrix elements on a 3d
>> tetrahedron
>> mesh. I have two different cases. In the first case I want to use
>> Whitney 1-form
>> edge elements, and in the second case Whitney 2-form face elements.
>> I'm given to
>> understand that these are the same as the NEDELEC and RT0 elements
>> offered by
>> GetFEM.
>>
>> The Whitney elements are defined as follows: For a point i, lambda_i
>> is the
>> usual affine function that is equal to 1 at i, and equal to zero at
>> neighbouring
>> points.
>>
>> For an edge {ij} going from point i to point j, the Whitney 1-form
>> element is
>>
>> lambda_ij := lambda_i grad lambda_j - lambda_j grad lambda_i
>>
>> The "2-form" element for a face {ijk} is
>>
>> lambda_{ijk} := 0.5*(lambda_i (grad lambda_j) x (grad lambda_k) + cyclic
>> permutations of {ijk}).
>>
>> When varying the "fineness"/lattice spacing of the mesh for a given
>> physical
>> domain [0,1]^3, these definitions lead to the following behaviour of
>> the mass
>> matrix components:
>>
>> m_ff' ~ 1/h
>>
>> m_ee' ~ h
>>
>> where m_ff' is a component of the mass matrix for the 2-form elements,
>> m_ee' is
>> the same for the 1-form elements, and h is the lattice spacing.
>>
>> However, in my case, the components m_ff' decreases when h decreases.
>> Does this
>> indicate a difference between the definition of the RT0 element
>> relative to the
>> Whitney 2-form element? If so, is it possible to get one from the
>> other by a
>> simple multiplication by h?
>>
>> Thanks
>> Torquil
>>
>>
>> _______________________________________________
>> Getfem-users mailing list
>> [email protected]
>> https://mail.gna.org/listinfo/getfem-users
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
