projection.form

利用C

CG1 = FiniteElement("Lagrange", "tetrahedron", 1)
DG1 = FiniteElement("Discontinuous Lagrange", "tetrahedron", 1)

v = TestFunction(CG1)
u = TrialFunction(CG1)
f = Function(DG1)

a = v*u*dx
L = v*f*dx