1
# Copyright (C) 2010 Anders Logg (logg@simula.no).
2
# Licensed under the GNU GPL version 3 or any later version.
4
scalar = FiniteElement("Lagrange", triangle, 3)
5
vector = VectorElement("Lagrange", triangle, 3)
7
v = TestFunction(vector)
8
u = TrialFunction(vector)
9
w = Coefficient(vector)
10
rho = Coefficient(scalar)
12
a = rho*inner(v, grad(w)*u)*dx