From f51a5069507227c9aa8657db0935a53cba294fcb Mon Sep 17 00:00:00 2001
From: Dominic Kempf <dominic.kempf@iwr.uni-heidelberg.de>
Date: Wed, 15 Jun 2016 17:21:45 +0200
Subject: [PATCH] Set up more (almost) full stokes dg example
---
test/stokes/stokes_dg.ufl | 18 ++++++++++++++----
1 file changed, 14 insertions(+), 4 deletions(-)
diff --git a/test/stokes/stokes_dg.ufl b/test/stokes/stokes_dg.ufl
index 5704c65a..71ff6674 100644
--- a/test/stokes/stokes_dg.ufl
+++ b/test/stokes/stokes_dg.ufl
@@ -8,10 +8,20 @@ u, p = TrialFunctions(TH)
n = FacetNormal(cell)('+')
eps = -1.0
+sigma = 1.0
-r = \
- + eps*inner(avg(grad(v))*n, jump(u))*dS
-# + inner(grad(u), grad(v))*dx \
-# - inner(avg(grad(u))*n, jump(v))*dS \
+r = inner(grad(u), grad(v))*dx \
+ - inner(avg(grad(u))*n, jump(v))*dS \
+ + eps * inner(avg(grad(v))*n, jump(u))*dS \
+ - inner(grad(u)*n, v)*ds \
+ + eps * inner(grad(v)*n, u)*ds \
+ + sigma * inner(jump(u), jump(v))*dS \
+ + sigma * inner(u, v)*ds \
+ - p*div(v)*dx \
+ - q*div(u)*dx \
+ + avg(p)*inner(jump(v), n)*dS \
+ + avg(q)*inner(jump(u), n)*dS \
+ + p*inner(v, n)*ds \
+ + q*inner(u, n)*ds
forms = [r]
--
GitLab