convectionmatrix3d.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", "tetrahedron", 1)

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

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