Merge branch 'bugfix/stokes_ufl_in_application' into 'master'

Divide through h_e for in Stokes application

See merge request dominic/dune-perftool!209

applications/stokes_dg/stokes_dg.ufl

@@ -16,20 +16,21 @@ ds = ds(subdomain_id=1, subdomain_data=bctype)
 n = FacetNormal(cell)('+')
 eps = -1.0
 sigma = 1.0
+h_e = Min(CellVolume(cell)('+'), CellVolume(cell)('-')) / FacetArea(cell)
 
 r = inner(grad(u), grad(v))*dx \
   - p*div(v)*dx \
   - q*div(u)*dx \
   + inner(avg(grad(u))*n, jump(v))*dS \
-  + sigma * inner(jump(u), jump(v))*dS \
+  + sigma / h_e * inner(jump(u), jump(v))*dS \
   - eps * inner(avg(grad(v))*n, jump(u))*dS \
   - avg(p)*inner(jump(v), n)*dS \
   - avg(q)*inner(jump(u), n)*dS \
   - inner(grad(u)*n, v)*ds \
-  + p*inner(v, n)*ds \
-  + q*inner(u, n)*ds \
+  + sigma / h_e * inner(u-g_v, v)*ds \
   + eps * inner(grad(v)*n, u-g_v)*ds \
-  + sigma * inner(u-g_v, v)*ds \
-  - q*inner(g_v, n)*ds
+  + p*inner(v, n)*ds \
+  + q*inner(u-g_v, n)*ds
 
 forms = [r]
+exact_solution = g_v, 8*(1.-x[0])
\ No newline at end of file