convectionmatrix2d.form

利用C

# Copyright (C) 2006 Anders Logg ([e-mail])
# Licensed under the GNU LGPL Version 2.1
#
# First added:  2006-08-21
# Last changed: 2006-08-21
#
# The bilinear form for a convection matrix

element = FiniteElement("Lagrange", "triangle", 1)

v  = TestFunction(element)
U  = TrialFunction(element)
cx = Constant("triangle")
cy = Constant("triangle")

a = v*(cx*U.dx(0) + cy*U.dx(1))*dx