1
# Copyright (C) 2006 Garth N. Wells.
2
# Licensed under the GNU GPL Version 2.
4
# First added: 2006-08-18
7
# Bilinear form a(v, U) for vector-valued Poisson equation
10
element = FiniteElement("Vector Lagrange", "triangle", 1)
12
v = TestFunction(element)
13
U = TrialFunction(element)
16
# Use symmetric gradient to increase coupling which is more demanding
20
return 0.5 * (grad(v) + transp(grad(v)))
22
a = dot(symgrad(v), symgrad(U))*dx