From: kanschat Date: Wed, 8 Aug 2012 14:36:32 +0000 (+0000) Subject: scalar boundary residuals X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=a0bea89a827f8166ae89cb49d035f889ed357205;p=dealii-svn.git 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 6ce7df4893..80b13e99b3 100644 --- a/tests/integrators/laplacian_01.cc +++ b/tests/integrators/laplacian_01.cc @@ -63,12 +63,20 @@ void test_boundary(const FEValuesBase& 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; }