Mercurial > fem-fenics-eugenio

view devel/example/Darcy-Stokes/P2P0/DS.ufl @ 161:7e922376b3e8

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Devel folder for examples which need to be improved
author gedeone-octave <marcovass89@hotmail.it>
date Thu, 12 Sep 2013 15:23:48 +0200
parents example/Darcy-Stokes/P2P0/DS.ufl@c15e5434a512
children
line wrap: on
line source

# P2-P0 finite element
P2 = VectorElement("CG", triangle, 2)
P0  = FiniteElement("DG", triangle, 0)

R = FiniteElement ("R", triangle, 0)

W = MixedElement ([P2, P0, R])

(u, p, c) = TrialFunctions(W)
(v, q, d) = TestFunctions(W)

f = Coefficient(P2)
g = Coefficient(P0)
ep  = Constant(triangle)

a = (dot(u, v) + ep*ep*div(u)*div(v) + ep*ep*dot(rot(u), rot(v)) + div(v)*p - div(u)*q + p*d + q*c)*dx
L = (dot (f, v) - g*q)*dx