Use tensor in poisson_dg_tensor test and disable test due to #122

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
 

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