From: kanschat Date: Wed, 8 Aug 2012 14:28:23 +0000 (+0000) Subject: boundary residuals X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=0fc0ef3298d4418e83676c9cef5da53c401d068e;p=dealii-svn.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