From: Guido Kanschat Date: Wed, 8 Aug 2012 14:28:23 +0000 (+0000) Subject: boundary residuals X-Git-Tag: v8.0.0~2376 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=2cdd5e702da032b1fa8e24077bddc280b4e8fc81;p=dealii.git boundary residuals git-svn-id: https://svn.dealii.org/trunk@25772 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/integrators/laplacian_01.cc b/tests/integrators/laplacian_01.cc index 4836a589d4..6ce7df4893 100644 --- a/tests/integrators/laplacian_01.cc +++ b/tests/integrators/laplacian_01.cc @@ -40,23 +40,37 @@ template void test_boundary(const FEValuesBase& fev) { const unsigned int n = fev.dofs_per_cell; + unsigned int d=fev.get_fe().n_components(); + FullMatrix M(n,n); + Vector u(n), v(n), w(n); + std::vector > + uval (d,std::vector(fev.n_quadrature_points)), + null_val(d,std::vector(fev.n_quadrature_points, 0.)); + std::vector > > + ugrad (d,std::vector >(fev.n_quadrature_points)); + + std::vector indices(n); + for (unsigned int i=0;i