scalar boundary residuals

git-svn-id: https://svn.dealii.org/trunk@25773 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/tests/integrators/laplacian_01.cc b/tests/integrators/laplacian_01.cc
index 6ce7df48931ebd248672751e8b91ecbe5d1f2837..80b13e99b300aedf246d42b91daaf40bfa03cfd3 100644 (file)
--- a/tests/integrators/laplacian_01.cc
+++ b/tests/integrators/laplacian_01.cc
@@ -63,12 +63,20 @@ void test_boundary(const FEValuesBase<dim>& fev)
     {
       u = 0.;
       u(i) = 1.;
+      w = 0.;
       fev.get_function_values(u, indices, uval, true);
       fev.get_function_gradients(u, indices, ugrad, true);
       nitsche_residual(w, fev, make_slice(uval), make_slice(ugrad), make_slice(null_val), 17);
       M.vmult(v,u);
       w.add(-1., v);
       deallog << ' ' << w.l2_norm();
+      if (d==1)
+       {
+         nitsche_residual(w, fev, uval[0], ugrad[0], null_val[0], 17);
+         M.vmult(v,u);
+         w.add(-1., v);
+         deallog << " e" << w.l2_norm();         
+       }
     }
   deallog << std::endl;
 }