28 Feb
2019
28 Feb
'19
6:31 p.m.
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