From 2cdd5e702da032b1fa8e24077bddc280b4e8fc81 Mon Sep 17 00:00:00 2001 From: Guido Kanschat Date: Wed, 8 Aug 2012 14:28:23 +0000 Subject: [PATCH] boundary residuals git-svn-id: https://svn.dealii.org/trunk@25772 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/integrators/laplacian_01.cc | 38 +++++++++++++++------- tests/integrators/laplacian_01/cmp/generic | 4 +++ 2 files changed, 30 insertions(+), 12 deletions(-) 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