From: kanschat Date: Wed, 8 Aug 2012 16:33:45 +0000 (+0000) Subject: add ip residuals X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=5463147462cdb93cecd63a738c86c6e1c458ab7b;p=dealii-svn.git add ip residuals git-svn-id: https://svn.dealii.org/trunk@25781 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/integrators/laplacian_01.cc b/tests/integrators/laplacian_01.cc index e5261e300d..10bbb9bb00 100644 --- a/tests/integrators/laplacian_01.cc +++ b/tests/integrators/laplacian_01.cc @@ -114,11 +114,11 @@ void test_face(const FEValuesBase& fev1, { const unsigned int n1 = fev1.dofs_per_cell; const unsigned int n2 = fev1.dofs_per_cell; + unsigned int d=fev1.get_fe().n_components(); FullMatrix M11(n1,n1); FullMatrix M12(n1,n2); FullMatrix M21(n2,n1); FullMatrix M22(n2,n2); -// Vector u(n), v(n); ip_matrix(M11, M12, M21, M22, fev1, fev2, 17); deallog << "M11" << std::endl; @@ -129,6 +129,55 @@ void test_face(const FEValuesBase& fev1, M12.print(deallog,8); deallog << "M21" << std::endl; M21.print(deallog,8); + + Vector u1(n1), v1(n1), w1(n1); + Vector u2(n2), v2(n2), w2(n2); + std::vector > + u1val (d,std::vector(fev1.n_quadrature_points)), + u2val (d,std::vector(fev2.n_quadrature_points)); + std::vector > > + u1grad(d,std::vector >(fev1.n_quadrature_points)), + u2grad(d,std::vector >(fev2.n_quadrature_points)); + + std::vector indices1(n1), indices2(n2); + for (unsigned int i=0;i& tr, FiniteElement& fe) cell1 = tr.begin(1); fef1.reinit(cell1, 1); fef2.reinit(cell2, 0); - test_face(fef1, fef2); - - + test_face(fef1, fef2); } diff --git a/tests/integrators/laplacian_01/cmp/generic b/tests/integrators/laplacian_01/cmp/generic index 6e0e9734f5..059b0eb97e 100644 --- a/tests/integrators/laplacian_01/cmp/generic +++ b/tests/integrators/laplacian_01/cmp/generic @@ -52,6 +52,7 @@ DEAL::-0.17 -5.3 -0.083 -2.7 DEAL::0 -0.17 0 -0.083 DEAL::-0.083 -2.7 -0.17 -5.3 DEAL::0 -0.083 0 -0.17 +DEAL::Residuals 0 e0 0 e0 0 e0 0 e0 0 e0 0 e0 0 e0 0 e0 0 e0 0 e0 0 e0 0 e0 0 e0 0 e0 0 e0 0 e0 DEAL::FE_Nedelec<2>(1) DEAL::cell matrix DEAL::1.0 0 -1.0 0 0 0 0 0 0 0 0 0 @@ -176,3 +177,4 @@ DEAL::0 0 0 0 2.5 4.0 2.5 4.0 -1.7 -2 DEAL::0 0 0 0 -4.0 -6.5 -4.0 -6.5 2.8 4.5 0 0 DEAL::0 0 0.87 0 0 0 0 0 0 0 0 0 DEAL::0 0 0 0.87 0 0 0 0 0 0 0 0 +DEAL::Residuals 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0