Mercurial > fem-fenics-eugenio

view devel/example/Darcy-Stokes/P2P0/Err_eps.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/Err_eps.ufl@065ae2bff9ef
children
line wrap: on
line source

P2 = VectorElement("CG", triangle, 2)

v = Coefficient(P2)
vh = Coefficient(P2)
ep  = Constant(triangle)

M = (dot(v-vh,v-vh) + div(v-vh)*div(v-vh) + ep*ep*div(v-vh)*div(v-vh) + ep*ep*dot(rot(v-vh), rot(v-vh)))*dx