From d4611ab6295d248791aa79ce34e157ac583304fe Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Ren=C3=A9=20He=C3=9F?= <rene.hess@iwr.uni-heidelberg.de>
Date: Thu, 2 Aug 2018 13:32:27 +0200
Subject: [PATCH] Use tensor in poisson_dg_tensor test and disable test due to
 #122

---
 test/sumfact/poisson/CMakeLists.txt        |  8 ++++----
 test/sumfact/poisson/poisson_dg_tensor.ufl | 14 +++++++-------
 2 files changed, 11 insertions(+), 11 deletions(-)

diff --git a/test/sumfact/poisson/CMakeLists.txt b/test/sumfact/poisson/CMakeLists.txt
index 824f2721..8a06e3f0 100644
--- a/test/sumfact/poisson/CMakeLists.txt
+++ b/test/sumfact/poisson/CMakeLists.txt
@@ -49,10 +49,10 @@ dune_add_formcompiler_system_test(UFLFILE poisson_dg_3d.ufl
 #================================
 # 'Poisson' DG with a full tensor
 #================================
-dune_add_formcompiler_system_test(UFLFILE poisson_dg_tensor.ufl
-                                  BASENAME sumfact_poisson_dg_tensor
-                                  INIFILE poisson_dg_tensor.mini
-                                  )
+# dune_add_formcompiler_system_test(UFLFILE poisson_dg_tensor.ufl
+#                                   BASENAME sumfact_poisson_dg_tensor
+#                                   INIFILE poisson_dg_tensor.mini
+#                                   )
 
 # Slicing vectorization stratgies
 
diff --git a/test/sumfact/poisson/poisson_dg_tensor.ufl b/test/sumfact/poisson/poisson_dg_tensor.ufl
index 2734383f..50af4173 100644
--- a/test/sumfact/poisson/poisson_dg_tensor.ufl
+++ b/test/sumfact/poisson/poisson_dg_tensor.ufl
@@ -25,15 +25,15 @@ gamma_int = (alpha * degree * (degree + dim - 1)) / h_int
 # SIPG: -1.0, IIPG: 0.0, NIPG: 1.0
 theta = 1.0
 
-r = inner(grad(u), grad(v))*dx \
-  - f*v*dx \
-  - inner(n, avg(grad(u)))*jump(v)*dS \
+r = inner(A*grad(u), grad(v))*dx \
+  - (c*u-f)*v*dx \
+  - inner(n, A*avg(grad(u)))*jump(v)*dS \
   + gamma_int*jump(u)*jump(v)*dS \
-  + theta*jump(u)*inner(avg(grad(v)), n)*dS \
-  - inner(n, grad(u))*v*ds \
+  + theta*jump(u)*inner(A*avg(grad(v)), n)*dS \
+  - inner(n, A*grad(u))*v*ds \
   + gamma_ext*u*v*ds \
-  + theta*u*inner(grad(v), n)*ds \
+  + theta*u*inner(A*grad(v), n)*ds \
   - gamma_ext*g*v*ds \
-  - theta*g*inner(grad(v), n)*ds
+  - theta*g*inner(A*grad(v), n)*ds
 
 exact_solution = g
-- 
GitLab