From 0fc0ef3298d4418e83676c9cef5da53c401d068e Mon Sep 17 00:00:00 2001
From: kanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
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 <int dim>
 void test_boundary(const FEValuesBase<dim>& fev)
 {
   const unsigned int n = fev.dofs_per_cell;
+  unsigned int d=fev.get_fe().n_components();
+  
   FullMatrix<double> M(n,n);
+  
   Vector<double> u(n), v(n), w(n);
+  std::vector<std::vector<double> >
+    uval    (d,std::vector<double>(fev.n_quadrature_points)),
+    null_val(d,std::vector<double>(fev.n_quadrature_points, 0.));
+  std::vector<std::vector<Tensor<1,dim> > >
+    ugrad   (d,std::vector<Tensor<1,dim> >(fev.n_quadrature_points));
+  
+  std::vector<unsigned int> indices(n);
+  for (unsigned int i=0;i<n;++i)
+    indices[i] = i;
   
   nitsche_matrix(M, fev, 17);
   M.print(deallog,8);
-
-  // deallog << "Residuals";
-  // for (unsigned int i=0;i<n;++i)
-  //   {
-  //     u = 0.;
-  //     u(i) = 1.;
-  //     M.vmult(v,u);
-  //     nitsche_residual(w, fev, 17);
-  //     w.add(-1., v);
-  //     deallog << ' ' << w.l2_norm();
-  //   }
-  // deallog << std::endl;
+  
+  deallog << "Residuals";
+  for (unsigned int i=0;i<n;++i)
+    {
+      u = 0.;
+      u(i) = 1.;
+      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();
+    }
+  deallog << std::endl;
 }
 
 
diff --git a/tests/integrators/laplacian_01/cmp/generic b/tests/integrators/laplacian_01/cmp/generic
index 476f291320..0208202d15 100644
--- a/tests/integrators/laplacian_01/cmp/generic
+++ b/tests/integrators/laplacian_01/cmp/generic
@@ -11,21 +11,25 @@ DEAL::11.     0.33    5.3     0.17
 DEAL::0.33    0       0.17    0       
 DEAL::5.3     0.17    11.     0.33    
 DEAL::0.17    0       0.33    0       
+DEAL::Residuals 0 0 0 0
 DEAL::boundary_matrix 1
 DEAL::0       0.33    0       0.17    
 DEAL::0.33    11.     0.17    5.3     
 DEAL::0       0.17    0       0.33    
 DEAL::0.17    5.3     0.33    11.     
+DEAL::Residuals 0 0 0 0
 DEAL::boundary_matrix 2
 DEAL::11.     5.3     0.33    0.17    
 DEAL::5.3     11.     0.17    0.33    
 DEAL::0.33    0.17    0       0       
 DEAL::0.17    0.33    0       0       
+DEAL::Residuals 0 0 0 0
 DEAL::boundary_matrix 3
 DEAL::0       0       0.33    0.17    
 DEAL::0       0       0.17    0.33    
 DEAL::0.33    0.17    11.     5.3     
 DEAL::0.17    0.33    5.3     11.     
+DEAL::Residuals 0 0 0 0
 DEAL::face_matrix 0
 DEAL::M11
 DEAL::0       0.17    0       0.083   
-- 
2.39.5