Hi,
I would like to implement the Raviart-Thomas element in libmesh.
At the moment, I have copied the Nedelec files (needle_one_shape, 2D, 3D),
renamed them as Raviart and change the others ones in order to make it compile.
Generally, to create a new element, I need to define its shape functions, the
transformation between the reference and the actual element and the assembly.
The shapes functions are defined in the raviart_one, raviart_2d, raviart_3D.
Let the bijectivee transformation be x=J X+ b, where x/X is the point on the
actual/reference element. By definition:
actual_phi(x) = 1/det(J) J reference_phi(X)
I should insert this kind of transformation in a file like the
hcurl_fe_transformation (so in hDIV_fe_transformation).
However I do not know where and when the transformation is it called and where
I can access to:
J=[ dx/dxi dx/deta dx/dzeta;
dy/dxi dy/deta dz/dzeta;
dz/dxi dy/deta dz/dzeta]
Furthermore, the new file hDIV_fe_transformation should be called during the
assembly and I would like to know how and where.
Any help would be really appreciated.
Thank you in advance
GR
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
