From: kayser-herold Date: Sat, 29 Jul 2006 21:59:02 +0000 (+0000) Subject: Check same with FEFaceValues. Results are presently wrong (no scaling, when it should... X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=2c7cf5b0b130871192da2bc1eb4e68619acc6d8e;p=dealii-svn.git Check same with FEFaceValues. Results are presently wrong (no scaling, when it should be 1/h git-svn-id: https://svn.dealii.org/trunk@13507 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/fe/rt_11.cc b/tests/fe/rt_11.cc new file mode 100644 index 0000000000..69a5a2d0c0 --- /dev/null +++ b/tests/fe/rt_11.cc @@ -0,0 +1,100 @@ +//---------------------------- rt_11.cc --------------------------- +// rt_11.cc,v 1.3 2003/06/09 16:00:38 wolf Exp +// Version: +// +// Copyright (C) 2003, 2005, 2006 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------- rt_11.cc --------------------------- + +// Observe how the values of the shape functions change as we make a +// cell smaller and smaller. Evaluate the values with FEFaceValues, to +// make sure the values scale as in rt_10 where we used FEValues + +#include "../tests.h" +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#include +#include +#include + +#define PRECISION 2 + + +std::ofstream logfile ("rt_11/output"); + +template +void +test (const unsigned int degree) +{ + FE_RaviartThomas fe_rt(degree); + + deallog << "Degree=" << degree + << std::endl; + + for (double h=1; h>1./128; h/=2) + { + deallog << " h=" << h + << std::endl; + + Triangulation tr; + GridGenerator::hyper_cube(tr, 0., h); + + DoFHandler dof(tr); + dof.distribute_dofs(fe_rt); + + QTrapez quadrature; + + FEFaceValues fe_values (fe_rt, quadrature, update_values); + fe_values.reinit (dof.begin_active(), 0); + for (unsigned int q=0; q(i); + + return 0; +} + + + diff --git a/tests/fe/rt_11/cmp/generic b/tests/fe/rt_11/cmp/generic new file mode 100644 index 0000000000..f971f84f17 --- /dev/null +++ b/tests/fe/rt_11/cmp/generic @@ -0,0 +1,145 @@ + +DEAL::Degree=0 +DEAL:: h=1.00 +DEAL:: Quadrature point 0: [1.00 0 ][0 0 ][0 1.00 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][1.00 0 ][0 1.00 ][0 0 ] +DEAL:: Quadrature point 2: [1.00 0 ][0 0 ][0 0 ][0 1.00 ] +DEAL:: Quadrature point 3: [0 0 ][1.00 0 ][0 0 ][0 1.00 ] +DEAL:: h=0.50 +DEAL:: Quadrature point 0: [2.00 0 ][0 0 ][0 2.00 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][2.00 0 ][0 2.00 ][0 0 ] +DEAL:: Quadrature point 2: [2.00 0 ][0 0 ][0 0 ][0 2.00 ] +DEAL:: Quadrature point 3: [0 0 ][2.00 0 ][0 0 ][0 2.00 ] +DEAL:: h=0.25 +DEAL:: Quadrature point 0: [4.00 0 ][0 0 ][0 4.00 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][4.00 0 ][0 4.00 ][0 0 ] +DEAL:: Quadrature point 2: [4.00 0 ][0 0 ][0 0 ][0 4.00 ] +DEAL:: Quadrature point 3: [0 0 ][4.00 0 ][0 0 ][0 4.00 ] +DEAL:: h=0.12 +DEAL:: Quadrature point 0: [8.00 0 ][0 0 ][0 8.00 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][8.00 0 ][0 8.00 ][0 0 ] +DEAL:: Quadrature point 2: [8.00 0 ][0 0 ][0 0 ][0 8.00 ] +DEAL:: Quadrature point 3: [0 0 ][8.00 0 ][0 0 ][0 8.00 ] +DEAL:: h=0.06 +DEAL:: Quadrature point 0: [16.00 0 ][0 0 ][0 16.00 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][16.00 0 ][0 16.00 ][0 0 ] +DEAL:: Quadrature point 2: [16.00 0 ][0 0 ][0 0 ][0 16.00 ] +DEAL:: Quadrature point 3: [0 0 ][16.00 0 ][0 0 ][0 16.00 ] +DEAL:: h=0.03 +DEAL:: Quadrature point 0: [32.00 0 ][0 0 ][0 32.00 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][32.00 0 ][0 32.00 ][0 0 ] +DEAL:: Quadrature point 2: [32.00 0 ][0 0 ][0 0 ][0 32.00 ] +DEAL:: Quadrature point 3: [0 0 ][32.00 0 ][0 0 ][0 32.00 ] +DEAL:: h=0.02 +DEAL:: Quadrature point 0: [64.00 0 ][0 0 ][0 64.00 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][64.00 0 ][0 64.00 ][0 0 ] +DEAL:: Quadrature point 2: [64.00 0 ][0 0 ][0 0 ][0 64.00 ] +DEAL:: Quadrature point 3: [0 0 ][64.00 0 ][0 0 ][0 64.00 ] +DEAL::Degree=1 +DEAL:: h=1.00 +DEAL:: Quadrature point 0: [1.00 0 ][-1.73 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][1.00 0 ][-1.73 0 ][0 1.00 ][0 1.73 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 2: [1.00 0 ][1.73 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][1.00 0 ][1.73 0 ][0 0 ][0 0 ][0 1.00 ][0 1.73 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: h=0.50 +DEAL:: Quadrature point 0: [2.00 0 ][-3.46 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][2.00 0 ][-3.46 0 ][0 2.00 ][0 3.46 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 2: [2.00 0 ][3.46 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][2.00 0 ][3.46 0 ][0 0 ][0 0 ][0 2.00 ][0 3.46 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: h=0.25 +DEAL:: Quadrature point 0: [4.00 0 ][-6.93 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][4.00 0 ][-6.93 0 ][0 4.00 ][0 6.93 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 2: [4.00 0 ][6.93 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][4.00 0 ][6.93 0 ][0 0 ][0 0 ][0 4.00 ][0 6.93 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: h=0.12 +DEAL:: Quadrature point 0: [8.00 0 ][-13.86 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][8.00 0 ][-13.86 0 ][0 8.00 ][0 13.86 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 2: [8.00 0 ][13.86 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][8.00 0 ][13.86 0 ][0 0 ][0 0 ][0 8.00 ][0 13.86 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: h=0.06 +DEAL:: Quadrature point 0: [16.00 0 ][-27.71 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][16.00 0 ][-27.71 0 ][0 16.00 ][0 27.71 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 2: [16.00 0 ][27.71 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][16.00 0 ][27.71 0 ][0 0 ][0 0 ][0 16.00 ][0 27.71 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: h=0.03 +DEAL:: Quadrature point 0: [32.00 0 ][-55.43 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][32.00 0 ][-55.43 0 ][0 32.00 ][0 55.43 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 2: [32.00 0 ][55.43 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][32.00 0 ][55.43 0 ][0 0 ][0 0 ][0 32.00 ][0 55.43 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: h=0.02 +DEAL:: Quadrature point 0: [64.00 0 ][-110.85 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][64.00 0 ][-110.85 0 ][0 64.00 ][0 110.85 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 2: [64.00 0 ][110.85 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][64.00 0 ][110.85 0 ][0 0 ][0 0 ][0 64.00 ][0 110.85 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL::Degree=2 +DEAL:: h=1.00 +DEAL:: Quadrature point 0: [1.00 0 ][-1.73 0 ][2.24 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 2.24 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][1.00 0 ][-1.73 0 ][2.24 0 ][0 1.00 ][0 1.73 ][0 2.24 ][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 ] +DEAL:: Quadrature point 2: [1.00 0 ][1.73 0 ][2.24 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 2.24 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][1.00 0 ][1.73 0 ][2.24 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 1.73 ][0 2.24 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: h=0.50 +DEAL:: Quadrature point 0: [2.00 0 ][-3.46 0 ][4.47 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 4.47 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][2.00 0 ][-3.46 0 ][4.47 0 ][0 2.00 ][0 3.46 ][0 4.47 ][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 ] +DEAL:: Quadrature point 2: [2.00 0 ][3.46 0 ][4.47 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 4.47 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][2.00 0 ][3.46 0 ][4.47 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 3.46 ][0 4.47 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: h=0.25 +DEAL:: Quadrature point 0: [4.00 0 ][-6.93 0 ][8.94 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 8.94 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][4.00 0 ][-6.93 0 ][8.94 0 ][0 4.00 ][0 6.93 ][0 8.94 ][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 ] +DEAL:: Quadrature point 2: [4.00 0 ][6.93 0 ][8.94 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 8.94 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][4.00 0 ][6.93 0 ][8.94 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 6.93 ][0 8.94 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: h=0.12 +DEAL:: Quadrature point 0: [8.00 0 ][-13.86 0 ][17.89 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 17.89 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][8.00 0 ][-13.86 0 ][17.89 0 ][0 8.00 ][0 13.86 ][0 17.89 ][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 ] +DEAL:: Quadrature point 2: [8.00 0 ][13.86 0 ][17.89 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 17.89 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][8.00 0 ][13.86 0 ][17.89 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 13.86 ][0 17.89 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: h=0.06 +DEAL:: Quadrature point 0: [16.00 0 ][-27.71 0 ][35.78 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 35.78 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][16.00 0 ][-27.71 0 ][35.78 0 ][0 16.00 ][0 27.71 ][0 35.78 ][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 ] +DEAL:: Quadrature point 2: [16.00 0 ][27.71 0 ][35.78 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 35.78 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][16.00 0 ][27.71 0 ][35.78 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 27.71 ][0 35.78 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: h=0.03 +DEAL:: Quadrature point 0: [32.00 0 ][-55.43 0 ][71.55 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 71.55 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][32.00 0 ][-55.43 0 ][71.55 0 ][0 32.00 ][0 55.43 ][0 71.55 ][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 ] +DEAL:: Quadrature point 2: [32.00 0 ][55.43 0 ][71.55 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 71.55 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][32.00 0 ][55.43 0 ][71.55 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 55.43 ][0 71.55 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: h=0.02 +DEAL:: Quadrature point 0: [64.00 0 ][-110.85 0 ][143.11 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 143.11 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][64.00 0 ][-110.85 0 ][143.11 0 ][0 64.00 ][0 110.85 ][0 143.11 ][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 ] +DEAL:: Quadrature point 2: [64.00 0 ][110.85 0 ][143.11 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 143.11 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][64.00 0 ][110.85 0 ][143.11 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 110.85 ][0 143.11 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ] +DEAL::Degree=3 +DEAL:: h=1.00 +DEAL:: Quadrature point 0: [1.00 0 ][-1.73 0 ][2.24 0 ][-2.65 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 2.24 ][0 -2.65 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][1.00 0 ][-1.73 0 ][2.24 0 ][-2.65 0 ][0 1.00 ][0 1.73 ][0 2.24 ][0 2.65 ][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 ] +DEAL:: Quadrature point 2: [1.00 0 ][1.73 0 ][2.24 0 ][2.65 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 -1.73 ][0 2.24 ][0 -2.65 ][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 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][1.00 0 ][1.73 0 ][2.24 0 ][2.65 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 1.00 ][0 1.73 ][0 2.24 ][0 2.65 ][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 ] +DEAL:: h=0.50 +DEAL:: Quadrature point 0: [2.00 0 ][-3.46 0 ][4.47 0 ][-5.29 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 4.47 ][0 -5.29 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][2.00 0 ][-3.46 0 ][4.47 0 ][-5.29 0 ][0 2.00 ][0 3.46 ][0 4.47 ][0 5.29 ][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 ] +DEAL:: Quadrature point 2: [2.00 0 ][3.46 0 ][4.47 0 ][5.29 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 -3.46 ][0 4.47 ][0 -5.29 ][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 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][2.00 0 ][3.46 0 ][4.47 0 ][5.29 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 2.00 ][0 3.46 ][0 4.47 ][0 5.29 ][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 ] +DEAL:: h=0.25 +DEAL:: Quadrature point 0: [4.00 0 ][-6.93 0 ][8.94 0 ][-10.58 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 8.94 ][0 -10.58 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][4.00 0 ][-6.93 0 ][8.94 0 ][-10.58 0 ][0 4.00 ][0 6.93 ][0 8.94 ][0 10.58 ][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 ] +DEAL:: Quadrature point 2: [4.00 0 ][6.93 0 ][8.94 0 ][10.58 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 -6.93 ][0 8.94 ][0 -10.58 ][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 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][4.00 0 ][6.93 0 ][8.94 0 ][10.58 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 4.00 ][0 6.93 ][0 8.94 ][0 10.58 ][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 ] +DEAL:: h=0.12 +DEAL:: Quadrature point 0: [8.00 0 ][-13.86 0 ][17.89 0 ][-21.17 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 17.89 ][0 -21.17 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][8.00 0 ][-13.86 0 ][17.89 0 ][-21.17 0 ][0 8.00 ][0 13.86 ][0 17.89 ][0 21.17 ][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 ] +DEAL:: Quadrature point 2: [8.00 0 ][13.86 0 ][17.89 0 ][21.17 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 -13.86 ][0 17.89 ][0 -21.17 ][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 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][8.00 0 ][13.86 0 ][17.89 0 ][21.17 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 8.00 ][0 13.86 ][0 17.89 ][0 21.17 ][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 ] +DEAL:: h=0.06 +DEAL:: Quadrature point 0: [16.00 0 ][-27.71 0 ][35.78 0 ][-42.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 35.78 ][0 -42.33 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][16.00 0 ][-27.71 0 ][35.78 0 ][-42.33 0 ][0 16.00 ][0 27.71 ][0 35.78 ][0 42.33 ][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 ] +DEAL:: Quadrature point 2: [16.00 0 ][27.71 0 ][35.78 0 ][42.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 -27.71 ][0 35.78 ][0 -42.33 ][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 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][16.00 0 ][27.71 0 ][35.78 0 ][42.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 16.00 ][0 27.71 ][0 35.78 ][0 42.33 ][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 ] +DEAL:: h=0.03 +DEAL:: Quadrature point 0: [32.00 0 ][-55.43 0 ][71.55 0 ][-84.66 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 71.55 ][0 -84.66 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][32.00 0 ][-55.43 0 ][71.55 0 ][-84.66 0 ][0 32.00 ][0 55.43 ][0 71.55 ][0 84.66 ][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 ] +DEAL:: Quadrature point 2: [32.00 0 ][55.43 0 ][71.55 0 ][84.66 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 -55.43 ][0 71.55 ][0 -84.66 ][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 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][32.00 0 ][55.43 0 ][71.55 0 ][84.66 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 32.00 ][0 55.43 ][0 71.55 ][0 84.66 ][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 ] +DEAL:: h=0.02 +DEAL:: Quadrature point 0: [64.00 0 ][-110.85 0 ][143.11 0 ][-169.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 143.11 ][0 -169.33 ][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 ] +DEAL:: Quadrature point 1: [0 0 ][0 0 ][0 0 ][0 0 ][64.00 0 ][-110.85 0 ][143.11 0 ][-169.33 0 ][0 64.00 ][0 110.85 ][0 143.11 ][0 169.33 ][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 ] +DEAL:: Quadrature point 2: [64.00 0 ][110.85 0 ][143.11 0 ][169.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 -110.85 ][0 143.11 ][0 -169.33 ][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 ] +DEAL:: Quadrature point 3: [0 0 ][0 0 ][0 0 ][0 0 ][64.00 0 ][110.85 0 ][143.11 0 ][169.33 0 ][0 0 ][0 0 ][0 0 ][0 0 ][0 64.00 ][0 110.85 ][0 143.11 ][0 169.33 ][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 ]