Dear Firedrakers, on an extruded mesh: meshbase = RectangleMesh(Nx, Ny, Lx, Ly, quadrilateral=True) mesh = ExtrudedMesh(meshbase, Nz, Delta_z) horiz_elt = FiniteElement("DG", quadrilateral, 0) vert_elt = FiniteElement("DG", interval, 0) elt = OuterProductElement(horiz_elt, vert_elt) How can I access the element center and the center of the element faces? Thank you! Henrik -- Dipl.-Math. Henrik Büsing Applied Geophysics and Geothermal Energy E.ON Energy Research Center RWTH Aachen University ----------------------------------------------- Mathieustr. 10 | Tel +49 (0)241 80 49907 52074 Aachen, Germany | Fax +49 (0)241 80 49889 ----------------------------------------------- http://www.eonerc.rwth-aachen.de/GGE hbuesing@eonerc.rwth-aachen.de -----------------------------------------------
-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Dear Henrik, On 22/09/15 17:52, Buesing, Henrik wrote:
Dear Firedrakers,
on an extruded mesh:
meshbase = RectangleMesh(Nx, Ny, Lx, Ly, quadrilateral=True) mesh = ExtrudedMesh(meshbase, Nz, Delta_z) horiz_elt = FiniteElement("DG", quadrilateral, 0) vert_elt = FiniteElement("DG", interval, 0) elt = OuterProductElement(horiz_elt, vert_elt)
How can I access the element center and the center of the element faces? Thank you!
I think we have no symbolic way of doing this. It is possible to do so by writing some C code, but before we try and do so, I wonder to what end you would like these values? Cheers, Lawrence -----BEGIN PGP SIGNATURE----- Version: GnuPG v1 iQEcBAEBAgAGBQJWAYjJAAoJECOc1kQ8PEYv5mAIAIb06bApIlIAFVfJ+ECLEtUn BjMOAk2puGq7tKHfN3eWD5l9Z5dWoPu7Sowbx07CXjpm6rW+XxFnNLk1QLIYr9n2 TBAq7BoKrJ8wLBTfSfoYzb7kxp5HOLbUMjz2Qht+mePaKcU/RLSr35WUinvtlUfP Pk+VugohUBIB6OxO1uK+qFCtMdZxrQq6dvWR7gZMfUVuiv9YNm6emdSGuK9MrUcK ifPsz/9D4/hRBkhCqhclihvbRVF7YXGnr4zZfHadh2ZRTm/RyCy5TLtqGSjY4rUv 79I/SMKbletXUDc2a8wlDyaCXBvEZvt5K39GUkLGs9f8MHBAxTAyhZgqUz477mI= =Zdn9 -----END PGP SIGNATURE-----
participants (2)
-
Buesing, Henrik
-
Lawrence Mitchell