Dear Firedrakers, With Toyo Sadare we are looking at implementing a new type of lumping scheme to remove the spectral gap from the RT2-DG1 (Norwegian numbering) spaces for linear shallow water on quadrilaterals. To implement this we need to be able to write some stuff in local coordinates instead of global, including defining coefficients that are functions of local coordinate (stuff like xi_1(1-xi_1) etc.). What is the best way to do this? all the best --cjc
On 28 Feb 2019, at 10:24, Cotter, Colin J <colin.cotter@imperial.ac.uk> wrote:
Dear Firedrakers, With Toyo Sadare we are looking at implementing a new type of lumping scheme to remove the spectral gap from the RT2-DG1 (Norwegian numbering) spaces for linear shallow water on quadrilaterals. To implement this we need to be able to write some stuff in local coordinates instead of global, including defining coefficients that are functions of local coordinate (stuff like xi_1(1-xi_1) etc.).
What is the best way to do this?
from ufl.geometry import CellCoordinate Xhat = CellCoordinate(mesh) Use just like SpatialCoordinate, it's just you get the local rather than global coordinates. Lawrence
On Thu, Feb 28, 2019, 19:25 Cotter, Colin J <colin.cotter@imperial.ac.uk> wrote:
Dear Firedrakers,
With Toyo Sadare we are looking at implementing a new type of lumping scheme to remove the spectral gap from the RT2-DG1 (Norwegian numbering)
Also known as Scandinavian numbering 😎 -- Marie spaces for linear shallow water on quadrilaterals. To implement this we
need to be able to write some stuff in local coordinates instead of global, including defining coefficients that are functions of local coordinate (stuff like xi_1(1-xi_1) etc.).
What is the best way to do this?
all the best
--cjc
_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake
participants (3)
-
Cotter, Colin J
-
Lawrence Mitchell
-
Marie Elisabeth Rognes