Added test for QGaussOneOverR

git-svn-id: https://svn.dealii.org/trunk@18495 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/tests/codim_one/Makefile b/tests/codim_one/Makefile
index 6efeeab122290604efa7d3158b16e0657e739ce4..13656b4874356ff772e269abca107e1ef86e70a4 100644 (file)
--- a/tests/codim_one/Makefile
+++ b/tests/codim_one/Makefile
@@ -27,10 +27,11 @@ tests_x = grid* \
        surface* \
        data_out* \
        projection*\
-       integrate_log*\
+       integrate*\
        gradients*\
        tensor_matrix_conversion*\
        bem*\
+       bspline*
        
        
 

diff --git a/tests/codim_one/integrate_one_over_r.cc b/tests/codim_one/integrate_one_over_r.cc
new file mode 100644 (file)
index 0000000..48f2616
--- /dev/null
+++ b/tests/codim_one/integrate_one_over_r.cc
@@ -0,0 +1,263 @@
+//----------------------------  integrate_one_over_r.cc  ---------------------------
+//    $Id: testsuite.html 13373 2006-07-13 13:12:08Z heltai $
+//    Version: $Name$ 
+//
+//    Copyright (C) 2005 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.
+//
+//----------------------------  integrate_one_over_r.cc  ---------------------------
+
+
+// integrates the function *f(x,y)/R, where f(x,y) is a power of x and
+// y on the set [0,1]x[0,1]. dim = 2 only.
+
+#include "../tests.h"
+#include <fstream>
+#include <base/logstream.h>
+#include <base/utilities.h>
+
+// all include files needed for the program
+#include <base/quadrature_lib.h>
+
+
+#include <fstream>
+#include <string>
+
+#include <math.h>
+
+using namespace std;
+
+// The exact integrals were computed with 20 digits precision using Maple. 
+double exact_integral_one_over_r(const unsigned int i,
+                                const unsigned int j,
+                                const unsigned int vertex_index);
+
+ofstream logfile("integrate_one_over_r/output");
+
+int main()
+{
+    deallog.attach(logfile);
+    deallog.depth_console(0);
+    deallog<< std::fixed;
+    
+    deallog << endl 
+           << "Calculation of the integral of f(x,y)*1/R on [0,1]x[0,1]" << endl
+           << "for f(x,y) = x^i y^j, with i,j ranging from 0 to 5, and R being" << endl
+           << "the distance from (x,y) to four vertices of the square." << endl 
+           << endl;
+
+    for(unsigned int m=1; m<7; ++m) {
+       deallog << " =========Quadrature Order: " << m 
+               << " =============================== " << endl;
+       deallog << " ============================================================ " << endl;
+       for(unsigned int index=0; index<4; ++index) {
+           deallog << " ===============Vertex Index: " << index 
+                   << " ============================= " << endl;
+           QGaussOneOverR<2> quad(m, index);
+           for (unsigned int i=0; i<6; ++i) {
+               for (unsigned int j=0; j<6; ++j) {
+                  
+                   double exact_integral  = exact_integral_one_over_r(index, i,j);
+                   double approx_integral = 0;
+                   
+                   for (unsigned int q=0; q< quad.size(); ++q) {
+                       double x = quad.point(q)[0];
+                       double y = quad.point(q)[1];
+                       approx_integral += ( pow(x, (double)i) * 
+                                            pow(y, (double)j) * 
+                                            quad.weight(q) );
+                   }
+
+                   deallog << "f(x,y) = x^" << i
+                           << " y^" << j 
+                           << ", Error = " 
+                           << approx_integral - exact_integral << endl;
+               }
+           }
+       }
+    }
+}
+
+double exact_integral_one_over_r(const unsigned int vertex_index,
+                                const unsigned int i,
+                                const unsigned int j) {
+    Assert(vertex_index < 4, ExcInternalError());
+    Assert(i<6, ExcNotImplemented());
+    Assert(j<6, ExcNotImplemented());
+    
+    // The integrals are computed using the following maple snippet of
+    // code:
+    // 
+    //  sing_int := proc(index, N, M) 
+    //     if index = 0 then 
+    //        return int(int(x^N *y^M/sqrt(x^2+y^2), x=0.0..1.0), y=0.0..1.0);
+    //     elif index = 1 then
+    //        return int(int(x^N *y^M/sqrt((x-1)^2+y^2), x=0.0..1.0), y=0.0..1.0);
+    //     elif index = 2 then 
+    //        return int(int(x^N *y^M/sqrt(x^2+(y-1)^2), x=0.0..1.0), y=0.0..1.0);
+    //     elif index = 3 then 
+    //        return int(int((1-x)^N *(1-y)^M/sqrt(x^2+y^2), x=0.0..1.0), y=0.0..1.0);
+    //     end if;
+    //  end proc;
+    //  Digits := 20;
+    //  for i from 3 to 3 do
+    //     for n from 0 to 5 do 
+    //     for m from 0 to 5 do
+    //           C( v[i+1][n+1][m+1] = sing_int(i, n, m), resultname="a");
+    //        end do;
+    //     end do;
+    //  end do;
+
+    static double v[4][6][6] = {{{0}}};
+    if(v[0][0][0] == 0) {
+       v[0][0][0] = 0.17627471740390860505e1;
+       v[0][0][1] = 0.64779357469631903702e0;
+       v[0][0][2] = 0.38259785823210634567e0;
+       v[0][0][3] = 0.26915893322379450224e0;
+       v[0][0][4] = 0.20702239737104695572e0;
+       v[0][0][5] = 0.16800109713227567467e0;
+       v[0][1][0] = 0.64779357469631903702e0;
+       v[0][1][1] = 0.27614237491539669920e0;
+       v[0][1][2] = 0.17015838751246776515e0;
+       v[0][1][3] = 0.12189514164974600651e0;
+       v[0][1][4] = 0.94658660368131133694e-1;
+       v[0][1][5] = 0.77263794021029438797e-1;
+       v[0][2][0] = 0.38259785823210634567e0;
+       v[0][2][1] = 0.17015838751246776515e0;
+       v[0][2][2] = 0.10656799507071040471e0;
+       v[0][2][3] = 0.76947022258735165920e-1;
+       v[0][2][4] = 0.60022626787495395021e-1;
+       v[0][2][5] = 0.49131622931360879320e-1;
+       v[0][3][0] = 0.26915893322379450224e0;
+       v[0][3][1] = 0.12189514164974600651e0;
+       v[0][3][2] = 0.76947022258735165919e-1;
+       v[0][3][3] = 0.55789184535895709637e-1;
+       v[0][3][4] = 0.43625068213915842136e-1;
+       v[0][3][5] = 0.35766126849971778500e-1;
+       v[0][4][0] = 0.20702239737104695572e0;
+       v[0][4][1] = 0.94658660368131133694e-1;
+       v[0][4][2] = 0.60022626787495395021e-1;
+       v[0][4][3] = 0.43625068213915842137e-1;
+       v[0][4][4] = 0.34164088852375945192e-1;
+       v[0][4][5] = 0.28037139560980277614e-1;
+       v[0][5][0] = 0.16800109713227567467e0;
+       v[0][5][1] = 0.77263794021029438797e-1;
+       v[0][5][2] = 0.49131622931360879320e-1;
+       v[0][5][3] = 0.35766126849971778501e-1;
+       v[0][5][4] = 0.28037139560980277614e-1;
+       v[0][5][5] = 0.23024181049838367777e-1;
+       v[1][0][0] = 0.17627471740390860505e1;
+       v[1][0][1] = 0.64779357469631903702e0;
+       v[1][0][2] = 0.38259785823210634567e0;
+       v[1][0][3] = 0.26915893322379450224e0;
+       v[1][0][4] = 0.20702239737104695572e0;
+       v[1][0][5] = 0.16800109713227567467e0;
+       v[1][1][0] = 0.11149535993427670134e1;
+       v[1][1][1] = 0.37165119978092233782e0;
+       v[1][1][2] = 0.21243947071963858053e0;
+       v[1][1][3] = 0.14726379157404849573e0;
+       v[1][1][4] = 0.11236373700291582202e0;
+       v[1][1][5] = 0.90737303111246235871e-1;
+       v[1][2][0] = 0.84975788287855432210e0;
+       v[1][2][1] = 0.26566721237799340376e0;
+       v[1][2][2] = 0.14884907827788122009e0;
+       v[1][2][3] = 0.10231567218303765515e0;
+       v[1][2][4] = 0.77727703422280083352e-1;
+       v[1][2][5] = 0.62605132021577676395e-1;
+       v[1][3][0] = 0.69800109142265347423e0;
+       v[1][3][1] = 0.20794647083778622837e0;
+       v[1][3][2] = 0.11487965864809909847e0;
+       v[1][3][3] = 0.78525390514866270852e-1;
+       v[1][3][4] = 0.59489228415223897572e-1;
+       v[1][3][5] = 0.47838457013298217744e-1;
+       v[1][4][0] = 0.59754668912231692323e0;
+       v[1][4][1] = 0.17125249387868593878e0;
+       v[1][4][2] = 0.93606816359052444729e-1;
+       v[1][4][3] = 0.63728830247554475330e-1;
+       v[1][4][4] = 0.48187332620207367724e-1;
+       v[1][4][5] = 0.38708290797416359020e-1;
+       v[1][5][0] = 0.52527944036356840363e0;
+       v[1][5][1] = 0.14574366656617935708e0;
+       v[1][5][2] = 0.78997159795636003667e-1;
+       v[1][5][3] = 0.53620816423066464705e-1;
+       v[1][5][4] = 0.40487985967086264433e-1;
+       v[1][5][5] = 0.32498604596082509165e-1;
+       v[2][0][0] = 0.17627471740390860505e1;
+       v[2][0][1] = 0.11149535993427670134e1;
+       v[2][0][2] = 0.84975788287855432210e0;
+       v[2][0][3] = 0.69800109142265347419e0;
+       v[2][0][4] = 0.59754668912231692318e0;
+       v[2][0][5] = 0.52527944036356840362e0;
+       v[2][1][0] = 0.64779357469631903702e0;
+       v[2][1][1] = 0.37165119978092233782e0;
+       v[2][1][2] = 0.26566721237799340376e0;
+       v[2][1][3] = 0.20794647083778622835e0;
+       v[2][1][4] = 0.17125249387868593876e0;
+       v[2][1][5] = 0.14574366656617935708e0;
+       v[2][2][0] = 0.38259785823210634567e0;
+       v[2][2][1] = 0.21243947071963858053e0;
+       v[2][2][2] = 0.14884907827788122009e0;
+       v[2][2][3] = 0.11487965864809909845e0;
+       v[2][2][4] = 0.93606816359052444712e-1;
+       v[2][2][5] = 0.78997159795636003667e-1;
+       v[2][3][0] = 0.26915893322379450223e0;
+       v[2][3][1] = 0.14726379157404849572e0;
+       v[2][3][2] = 0.10231567218303765514e0;
+       v[2][3][3] = 0.78525390514866270835e-1;
+       v[2][3][4] = 0.63728830247554475311e-1;
+       v[2][3][5] = 0.53620816423066464702e-1;
+       v[2][4][0] = 0.20702239737104695572e0;
+       v[2][4][1] = 0.11236373700291582202e0;
+       v[2][4][2] = 0.77727703422280083352e-1;
+       v[2][4][3] = 0.59489228415223897563e-1;
+       v[2][4][4] = 0.48187332620207367713e-1;
+       v[2][4][5] = 0.40487985967086264434e-1;
+       v[2][5][0] = 0.16800109713227567468e0;
+       v[2][5][1] = 0.90737303111246235879e-1;
+       v[2][5][2] = 0.62605132021577676399e-1;
+       v[2][5][3] = 0.47838457013298217740e-1;
+       v[2][5][4] = 0.38708290797416359014e-1;
+       v[2][5][5] = 0.32498604596082509169e-1;
+       v[3][0][0] = 0.17627471740390860505e1;
+       v[3][0][1] = 0.11149535993427670134e1;
+       v[3][0][2] = 0.84975788287855432210e0;
+       v[3][0][3] = 0.69800109142265347419e0;
+       v[3][0][4] = 0.59754668912231692318e0;
+       v[3][0][5] = 0.52527944036356840362e0;
+       v[3][1][0] = 0.11149535993427670134e1;
+       v[3][1][1] = 0.74330239956184467563e0;
+       v[3][1][2] = 0.58409067050056091834e0;
+       v[3][1][3] = 0.49005462058486724584e0;
+       v[3][1][4] = 0.42629419524363098443e0;
+       v[3][1][5] = 0.37953577379738904654e0;
+       v[3][2][0] = 0.84975788287855432210e0;
+       v[3][2][1] = 0.58409067050056091834e0;
+       v[3][2][2] = 0.46727253640044873467e0;
+       v[3][2][3] = 0.39698780839518011595e0;
+       v[3][2][4] = 0.34864851772399749038e0;
+       v[3][2][5] = 0.31278926702684569312e0;
+       v[3][3][0] = 0.69800109142265347423e0;
+       v[3][3][1] = 0.49005462058486724586e0;
+       v[3][3][2] = 0.39698780839518011599e0;
+       v[3][3][3] = 0.34027526433872581371e0;
+       v[3][3][4] = 0.30088082631586196583e0;
+       v[3][3][5] = 0.27141910362887187844e0;
+       v[3][4][0] = 0.59754668912231692323e0;
+       v[3][4][1] = 0.42629419524363098445e0;
+       v[3][4][2] = 0.34864851772399749044e0;
+       v[3][4][3] = 0.30088082631586196576e0;
+       v[3][4][4] = 0.26744962339187730308e0;
+       v[3][4][5] = 0.24229245314748740295e0;
+       v[3][5][0] = 0.52527944036356840363e0;
+       v[3][5][1] = 0.37953577379738904655e0;
+       v[3][5][2] = 0.31278926702684569301e0;
+       v[3][5][3] = 0.27141910362887187862e0;
+       v[3][5][4] = 0.24229245314748740263e0;
+       v[3][5][5] = 0.22026586649771582089e0;
+    }
+    return v[vertex_index][i][j];
+}

diff --git a/tests/codim_one/integrate_one_over_r/cmp/generic b/tests/codim_one/integrate_one_over_r/cmp/generic
new file mode 100644 (file)
index 0000000..876b05f
--- /dev/null
+++ b/tests/codim_one/integrate_one_over_r/cmp/generic
@@ -0,0 +1,906 @@
+
+DEAL::
+DEAL::Calculation of the integral of f(x,y)*1/R on [0,1]x[0,1]
+DEAL::for f(x,y) = x^i y^j, with i,j ranging from 0 to 5, and R being
+DEAL::the distance from (x,y) to four vertices of the square.
+DEAL::
+DEAL:: =========Quadrature Order: 1 =============================== 
+DEAL:: ============================================================ 
+DEAL:: ===============Vertex Index: 0 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = -0.062529
+DEAL::f(x,y) = x^0 y^1, Error = -0.046676
+DEAL::f(x,y) = x^0 y^2, Error = -0.133607
+DEAL::f(x,y) = x^0 y^3, Error = -0.155343
+DEAL::f(x,y) = x^0 y^4, Error = -0.152327
+DEAL::f(x,y) = x^0 y^5, Error = -0.141111
+DEAL::f(x,y) = x^1 y^0, Error = -0.046676
+DEAL::f(x,y) = x^1 y^1, Error = -0.100079
+DEAL::f(x,y) = x^1 y^2, Error = -0.107911
+DEAL::f(x,y) = x^1 y^3, Error = -0.096111
+DEAL::f(x,y) = x^1 y^4, Error = -0.082873
+DEAL::f(x,y) = x^1 y^5, Error = -0.071600
+DEAL::f(x,y) = x^2 y^0, Error = -0.133607
+DEAL::f(x,y) = x^2 y^1, Error = -0.107911
+DEAL::f(x,y) = x^2 y^2, Error = -0.088336
+DEAL::f(x,y) = x^2 y^3, Error = -0.070501
+DEAL::f(x,y) = x^2 y^4, Error = -0.057353
+DEAL::f(x,y) = x^2 y^5, Error = -0.047911
+DEAL::f(x,y) = x^3 y^0, Error = -0.155343
+DEAL::f(x,y) = x^3 y^1, Error = -0.096111
+DEAL::f(x,y) = x^3 y^2, Error = -0.070501
+DEAL::f(x,y) = x^3 y^3, Error = -0.053901
+DEAL::f(x,y) = x^3 y^4, Error = -0.042958
+DEAL::f(x,y) = x^3 y^5, Error = -0.035490
+DEAL::f(x,y) = x^4 y^0, Error = -0.152327
+DEAL::f(x,y) = x^4 y^1, Error = -0.082873
+DEAL::f(x,y) = x^4 y^2, Error = -0.057353
+DEAL::f(x,y) = x^4 y^3, Error = -0.042958
+DEAL::f(x,y) = x^4 y^4, Error = -0.033969
+DEAL::f(x,y) = x^4 y^5, Error = -0.027968
+DEAL::f(x,y) = x^5 y^0, Error = -0.141111
+DEAL::f(x,y) = x^5 y^1, Error = -0.071600
+DEAL::f(x,y) = x^5 y^2, Error = -0.047911
+DEAL::f(x,y) = x^5 y^3, Error = -0.035490
+DEAL::f(x,y) = x^5 y^4, Error = -0.027968
+DEAL::f(x,y) = x^5 y^5, Error = -0.023004
+DEAL:: ===============Vertex Index: 1 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = -0.062529
+DEAL::f(x,y) = x^0 y^1, Error = -0.046676
+DEAL::f(x,y) = x^0 y^2, Error = -0.133607
+DEAL::f(x,y) = x^0 y^3, Error = -0.155343
+DEAL::f(x,y) = x^0 y^4, Error = -0.152327
+DEAL::f(x,y) = x^0 y^5, Error = -0.141111
+DEAL::f(x,y) = x^1 y^0, Error = -0.015854
+DEAL::f(x,y) = x^1 y^1, Error = 0.053403
+DEAL::f(x,y) = x^1 y^2, Error = -0.025696
+DEAL::f(x,y) = x^1 y^3, Error = -0.059232
+DEAL::f(x,y) = x^1 y^4, Error = -0.069454
+DEAL::f(x,y) = x^1 y^5, Error = -0.069511
+DEAL::f(x,y) = x^2 y^0, Error = -0.102785
+DEAL::f(x,y) = x^2 y^1, Error = 0.045572
+DEAL::f(x,y) = x^2 y^2, Error = 0.006122
+DEAL::f(x,y) = x^2 y^3, Error = -0.033622
+DEAL::f(x,y) = x^2 y^4, Error = -0.043934
+DEAL::f(x,y) = x^2 y^5, Error = -0.045823
+DEAL::f(x,y) = x^3 y^0, Error = -0.167979
+DEAL::f(x,y) = x^3 y^1, Error = 0.025941
+DEAL::f(x,y) = x^3 y^2, Error = 0.004382
+DEAL::f(x,y) = x^3 y^3, Error = -0.024612
+DEAL::f(x,y) = x^3 y^4, Error = -0.032809
+DEAL::f(x,y) = x^3 y^5, Error = -0.034556
+DEAL::f(x,y) = x^4 y^0, Error = -0.208419
+DEAL::f(x,y) = x^4 y^1, Error = 0.007749
+DEAL::f(x,y) = x^4 y^2, Error = 0.007329
+DEAL::f(x,y) = x^4 y^3, Error = -0.021257
+DEAL::f(x,y) = x^4 y^4, Error = -0.027090
+DEAL::f(x,y) = x^4 y^5, Error = -0.028188
+DEAL::f(x,y) = x^5 y^0, Error = -0.232305
+DEAL::f(x,y) = x^5 y^1, Error = 0.007037
+DEAL::f(x,y) = x^5 y^2, Error = -0.011256
+DEAL::f(x,y) = x^5 y^3, Error = -0.020084
+DEAL::f(x,y) = x^5 y^4, Error = -0.023789
+DEAL::f(x,y) = x^5 y^5, Error = -0.024163
+DEAL:: ===============Vertex Index: 2 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = -0.062529
+DEAL::f(x,y) = x^0 y^1, Error = -0.015854
+DEAL::f(x,y) = x^0 y^2, Error = -0.102785
+DEAL::f(x,y) = x^0 y^3, Error = -0.167979
+DEAL::f(x,y) = x^0 y^4, Error = -0.208419
+DEAL::f(x,y) = x^0 y^5, Error = -0.232305
+DEAL::f(x,y) = x^1 y^0, Error = -0.046676
+DEAL::f(x,y) = x^1 y^1, Error = 0.053403
+DEAL::f(x,y) = x^1 y^2, Error = 0.045572
+DEAL::f(x,y) = x^1 y^3, Error = 0.025941
+DEAL::f(x,y) = x^1 y^4, Error = 0.007749
+DEAL::f(x,y) = x^1 y^5, Error = 0.007037
+DEAL::f(x,y) = x^2 y^0, Error = -0.133607
+DEAL::f(x,y) = x^2 y^1, Error = -0.025696
+DEAL::f(x,y) = x^2 y^2, Error = 0.006122
+DEAL::f(x,y) = x^2 y^3, Error = 0.004382
+DEAL::f(x,y) = x^2 y^4, Error = 0.007329
+DEAL::f(x,y) = x^2 y^5, Error = -0.011256
+DEAL::f(x,y) = x^3 y^0, Error = -0.155343
+DEAL::f(x,y) = x^3 y^1, Error = -0.059232
+DEAL::f(x,y) = x^3 y^2, Error = -0.033622
+DEAL::f(x,y) = x^3 y^3, Error = -0.024612
+DEAL::f(x,y) = x^3 y^4, Error = -0.021257
+DEAL::f(x,y) = x^3 y^5, Error = -0.020084
+DEAL::f(x,y) = x^4 y^0, Error = -0.152327
+DEAL::f(x,y) = x^4 y^1, Error = -0.069454
+DEAL::f(x,y) = x^4 y^2, Error = -0.043934
+DEAL::f(x,y) = x^4 y^3, Error = -0.032809
+DEAL::f(x,y) = x^4 y^4, Error = -0.027090
+DEAL::f(x,y) = x^4 y^5, Error = -0.023789
+DEAL::f(x,y) = x^5 y^0, Error = -0.141111
+DEAL::f(x,y) = x^5 y^1, Error = -0.069511
+DEAL::f(x,y) = x^5 y^2, Error = -0.045823
+DEAL::f(x,y) = x^5 y^3, Error = -0.034556
+DEAL::f(x,y) = x^5 y^4, Error = -0.028188
+DEAL::f(x,y) = x^5 y^5, Error = -0.024163
+DEAL:: ===============Vertex Index: 3 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = -0.062529
+DEAL::f(x,y) = x^0 y^1, Error = -0.015854
+DEAL::f(x,y) = x^0 y^2, Error = -0.102785
+DEAL::f(x,y) = x^0 y^3, Error = -0.167979
+DEAL::f(x,y) = x^0 y^4, Error = -0.208419
+DEAL::f(x,y) = x^0 y^5, Error = -0.232305
+DEAL::f(x,y) = x^1 y^0, Error = -0.015854
+DEAL::f(x,y) = x^1 y^1, Error = -0.069257
+DEAL::f(x,y) = x^1 y^2, Error = -0.148356
+DEAL::f(x,y) = x^1 y^3, Error = -0.193920
+DEAL::f(x,y) = x^1 y^4, Error = -0.216169
+DEAL::f(x,y) = x^1 y^5, Error = -0.225268
+DEAL::f(x,y) = x^2 y^0, Error = -0.102785
+DEAL::f(x,y) = x^2 y^1, Error = -0.148356
+DEAL::f(x,y) = x^2 y^2, Error = -0.200049
+DEAL::f(x,y) = x^2 y^3, Error = -0.224242
+DEAL::f(x,y) = x^2 y^4, Error = -0.231247
+DEAL::f(x,y) = x^2 y^5, Error = -0.229486
+DEAL::f(x,y) = x^3 y^0, Error = -0.167979
+DEAL::f(x,y) = x^3 y^1, Error = -0.193920
+DEAL::f(x,y) = x^3 y^2, Error = -0.224242
+DEAL::f(x,y) = x^3 y^3, Error = -0.234336
+DEAL::f(x,y) = x^3 y^4, Error = -0.232396
+DEAL::f(x,y) = x^3 y^5, Error = -0.224876
+DEAL::f(x,y) = x^4 y^0, Error = -0.208419
+DEAL::f(x,y) = x^4 y^1, Error = -0.216169
+DEAL::f(x,y) = x^4 y^2, Error = -0.231247
+DEAL::f(x,y) = x^4 y^3, Error = -0.232396
+DEAL::f(x,y) = x^4 y^4, Error = -0.225450
+DEAL::f(x,y) = x^4 y^5, Error = -0.215142
+DEAL::f(x,y) = x^5 y^0, Error = -0.232305
+DEAL::f(x,y) = x^5 y^1, Error = -0.225268
+DEAL::f(x,y) = x^5 y^2, Error = -0.229486
+DEAL::f(x,y) = x^5 y^3, Error = -0.224876
+DEAL::f(x,y) = x^5 y^4, Error = -0.215142
+DEAL::f(x,y) = x^5 y^5, Error = -0.203615
+DEAL:: =========Quadrature Order: 2 =============================== 
+DEAL:: ============================================================ 
+DEAL:: ===============Vertex Index: 0 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.001795
+DEAL::f(x,y) = x^0 y^1, Error = 0.002254
+DEAL::f(x,y) = x^0 y^2, Error = 0.003643
+DEAL::f(x,y) = x^0 y^3, Error = 0.004860
+DEAL::f(x,y) = x^0 y^4, Error = -0.011516
+DEAL::f(x,y) = x^0 y^5, Error = -0.019894
+DEAL::f(x,y) = x^1 y^0, Error = 0.002254
+DEAL::f(x,y) = x^1 y^1, Error = 0.002407
+DEAL::f(x,y) = x^1 y^2, Error = 0.003410
+DEAL::f(x,y) = x^1 y^3, Error = 0.007693
+DEAL::f(x,y) = x^1 y^4, Error = -0.012832
+DEAL::f(x,y) = x^1 y^5, Error = -0.017132
+DEAL::f(x,y) = x^2 y^0, Error = 0.003643
+DEAL::f(x,y) = x^2 y^1, Error = 0.003410
+DEAL::f(x,y) = x^2 y^2, Error = 0.006861
+DEAL::f(x,y) = x^2 y^3, Error = -0.010777
+DEAL::f(x,y) = x^2 y^4, Error = -0.014116
+DEAL::f(x,y) = x^2 y^5, Error = -0.016356
+DEAL::f(x,y) = x^3 y^0, Error = 0.004860
+DEAL::f(x,y) = x^3 y^1, Error = 0.007693
+DEAL::f(x,y) = x^3 y^2, Error = -0.010777
+DEAL::f(x,y) = x^3 y^3, Error = -0.013245
+DEAL::f(x,y) = x^3 y^4, Error = -0.014991
+DEAL::f(x,y) = x^3 y^5, Error = -0.015882
+DEAL::f(x,y) = x^4 y^0, Error = -0.011516
+DEAL::f(x,y) = x^4 y^1, Error = -0.012832
+DEAL::f(x,y) = x^4 y^2, Error = -0.014116
+DEAL::f(x,y) = x^4 y^3, Error = -0.014991
+DEAL::f(x,y) = x^4 y^4, Error = -0.015454
+DEAL::f(x,y) = x^4 y^5, Error = -0.015409
+DEAL::f(x,y) = x^5 y^0, Error = -0.019894
+DEAL::f(x,y) = x^5 y^1, Error = -0.017132
+DEAL::f(x,y) = x^5 y^2, Error = -0.016356
+DEAL::f(x,y) = x^5 y^3, Error = -0.015882
+DEAL::f(x,y) = x^5 y^4, Error = -0.015409
+DEAL::f(x,y) = x^5 y^5, Error = -0.014742
+DEAL:: ===============Vertex Index: 1 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.001795
+DEAL::f(x,y) = x^0 y^1, Error = 0.002254
+DEAL::f(x,y) = x^0 y^2, Error = 0.003643
+DEAL::f(x,y) = x^0 y^3, Error = 0.004860
+DEAL::f(x,y) = x^0 y^4, Error = -0.011516
+DEAL::f(x,y) = x^0 y^5, Error = -0.019894
+DEAL::f(x,y) = x^1 y^0, Error = 0.000459
+DEAL::f(x,y) = x^1 y^1, Error = 0.000153
+DEAL::f(x,y) = x^1 y^2, Error = 0.000232
+DEAL::f(x,y) = x^1 y^3, Error = 0.002834
+DEAL::f(x,y) = x^1 y^4, Error = 0.001317
+DEAL::f(x,y) = x^1 y^5, Error = 0.002762
+DEAL::f(x,y) = x^2 y^0, Error = 0.000929
+DEAL::f(x,y) = x^2 y^1, Error = 0.000850
+DEAL::f(x,y) = x^2 y^2, Error = 0.003683
+DEAL::f(x,y) = x^2 y^3, Error = 0.000250
+DEAL::f(x,y) = x^2 y^4, Error = 0.000033
+DEAL::f(x,y) = x^2 y^5, Error = 0.001986
+DEAL::f(x,y) = x^3 y^0, Error = 0.001100
+DEAL::f(x,y) = x^3 y^1, Error = 0.002430
+DEAL::f(x,y) = x^3 y^2, Error = 0.003217
+DEAL::f(x,y) = x^3 y^3, Error = 0.000866
+DEAL::f(x,y) = x^3 y^4, Error = 0.000376
+DEAL::f(x,y) = x^3 y^5, Error = 0.001684
+DEAL::f(x,y) = x^4 y^0, Error = 0.006710
+DEAL::f(x,y) = x^4 y^1, Error = 0.004854
+DEAL::f(x,y) = x^4 y^2, Error = 0.002173
+DEAL::f(x,y) = x^4 y^3, Error = 0.000760
+DEAL::f(x,y) = x^4 y^4, Error = 0.000374
+DEAL::f(x,y) = x^4 y^5, Error = 0.001384
+DEAL::f(x,y) = x^5 y^0, Error = -0.016036
+DEAL::f(x,y) = x^5 y^1, Error = 0.005584
+DEAL::f(x,y) = x^5 y^2, Error = 0.001650
+DEAL::f(x,y) = x^5 y^3, Error = 0.000785
+DEAL::f(x,y) = x^5 y^4, Error = 0.000471
+DEAL::f(x,y) = x^5 y^5, Error = 0.001279
+DEAL:: ===============Vertex Index: 2 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.001795
+DEAL::f(x,y) = x^0 y^1, Error = 0.000459
+DEAL::f(x,y) = x^0 y^2, Error = 0.000929
+DEAL::f(x,y) = x^0 y^3, Error = 0.001100
+DEAL::f(x,y) = x^0 y^4, Error = 0.006710
+DEAL::f(x,y) = x^0 y^5, Error = -0.016036
+DEAL::f(x,y) = x^1 y^0, Error = 0.002254
+DEAL::f(x,y) = x^1 y^1, Error = 0.000153
+DEAL::f(x,y) = x^1 y^2, Error = 0.000850
+DEAL::f(x,y) = x^1 y^3, Error = 0.002430
+DEAL::f(x,y) = x^1 y^4, Error = 0.004854
+DEAL::f(x,y) = x^1 y^5, Error = 0.005584
+DEAL::f(x,y) = x^2 y^0, Error = 0.003643
+DEAL::f(x,y) = x^2 y^1, Error = 0.000232
+DEAL::f(x,y) = x^2 y^2, Error = 0.003683
+DEAL::f(x,y) = x^2 y^3, Error = 0.003217
+DEAL::f(x,y) = x^2 y^4, Error = 0.002173
+DEAL::f(x,y) = x^2 y^5, Error = 0.001650
+DEAL::f(x,y) = x^3 y^0, Error = 0.004860
+DEAL::f(x,y) = x^3 y^1, Error = 0.002834
+DEAL::f(x,y) = x^3 y^2, Error = 0.000250
+DEAL::f(x,y) = x^3 y^3, Error = 0.000866
+DEAL::f(x,y) = x^3 y^4, Error = 0.000760
+DEAL::f(x,y) = x^3 y^5, Error = 0.000785
+DEAL::f(x,y) = x^4 y^0, Error = -0.011516
+DEAL::f(x,y) = x^4 y^1, Error = 0.001317
+DEAL::f(x,y) = x^4 y^2, Error = 0.000033
+DEAL::f(x,y) = x^4 y^3, Error = 0.000376
+DEAL::f(x,y) = x^4 y^4, Error = 0.000374
+DEAL::f(x,y) = x^4 y^5, Error = 0.000471
+DEAL::f(x,y) = x^5 y^0, Error = -0.019894
+DEAL::f(x,y) = x^5 y^1, Error = 0.002762
+DEAL::f(x,y) = x^5 y^2, Error = 0.001986
+DEAL::f(x,y) = x^5 y^3, Error = 0.001684
+DEAL::f(x,y) = x^5 y^4, Error = 0.001384
+DEAL::f(x,y) = x^5 y^5, Error = 0.001279
+DEAL:: ===============Vertex Index: 3 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.001795
+DEAL::f(x,y) = x^0 y^1, Error = 0.000459
+DEAL::f(x,y) = x^0 y^2, Error = 0.000929
+DEAL::f(x,y) = x^0 y^3, Error = 0.001100
+DEAL::f(x,y) = x^0 y^4, Error = 0.006710
+DEAL::f(x,y) = x^0 y^5, Error = -0.016036
+DEAL::f(x,y) = x^1 y^0, Error = 0.000459
+DEAL::f(x,y) = x^1 y^1, Error = 0.000306
+DEAL::f(x,y) = x^1 y^2, Error = 0.000079
+DEAL::f(x,y) = x^1 y^3, Error = 0.003531
+DEAL::f(x,y) = x^1 y^4, Error = -0.011565
+DEAL::f(x,y) = x^1 y^5, Error = -0.021620
+DEAL::f(x,y) = x^2 y^0, Error = 0.000929
+DEAL::f(x,y) = x^2 y^1, Error = 0.000079
+DEAL::f(x,y) = x^2 y^2, Error = 0.002912
+DEAL::f(x,y) = x^2 y^3, Error = 0.009178
+DEAL::f(x,y) = x^2 y^4, Error = -0.018592
+DEAL::f(x,y) = x^2 y^5, Error = -0.028854
+DEAL::f(x,y) = x^3 y^0, Error = 0.001100
+DEAL::f(x,y) = x^3 y^1, Error = 0.003531
+DEAL::f(x,y) = x^3 y^2, Error = 0.009178
+DEAL::f(x,y) = x^3 y^3, Error = -0.017175
+DEAL::f(x,y) = x^3 y^4, Error = -0.027033
+DEAL::f(x,y) = x^3 y^5, Error = -0.036953
+DEAL::f(x,y) = x^4 y^0, Error = 0.006710
+DEAL::f(x,y) = x^4 y^1, Error = -0.011565
+DEAL::f(x,y) = x^4 y^2, Error = -0.018592
+DEAL::f(x,y) = x^4 y^3, Error = -0.027033
+DEAL::f(x,y) = x^4 y^4, Error = -0.036502
+DEAL::f(x,y) = x^4 y^5, Error = -0.045603
+DEAL::f(x,y) = x^5 y^0, Error = -0.016036
+DEAL::f(x,y) = x^5 y^1, Error = -0.021620
+DEAL::f(x,y) = x^5 y^2, Error = -0.028854
+DEAL::f(x,y) = x^5 y^3, Error = -0.036953
+DEAL::f(x,y) = x^5 y^4, Error = -0.045603
+DEAL::f(x,y) = x^5 y^5, Error = -0.053683
+DEAL:: =========Quadrature Order: 3 =============================== 
+DEAL:: ============================================================ 
+DEAL:: ===============Vertex Index: 0 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000052
+DEAL::f(x,y) = x^0 y^1, Error = 0.000092
+DEAL::f(x,y) = x^0 y^2, Error = 0.000200
+DEAL::f(x,y) = x^0 y^3, Error = 0.000369
+DEAL::f(x,y) = x^0 y^4, Error = 0.000603
+DEAL::f(x,y) = x^0 y^5, Error = 0.000866
+DEAL::f(x,y) = x^1 y^0, Error = 0.000092
+DEAL::f(x,y) = x^1 y^1, Error = 0.000105
+DEAL::f(x,y) = x^1 y^2, Error = 0.000183
+DEAL::f(x,y) = x^1 y^3, Error = 0.000322
+DEAL::f(x,y) = x^1 y^4, Error = 0.000525
+DEAL::f(x,y) = x^1 y^5, Error = 0.000953
+DEAL::f(x,y) = x^2 y^0, Error = 0.000200
+DEAL::f(x,y) = x^2 y^1, Error = 0.000183
+DEAL::f(x,y) = x^2 y^2, Error = 0.000229
+DEAL::f(x,y) = x^2 y^3, Error = 0.000337
+DEAL::f(x,y) = x^2 y^4, Error = 0.000658
+DEAL::f(x,y) = x^2 y^5, Error = 0.001202
+DEAL::f(x,y) = x^3 y^0, Error = 0.000369
+DEAL::f(x,y) = x^3 y^1, Error = 0.000322
+DEAL::f(x,y) = x^3 y^2, Error = 0.000337
+DEAL::f(x,y) = x^3 y^3, Error = 0.000553
+DEAL::f(x,y) = x^3 y^4, Error = 0.000986
+DEAL::f(x,y) = x^3 y^5, Error = 0.001572
+DEAL::f(x,y) = x^4 y^0, Error = 0.000603
+DEAL::f(x,y) = x^4 y^1, Error = 0.000525
+DEAL::f(x,y) = x^4 y^2, Error = 0.000658
+DEAL::f(x,y) = x^4 y^3, Error = 0.000986
+DEAL::f(x,y) = x^4 y^4, Error = 0.001464
+DEAL::f(x,y) = x^4 y^5, Error = 0.002024
+DEAL::f(x,y) = x^5 y^0, Error = 0.000866
+DEAL::f(x,y) = x^5 y^1, Error = 0.000953
+DEAL::f(x,y) = x^5 y^2, Error = 0.001202
+DEAL::f(x,y) = x^5 y^3, Error = 0.001572
+DEAL::f(x,y) = x^5 y^4, Error = 0.002024
+DEAL::f(x,y) = x^5 y^5, Error = 0.002508
+DEAL:: ===============Vertex Index: 1 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000052
+DEAL::f(x,y) = x^0 y^1, Error = 0.000092
+DEAL::f(x,y) = x^0 y^2, Error = 0.000200
+DEAL::f(x,y) = x^0 y^3, Error = 0.000369
+DEAL::f(x,y) = x^0 y^4, Error = 0.000603
+DEAL::f(x,y) = x^0 y^5, Error = 0.000866
+DEAL::f(x,y) = x^1 y^0, Error = 0.000040
+DEAL::f(x,y) = x^1 y^1, Error = 0.000013
+DEAL::f(x,y) = x^1 y^2, Error = 0.000017
+DEAL::f(x,y) = x^1 y^3, Error = 0.000047
+DEAL::f(x,y) = x^1 y^4, Error = 0.000079
+DEAL::f(x,y) = x^1 y^5, Error = 0.000086
+DEAL::f(x,y) = x^2 y^0, Error = 0.000068
+DEAL::f(x,y) = x^2 y^1, Error = 0.000064
+DEAL::f(x,y) = x^2 y^2, Error = 0.000063
+DEAL::f(x,y) = x^2 y^3, Error = 0.000063
+DEAL::f(x,y) = x^2 y^4, Error = 0.000212
+DEAL::f(x,y) = x^2 y^5, Error = 0.000163
+DEAL::f(x,y) = x^3 y^0, Error = 0.000005
+DEAL::f(x,y) = x^3 y^1, Error = 0.000002
+DEAL::f(x,y) = x^3 y^2, Error = 0.000001
+DEAL::f(x,y) = x^3 y^3, Error = 0.000137
+DEAL::f(x,y) = x^3 y^4, Error = 0.000017
+DEAL::f(x,y) = x^3 y^5, Error = 0.000043
+DEAL::f(x,y) = x^4 y^0, Error = 0.000007
+DEAL::f(x,y) = x^4 y^1, Error = 0.000004
+DEAL::f(x,y) = x^4 y^2, Error = 0.000152
+DEAL::f(x,y) = x^4 y^3, Error = 0.000121
+DEAL::f(x,y) = x^4 y^4, Error = 0.000028
+DEAL::f(x,y) = x^4 y^5, Error = 0.000004
+DEAL::f(x,y) = x^5 y^0, Error = 0.000044
+DEAL::f(x,y) = x^5 y^1, Error = 0.000156
+DEAL::f(x,y) = x^5 y^2, Error = 0.000290
+DEAL::f(x,y) = x^5 y^3, Error = 0.000040
+DEAL::f(x,y) = x^5 y^4, Error = 0.000004
+DEAL::f(x,y) = x^5 y^5, Error = 0.000015
+DEAL:: ===============Vertex Index: 2 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000052
+DEAL::f(x,y) = x^0 y^1, Error = 0.000040
+DEAL::f(x,y) = x^0 y^2, Error = 0.000068
+DEAL::f(x,y) = x^0 y^3, Error = 0.000005
+DEAL::f(x,y) = x^0 y^4, Error = 0.000007
+DEAL::f(x,y) = x^0 y^5, Error = 0.000044
+DEAL::f(x,y) = x^1 y^0, Error = 0.000092
+DEAL::f(x,y) = x^1 y^1, Error = 0.000013
+DEAL::f(x,y) = x^1 y^2, Error = 0.000064
+DEAL::f(x,y) = x^1 y^3, Error = 0.000002
+DEAL::f(x,y) = x^1 y^4, Error = 0.000004
+DEAL::f(x,y) = x^1 y^5, Error = 0.000156
+DEAL::f(x,y) = x^2 y^0, Error = 0.000200
+DEAL::f(x,y) = x^2 y^1, Error = 0.000017
+DEAL::f(x,y) = x^2 y^2, Error = 0.000063
+DEAL::f(x,y) = x^2 y^3, Error = 0.000001
+DEAL::f(x,y) = x^2 y^4, Error = 0.000152
+DEAL::f(x,y) = x^2 y^5, Error = 0.000290
+DEAL::f(x,y) = x^3 y^0, Error = 0.000369
+DEAL::f(x,y) = x^3 y^1, Error = 0.000047
+DEAL::f(x,y) = x^3 y^2, Error = 0.000063
+DEAL::f(x,y) = x^3 y^3, Error = 0.000137
+DEAL::f(x,y) = x^3 y^4, Error = 0.000121
+DEAL::f(x,y) = x^3 y^5, Error = 0.000040
+DEAL::f(x,y) = x^4 y^0, Error = 0.000603
+DEAL::f(x,y) = x^4 y^1, Error = 0.000079
+DEAL::f(x,y) = x^4 y^2, Error = 0.000212
+DEAL::f(x,y) = x^4 y^3, Error = 0.000017
+DEAL::f(x,y) = x^4 y^4, Error = 0.000028
+DEAL::f(x,y) = x^4 y^5, Error = 0.000004
+DEAL::f(x,y) = x^5 y^0, Error = 0.000866
+DEAL::f(x,y) = x^5 y^1, Error = 0.000086
+DEAL::f(x,y) = x^5 y^2, Error = 0.000163
+DEAL::f(x,y) = x^5 y^3, Error = 0.000043
+DEAL::f(x,y) = x^5 y^4, Error = 0.000004
+DEAL::f(x,y) = x^5 y^5, Error = 0.000015
+DEAL:: ===============Vertex Index: 3 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000052
+DEAL::f(x,y) = x^0 y^1, Error = 0.000040
+DEAL::f(x,y) = x^0 y^2, Error = 0.000068
+DEAL::f(x,y) = x^0 y^3, Error = 0.000005
+DEAL::f(x,y) = x^0 y^4, Error = 0.000007
+DEAL::f(x,y) = x^0 y^5, Error = 0.000044
+DEAL::f(x,y) = x^1 y^0, Error = 0.000040
+DEAL::f(x,y) = x^1 y^1, Error = 0.000027
+DEAL::f(x,y) = x^1 y^2, Error = 0.000004
+DEAL::f(x,y) = x^1 y^3, Error = 0.000003
+DEAL::f(x,y) = x^1 y^4, Error = 0.000003
+DEAL::f(x,y) = x^1 y^5, Error = 0.000200
+DEAL::f(x,y) = x^2 y^0, Error = 0.000068
+DEAL::f(x,y) = x^2 y^1, Error = 0.000004
+DEAL::f(x,y) = x^2 y^2, Error = 0.000003
+DEAL::f(x,y) = x^2 y^3, Error = 0.000003
+DEAL::f(x,y) = x^2 y^4, Error = 0.000151
+DEAL::f(x,y) = x^2 y^5, Error = 0.000646
+DEAL::f(x,y) = x^3 y^0, Error = 0.000005
+DEAL::f(x,y) = x^3 y^1, Error = 0.000003
+DEAL::f(x,y) = x^3 y^2, Error = 0.000003
+DEAL::f(x,y) = x^3 y^3, Error = 0.000141
+DEAL::f(x,y) = x^3 y^4, Error = 0.000571
+DEAL::f(x,y) = x^3 y^5, Error = 0.001423
+DEAL::f(x,y) = x^4 y^0, Error = 0.000007
+DEAL::f(x,y) = x^4 y^1, Error = 0.000003
+DEAL::f(x,y) = x^4 y^2, Error = 0.000151
+DEAL::f(x,y) = x^4 y^3, Error = 0.000571
+DEAL::f(x,y) = x^4 y^4, Error = 0.001356
+DEAL::f(x,y) = x^4 y^5, Error = 0.002566
+DEAL::f(x,y) = x^5 y^0, Error = 0.000044
+DEAL::f(x,y) = x^5 y^1, Error = 0.000200
+DEAL::f(x,y) = x^5 y^2, Error = 0.000646
+DEAL::f(x,y) = x^5 y^3, Error = 0.001423
+DEAL::f(x,y) = x^5 y^4, Error = 0.002566
+DEAL::f(x,y) = x^5 y^5, Error = 0.004094
+DEAL:: =========Quadrature Order: 4 =============================== 
+DEAL:: ============================================================ 
+DEAL:: ===============Vertex Index: 0 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000002
+DEAL::f(x,y) = x^0 y^1, Error = 0.000003
+DEAL::f(x,y) = x^0 y^2, Error = 0.000009
+DEAL::f(x,y) = x^0 y^3, Error = 0.000022
+DEAL::f(x,y) = x^0 y^4, Error = 0.000045
+DEAL::f(x,y) = x^0 y^5, Error = 0.000080
+DEAL::f(x,y) = x^1 y^0, Error = 0.000003
+DEAL::f(x,y) = x^1 y^1, Error = 0.000004
+DEAL::f(x,y) = x^1 y^2, Error = 0.000008
+DEAL::f(x,y) = x^1 y^3, Error = 0.000019
+DEAL::f(x,y) = x^1 y^4, Error = 0.000038
+DEAL::f(x,y) = x^1 y^5, Error = 0.000069
+DEAL::f(x,y) = x^2 y^0, Error = 0.000009
+DEAL::f(x,y) = x^2 y^1, Error = 0.000008
+DEAL::f(x,y) = x^2 y^2, Error = 0.000011
+DEAL::f(x,y) = x^2 y^3, Error = 0.000019
+DEAL::f(x,y) = x^2 y^4, Error = 0.000036
+DEAL::f(x,y) = x^2 y^5, Error = 0.000063
+DEAL::f(x,y) = x^3 y^0, Error = 0.000022
+DEAL::f(x,y) = x^3 y^1, Error = 0.000019
+DEAL::f(x,y) = x^3 y^2, Error = 0.000019
+DEAL::f(x,y) = x^3 y^3, Error = 0.000025
+DEAL::f(x,y) = x^3 y^4, Error = 0.000039
+DEAL::f(x,y) = x^3 y^5, Error = 0.000070
+DEAL::f(x,y) = x^4 y^0, Error = 0.000045
+DEAL::f(x,y) = x^4 y^1, Error = 0.000038
+DEAL::f(x,y) = x^4 y^2, Error = 0.000036
+DEAL::f(x,y) = x^4 y^3, Error = 0.000039
+DEAL::f(x,y) = x^4 y^4, Error = 0.000057
+DEAL::f(x,y) = x^4 y^5, Error = 0.000099
+DEAL::f(x,y) = x^5 y^0, Error = 0.000080
+DEAL::f(x,y) = x^5 y^1, Error = 0.000069
+DEAL::f(x,y) = x^5 y^2, Error = 0.000063
+DEAL::f(x,y) = x^5 y^3, Error = 0.000070
+DEAL::f(x,y) = x^5 y^4, Error = 0.000099
+DEAL::f(x,y) = x^5 y^5, Error = 0.000155
+DEAL:: ===============Vertex Index: 1 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000002
+DEAL::f(x,y) = x^0 y^1, Error = 0.000003
+DEAL::f(x,y) = x^0 y^2, Error = 0.000009
+DEAL::f(x,y) = x^0 y^3, Error = 0.000022
+DEAL::f(x,y) = x^0 y^4, Error = 0.000045
+DEAL::f(x,y) = x^0 y^5, Error = 0.000080
+DEAL::f(x,y) = x^1 y^0, Error = 0.000002
+DEAL::f(x,y) = x^1 y^1, Error = 0.000001
+DEAL::f(x,y) = x^1 y^2, Error = 0.000001
+DEAL::f(x,y) = x^1 y^3, Error = 0.000003
+DEAL::f(x,y) = x^1 y^4, Error = 0.000007
+DEAL::f(x,y) = x^1 y^5, Error = 0.000011
+DEAL::f(x,y) = x^2 y^0, Error = 0.000004
+DEAL::f(x,y) = x^2 y^1, Error = 0.000004
+DEAL::f(x,y) = x^2 y^2, Error = 0.000004
+DEAL::f(x,y) = x^2 y^3, Error = 0.000004
+DEAL::f(x,y) = x^2 y^4, Error = 0.000004
+DEAL::f(x,y) = x^2 y^5, Error = 0.000005
+DEAL::f(x,y) = x^3 y^0, Error = 0.000003
+DEAL::f(x,y) = x^3 y^1, Error = 0.000002
+DEAL::f(x,y) = x^3 y^2, Error = 0.000002
+DEAL::f(x,y) = x^3 y^3, Error = 0.000002
+DEAL::f(x,y) = x^3 y^4, Error = 0.000001
+DEAL::f(x,y) = x^3 y^5, Error = 0.000008
+DEAL::f(x,y) = x^4 y^0, Error = 0.000001
+DEAL::f(x,y) = x^4 y^1, Error = 0.000001
+DEAL::f(x,y) = x^4 y^2, Error = 0.000001
+DEAL::f(x,y) = x^4 y^3, Error = 0.000001
+DEAL::f(x,y) = x^4 y^4, Error = 0.000008
+DEAL::f(x,y) = x^4 y^5, Error = 0.000000
+DEAL::f(x,y) = x^5 y^0, Error = 0.000002
+DEAL::f(x,y) = x^5 y^1, Error = 0.000002
+DEAL::f(x,y) = x^5 y^2, Error = 0.000002
+DEAL::f(x,y) = x^5 y^3, Error = 0.000006
+DEAL::f(x,y) = x^5 y^4, Error = 0.000008
+DEAL::f(x,y) = x^5 y^5, Error = 0.000004
+DEAL:: ===============Vertex Index: 2 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000002
+DEAL::f(x,y) = x^0 y^1, Error = 0.000002
+DEAL::f(x,y) = x^0 y^2, Error = 0.000004
+DEAL::f(x,y) = x^0 y^3, Error = 0.000003
+DEAL::f(x,y) = x^0 y^4, Error = 0.000001
+DEAL::f(x,y) = x^0 y^5, Error = 0.000002
+DEAL::f(x,y) = x^1 y^0, Error = 0.000003
+DEAL::f(x,y) = x^1 y^1, Error = 0.000001
+DEAL::f(x,y) = x^1 y^2, Error = 0.000004
+DEAL::f(x,y) = x^1 y^3, Error = 0.000002
+DEAL::f(x,y) = x^1 y^4, Error = 0.000001
+DEAL::f(x,y) = x^1 y^5, Error = 0.000002
+DEAL::f(x,y) = x^2 y^0, Error = 0.000009
+DEAL::f(x,y) = x^2 y^1, Error = 0.000001
+DEAL::f(x,y) = x^2 y^2, Error = 0.000004
+DEAL::f(x,y) = x^2 y^3, Error = 0.000002
+DEAL::f(x,y) = x^2 y^4, Error = 0.000001
+DEAL::f(x,y) = x^2 y^5, Error = 0.000002
+DEAL::f(x,y) = x^3 y^0, Error = 0.000022
+DEAL::f(x,y) = x^3 y^1, Error = 0.000003
+DEAL::f(x,y) = x^3 y^2, Error = 0.000004
+DEAL::f(x,y) = x^3 y^3, Error = 0.000002
+DEAL::f(x,y) = x^3 y^4, Error = 0.000001
+DEAL::f(x,y) = x^3 y^5, Error = 0.000006
+DEAL::f(x,y) = x^4 y^0, Error = 0.000045
+DEAL::f(x,y) = x^4 y^1, Error = 0.000007
+DEAL::f(x,y) = x^4 y^2, Error = 0.000004
+DEAL::f(x,y) = x^4 y^3, Error = 0.000001
+DEAL::f(x,y) = x^4 y^4, Error = 0.000008
+DEAL::f(x,y) = x^4 y^5, Error = 0.000008
+DEAL::f(x,y) = x^5 y^0, Error = 0.000080
+DEAL::f(x,y) = x^5 y^1, Error = 0.000011
+DEAL::f(x,y) = x^5 y^2, Error = 0.000005
+DEAL::f(x,y) = x^5 y^3, Error = 0.000008
+DEAL::f(x,y) = x^5 y^4, Error = 0.000000
+DEAL::f(x,y) = x^5 y^5, Error = 0.000004
+DEAL:: ===============Vertex Index: 3 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000002
+DEAL::f(x,y) = x^0 y^1, Error = 0.000002
+DEAL::f(x,y) = x^0 y^2, Error = 0.000004
+DEAL::f(x,y) = x^0 y^3, Error = 0.000003
+DEAL::f(x,y) = x^0 y^4, Error = 0.000001
+DEAL::f(x,y) = x^0 y^5, Error = 0.000002
+DEAL::f(x,y) = x^1 y^0, Error = 0.000002
+DEAL::f(x,y) = x^1 y^1, Error = 0.000001
+DEAL::f(x,y) = x^1 y^2, Error = 0.000000
+DEAL::f(x,y) = x^1 y^3, Error = 0.000000
+DEAL::f(x,y) = x^1 y^4, Error = 0.000000
+DEAL::f(x,y) = x^1 y^5, Error = 0.000000
+DEAL::f(x,y) = x^2 y^0, Error = 0.000004
+DEAL::f(x,y) = x^2 y^1, Error = 0.000000
+DEAL::f(x,y) = x^2 y^2, Error = 0.000000
+DEAL::f(x,y) = x^2 y^3, Error = 0.000000
+DEAL::f(x,y) = x^2 y^4, Error = 0.000000
+DEAL::f(x,y) = x^2 y^5, Error = 0.000000
+DEAL::f(x,y) = x^3 y^0, Error = 0.000003
+DEAL::f(x,y) = x^3 y^1, Error = 0.000000
+DEAL::f(x,y) = x^3 y^2, Error = 0.000000
+DEAL::f(x,y) = x^3 y^3, Error = 0.000000
+DEAL::f(x,y) = x^3 y^4, Error = 0.000000
+DEAL::f(x,y) = x^3 y^5, Error = 0.000007
+DEAL::f(x,y) = x^4 y^0, Error = 0.000001
+DEAL::f(x,y) = x^4 y^1, Error = 0.000000
+DEAL::f(x,y) = x^4 y^2, Error = 0.000000
+DEAL::f(x,y) = x^4 y^3, Error = 0.000000
+DEAL::f(x,y) = x^4 y^4, Error = 0.000007
+DEAL::f(x,y) = x^4 y^5, Error = 0.000035
+DEAL::f(x,y) = x^5 y^0, Error = 0.000002
+DEAL::f(x,y) = x^5 y^1, Error = 0.000000
+DEAL::f(x,y) = x^5 y^2, Error = 0.000000
+DEAL::f(x,y) = x^5 y^3, Error = 0.000007
+DEAL::f(x,y) = x^5 y^4, Error = 0.000035
+DEAL::f(x,y) = x^5 y^5, Error = 0.000102
+DEAL:: =========Quadrature Order: 5 =============================== 
+DEAL:: ============================================================ 
+DEAL:: ===============Vertex Index: 0 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000000
+DEAL::f(x,y) = x^0 y^1, Error = 0.000000
+DEAL::f(x,y) = x^0 y^2, Error = 0.000000
+DEAL::f(x,y) = x^0 y^3, Error = 0.000001
+DEAL::f(x,y) = x^0 y^4, Error = 0.000003
+DEAL::f(x,y) = x^0 y^5, Error = 0.000006
+DEAL::f(x,y) = x^1 y^0, Error = 0.000000
+DEAL::f(x,y) = x^1 y^1, Error = 0.000000
+DEAL::f(x,y) = x^1 y^2, Error = 0.000000
+DEAL::f(x,y) = x^1 y^3, Error = 0.000001
+DEAL::f(x,y) = x^1 y^4, Error = 0.000002
+DEAL::f(x,y) = x^1 y^5, Error = 0.000005
+DEAL::f(x,y) = x^2 y^0, Error = 0.000000
+DEAL::f(x,y) = x^2 y^1, Error = 0.000000
+DEAL::f(x,y) = x^2 y^2, Error = 0.000000
+DEAL::f(x,y) = x^2 y^3, Error = 0.000001
+DEAL::f(x,y) = x^2 y^4, Error = 0.000002
+DEAL::f(x,y) = x^2 y^5, Error = 0.000005
+DEAL::f(x,y) = x^3 y^0, Error = 0.000001
+DEAL::f(x,y) = x^3 y^1, Error = 0.000001
+DEAL::f(x,y) = x^3 y^2, Error = 0.000001
+DEAL::f(x,y) = x^3 y^3, Error = 0.000001
+DEAL::f(x,y) = x^3 y^4, Error = 0.000002
+DEAL::f(x,y) = x^3 y^5, Error = 0.000004
+DEAL::f(x,y) = x^4 y^0, Error = 0.000003
+DEAL::f(x,y) = x^4 y^1, Error = 0.000002
+DEAL::f(x,y) = x^4 y^2, Error = 0.000002
+DEAL::f(x,y) = x^4 y^3, Error = 0.000002
+DEAL::f(x,y) = x^4 y^4, Error = 0.000003
+DEAL::f(x,y) = x^4 y^5, Error = 0.000005
+DEAL::f(x,y) = x^5 y^0, Error = 0.000006
+DEAL::f(x,y) = x^5 y^1, Error = 0.000005
+DEAL::f(x,y) = x^5 y^2, Error = 0.000005
+DEAL::f(x,y) = x^5 y^3, Error = 0.000004
+DEAL::f(x,y) = x^5 y^4, Error = 0.000005
+DEAL::f(x,y) = x^5 y^5, Error = 0.000007
+DEAL:: ===============Vertex Index: 1 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000000
+DEAL::f(x,y) = x^0 y^1, Error = 0.000000
+DEAL::f(x,y) = x^0 y^2, Error = 0.000000
+DEAL::f(x,y) = x^0 y^3, Error = 0.000001
+DEAL::f(x,y) = x^0 y^4, Error = 0.000003
+DEAL::f(x,y) = x^0 y^5, Error = 0.000006
+DEAL::f(x,y) = x^1 y^0, Error = 0.000000
+DEAL::f(x,y) = x^1 y^1, Error = 0.000000
+DEAL::f(x,y) = x^1 y^2, Error = 0.000000
+DEAL::f(x,y) = x^1 y^3, Error = 0.000000
+DEAL::f(x,y) = x^1 y^4, Error = 0.000000
+DEAL::f(x,y) = x^1 y^5, Error = 0.000001
+DEAL::f(x,y) = x^2 y^0, Error = 0.000000
+DEAL::f(x,y) = x^2 y^1, Error = 0.000000
+DEAL::f(x,y) = x^2 y^2, Error = 0.000000
+DEAL::f(x,y) = x^2 y^3, Error = 0.000000
+DEAL::f(x,y) = x^2 y^4, Error = 0.000000
+DEAL::f(x,y) = x^2 y^5, Error = 0.000000
+DEAL::f(x,y) = x^3 y^0, Error = 0.000000
+DEAL::f(x,y) = x^3 y^1, Error = 0.000000
+DEAL::f(x,y) = x^3 y^2, Error = 0.000000
+DEAL::f(x,y) = x^3 y^3, Error = 0.000000
+DEAL::f(x,y) = x^3 y^4, Error = 0.000000
+DEAL::f(x,y) = x^3 y^5, Error = 0.000000
+DEAL::f(x,y) = x^4 y^0, Error = 0.000000
+DEAL::f(x,y) = x^4 y^1, Error = 0.000000
+DEAL::f(x,y) = x^4 y^2, Error = 0.000000
+DEAL::f(x,y) = x^4 y^3, Error = 0.000000
+DEAL::f(x,y) = x^4 y^4, Error = 0.000000
+DEAL::f(x,y) = x^4 y^5, Error = 0.000000
+DEAL::f(x,y) = x^5 y^0, Error = 0.000000
+DEAL::f(x,y) = x^5 y^1, Error = 0.000000
+DEAL::f(x,y) = x^5 y^2, Error = 0.000000
+DEAL::f(x,y) = x^5 y^3, Error = 0.000000
+DEAL::f(x,y) = x^5 y^4, Error = 0.000000
+DEAL::f(x,y) = x^5 y^5, Error = 0.000000
+DEAL:: ===============Vertex Index: 2 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000000
+DEAL::f(x,y) = x^0 y^1, Error = 0.000000
+DEAL::f(x,y) = x^0 y^2, Error = 0.000000
+DEAL::f(x,y) = x^0 y^3, Error = 0.000000
+DEAL::f(x,y) = x^0 y^4, Error = 0.000000
+DEAL::f(x,y) = x^0 y^5, Error = 0.000000
+DEAL::f(x,y) = x^1 y^0, Error = 0.000000
+DEAL::f(x,y) = x^1 y^1, Error = 0.000000
+DEAL::f(x,y) = x^1 y^2, Error = 0.000000
+DEAL::f(x,y) = x^1 y^3, Error = 0.000000
+DEAL::f(x,y) = x^1 y^4, Error = 0.000000
+DEAL::f(x,y) = x^1 y^5, Error = 0.000000
+DEAL::f(x,y) = x^2 y^0, Error = 0.000000
+DEAL::f(x,y) = x^2 y^1, Error = 0.000000
+DEAL::f(x,y) = x^2 y^2, Error = 0.000000
+DEAL::f(x,y) = x^2 y^3, Error = 0.000000
+DEAL::f(x,y) = x^2 y^4, Error = 0.000000
+DEAL::f(x,y) = x^2 y^5, Error = 0.000000
+DEAL::f(x,y) = x^3 y^0, Error = 0.000001
+DEAL::f(x,y) = x^3 y^1, Error = 0.000000
+DEAL::f(x,y) = x^3 y^2, Error = 0.000000
+DEAL::f(x,y) = x^3 y^3, Error = 0.000000
+DEAL::f(x,y) = x^3 y^4, Error = 0.000000
+DEAL::f(x,y) = x^3 y^5, Error = 0.000000
+DEAL::f(x,y) = x^4 y^0, Error = 0.000003
+DEAL::f(x,y) = x^4 y^1, Error = 0.000000
+DEAL::f(x,y) = x^4 y^2, Error = 0.000000
+DEAL::f(x,y) = x^4 y^3, Error = 0.000000
+DEAL::f(x,y) = x^4 y^4, Error = 0.000000
+DEAL::f(x,y) = x^4 y^5, Error = 0.000000
+DEAL::f(x,y) = x^5 y^0, Error = 0.000006
+DEAL::f(x,y) = x^5 y^1, Error = 0.000001
+DEAL::f(x,y) = x^5 y^2, Error = 0.000000
+DEAL::f(x,y) = x^5 y^3, Error = 0.000000
+DEAL::f(x,y) = x^5 y^4, Error = 0.000000
+DEAL::f(x,y) = x^5 y^5, Error = 0.000000
+DEAL:: ===============Vertex Index: 3 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000000
+DEAL::f(x,y) = x^0 y^1, Error = 0.000000
+DEAL::f(x,y) = x^0 y^2, Error = 0.000000
+DEAL::f(x,y) = x^0 y^3, Error = 0.000000
+DEAL::f(x,y) = x^0 y^4, Error = 0.000000
+DEAL::f(x,y) = x^0 y^5, Error = 0.000000
+DEAL::f(x,y) = x^1 y^0, Error = 0.000000
+DEAL::f(x,y) = x^1 y^1, Error = 0.000000
+DEAL::f(x,y) = x^1 y^2, Error = 0.000000
+DEAL::f(x,y) = x^1 y^3, Error = 0.000000
+DEAL::f(x,y) = x^1 y^4, Error = 0.000000
+DEAL::f(x,y) = x^1 y^5, Error = 0.000000
+DEAL::f(x,y) = x^2 y^0, Error = 0.000000
+DEAL::f(x,y) = x^2 y^1, Error = 0.000000
+DEAL::f(x,y) = x^2 y^2, Error = 0.000000
+DEAL::f(x,y) = x^2 y^3, Error = 0.000000
+DEAL::f(x,y) = x^2 y^4, Error = 0.000000
+DEAL::f(x,y) = x^2 y^5, Error = 0.000000
+DEAL::f(x,y) = x^3 y^0, Error = 0.000000
+DEAL::f(x,y) = x^3 y^1, Error = 0.000000
+DEAL::f(x,y) = x^3 y^2, Error = 0.000000
+DEAL::f(x,y) = x^3 y^3, Error = 0.000000
+DEAL::f(x,y) = x^3 y^4, Error = 0.000000
+DEAL::f(x,y) = x^3 y^5, Error = 0.000000
+DEAL::f(x,y) = x^4 y^0, Error = 0.000000
+DEAL::f(x,y) = x^4 y^1, Error = 0.000000
+DEAL::f(x,y) = x^4 y^2, Error = 0.000000
+DEAL::f(x,y) = x^4 y^3, Error = 0.000000
+DEAL::f(x,y) = x^4 y^4, Error = 0.000000
+DEAL::f(x,y) = x^4 y^5, Error = 0.000000
+DEAL::f(x,y) = x^5 y^0, Error = 0.000000
+DEAL::f(x,y) = x^5 y^1, Error = 0.000000
+DEAL::f(x,y) = x^5 y^2, Error = 0.000000
+DEAL::f(x,y) = x^5 y^3, Error = 0.000000
+DEAL::f(x,y) = x^5 y^4, Error = 0.000000
+DEAL::f(x,y) = x^5 y^5, Error = 0.000000
+DEAL:: =========Quadrature Order: 6 =============================== 
+DEAL:: ============================================================ 
+DEAL:: ===============Vertex Index: 0 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000000
+DEAL::f(x,y) = x^0 y^1, Error = 0.000000
+DEAL::f(x,y) = x^0 y^2, Error = 0.000000
+DEAL::f(x,y) = x^0 y^3, Error = 0.000000
+DEAL::f(x,y) = x^0 y^4, Error = 0.000000
+DEAL::f(x,y) = x^0 y^5, Error = 0.000000
+DEAL::f(x,y) = x^1 y^0, Error = 0.000000
+DEAL::f(x,y) = x^1 y^1, Error = 0.000000
+DEAL::f(x,y) = x^1 y^2, Error = 0.000000
+DEAL::f(x,y) = x^1 y^3, Error = 0.000000
+DEAL::f(x,y) = x^1 y^4, Error = 0.000000
+DEAL::f(x,y) = x^1 y^5, Error = 0.000000
+DEAL::f(x,y) = x^2 y^0, Error = 0.000000
+DEAL::f(x,y) = x^2 y^1, Error = 0.000000
+DEAL::f(x,y) = x^2 y^2, Error = 0.000000
+DEAL::f(x,y) = x^2 y^3, Error = 0.000000
+DEAL::f(x,y) = x^2 y^4, Error = 0.000000
+DEAL::f(x,y) = x^2 y^5, Error = 0.000000
+DEAL::f(x,y) = x^3 y^0, Error = 0.000000
+DEAL::f(x,y) = x^3 y^1, Error = 0.000000
+DEAL::f(x,y) = x^3 y^2, Error = 0.000000
+DEAL::f(x,y) = x^3 y^3, Error = 0.000000
+DEAL::f(x,y) = x^3 y^4, Error = 0.000000
+DEAL::f(x,y) = x^3 y^5, Error = 0.000000
+DEAL::f(x,y) = x^4 y^0, Error = 0.000000
+DEAL::f(x,y) = x^4 y^1, Error = 0.000000
+DEAL::f(x,y) = x^4 y^2, Error = 0.000000
+DEAL::f(x,y) = x^4 y^3, Error = 0.000000
+DEAL::f(x,y) = x^4 y^4, Error = 0.000000
+DEAL::f(x,y) = x^4 y^5, Error = 0.000000
+DEAL::f(x,y) = x^5 y^0, Error = 0.000000
+DEAL::f(x,y) = x^5 y^1, Error = 0.000000
+DEAL::f(x,y) = x^5 y^2, Error = 0.000000
+DEAL::f(x,y) = x^5 y^3, Error = 0.000000
+DEAL::f(x,y) = x^5 y^4, Error = 0.000000
+DEAL::f(x,y) = x^5 y^5, Error = 0.000000
+DEAL:: ===============Vertex Index: 1 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000000
+DEAL::f(x,y) = x^0 y^1, Error = 0.000000
+DEAL::f(x,y) = x^0 y^2, Error = 0.000000
+DEAL::f(x,y) = x^0 y^3, Error = 0.000000
+DEAL::f(x,y) = x^0 y^4, Error = 0.000000
+DEAL::f(x,y) = x^0 y^5, Error = 0.000000
+DEAL::f(x,y) = x^1 y^0, Error = 0.000000
+DEAL::f(x,y) = x^1 y^1, Error = 0.000000
+DEAL::f(x,y) = x^1 y^2, Error = 0.000000
+DEAL::f(x,y) = x^1 y^3, Error = 0.000000
+DEAL::f(x,y) = x^1 y^4, Error = 0.000000
+DEAL::f(x,y) = x^1 y^5, Error = 0.000000
+DEAL::f(x,y) = x^2 y^0, Error = 0.000000
+DEAL::f(x,y) = x^2 y^1, Error = 0.000000
+DEAL::f(x,y) = x^2 y^2, Error = 0.000000
+DEAL::f(x,y) = x^2 y^3, Error = 0.000000
+DEAL::f(x,y) = x^2 y^4, Error = 0.000000
+DEAL::f(x,y) = x^2 y^5, Error = 0.000000
+DEAL::f(x,y) = x^3 y^0, Error = 0.000000
+DEAL::f(x,y) = x^3 y^1, Error = 0.000000
+DEAL::f(x,y) = x^3 y^2, Error = 0.000000
+DEAL::f(x,y) = x^3 y^3, Error = 0.000000
+DEAL::f(x,y) = x^3 y^4, Error = 0.000000
+DEAL::f(x,y) = x^3 y^5, Error = 0.000000
+DEAL::f(x,y) = x^4 y^0, Error = 0.000000
+DEAL::f(x,y) = x^4 y^1, Error = 0.000000
+DEAL::f(x,y) = x^4 y^2, Error = 0.000000
+DEAL::f(x,y) = x^4 y^3, Error = 0.000000
+DEAL::f(x,y) = x^4 y^4, Error = 0.000000
+DEAL::f(x,y) = x^4 y^5, Error = 0.000000
+DEAL::f(x,y) = x^5 y^0, Error = 0.000000
+DEAL::f(x,y) = x^5 y^1, Error = 0.000000
+DEAL::f(x,y) = x^5 y^2, Error = 0.000000
+DEAL::f(x,y) = x^5 y^3, Error = 0.000000
+DEAL::f(x,y) = x^5 y^4, Error = 0.000000
+DEAL::f(x,y) = x^5 y^5, Error = 0.000000
+DEAL:: ===============Vertex Index: 2 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000000
+DEAL::f(x,y) = x^0 y^1, Error = 0.000000
+DEAL::f(x,y) = x^0 y^2, Error = 0.000000
+DEAL::f(x,y) = x^0 y^3, Error = 0.000000
+DEAL::f(x,y) = x^0 y^4, Error = 0.000000
+DEAL::f(x,y) = x^0 y^5, Error = 0.000000
+DEAL::f(x,y) = x^1 y^0, Error = 0.000000
+DEAL::f(x,y) = x^1 y^1, Error = 0.000000
+DEAL::f(x,y) = x^1 y^2, Error = 0.000000
+DEAL::f(x,y) = x^1 y^3, Error = 0.000000
+DEAL::f(x,y) = x^1 y^4, Error = 0.000000
+DEAL::f(x,y) = x^1 y^5, Error = 0.000000
+DEAL::f(x,y) = x^2 y^0, Error = 0.000000
+DEAL::f(x,y) = x^2 y^1, Error = 0.000000
+DEAL::f(x,y) = x^2 y^2, Error = 0.000000
+DEAL::f(x,y) = x^2 y^3, Error = 0.000000
+DEAL::f(x,y) = x^2 y^4, Error = 0.000000
+DEAL::f(x,y) = x^2 y^5, Error = 0.000000
+DEAL::f(x,y) = x^3 y^0, Error = 0.000000
+DEAL::f(x,y) = x^3 y^1, Error = 0.000000
+DEAL::f(x,y) = x^3 y^2, Error = 0.000000
+DEAL::f(x,y) = x^3 y^3, Error = 0.000000
+DEAL::f(x,y) = x^3 y^4, Error = 0.000000
+DEAL::f(x,y) = x^3 y^5, Error = 0.000000
+DEAL::f(x,y) = x^4 y^0, Error = 0.000000
+DEAL::f(x,y) = x^4 y^1, Error = 0.000000
+DEAL::f(x,y) = x^4 y^2, Error = 0.000000
+DEAL::f(x,y) = x^4 y^3, Error = 0.000000
+DEAL::f(x,y) = x^4 y^4, Error = 0.000000
+DEAL::f(x,y) = x^4 y^5, Error = 0.000000
+DEAL::f(x,y) = x^5 y^0, Error = 0.000000
+DEAL::f(x,y) = x^5 y^1, Error = 0.000000
+DEAL::f(x,y) = x^5 y^2, Error = 0.000000
+DEAL::f(x,y) = x^5 y^3, Error = 0.000000
+DEAL::f(x,y) = x^5 y^4, Error = 0.000000
+DEAL::f(x,y) = x^5 y^5, Error = 0.000000
+DEAL:: ===============Vertex Index: 3 ============================= 
+DEAL::f(x,y) = x^0 y^0, Error = 0.000000
+DEAL::f(x,y) = x^0 y^1, Error = 0.000000
+DEAL::f(x,y) = x^0 y^2, Error = 0.000000
+DEAL::f(x,y) = x^0 y^3, Error = 0.000000
+DEAL::f(x,y) = x^0 y^4, Error = 0.000000
+DEAL::f(x,y) = x^0 y^5, Error = 0.000000
+DEAL::f(x,y) = x^1 y^0, Error = 0.000000
+DEAL::f(x,y) = x^1 y^1, Error = 0.000000
+DEAL::f(x,y) = x^1 y^2, Error = 0.000000
+DEAL::f(x,y) = x^1 y^3, Error = 0.000000
+DEAL::f(x,y) = x^1 y^4, Error = 0.000000
+DEAL::f(x,y) = x^1 y^5, Error = 0.000000
+DEAL::f(x,y) = x^2 y^0, Error = 0.000000
+DEAL::f(x,y) = x^2 y^1, Error = 0.000000
+DEAL::f(x,y) = x^2 y^2, Error = 0.000000
+DEAL::f(x,y) = x^2 y^3, Error = 0.000000
+DEAL::f(x,y) = x^2 y^4, Error = 0.000000
+DEAL::f(x,y) = x^2 y^5, Error = 0.000000
+DEAL::f(x,y) = x^3 y^0, Error = 0.000000
+DEAL::f(x,y) = x^3 y^1, Error = 0.000000
+DEAL::f(x,y) = x^3 y^2, Error = 0.000000
+DEAL::f(x,y) = x^3 y^3, Error = 0.000000
+DEAL::f(x,y) = x^3 y^4, Error = 0.000000
+DEAL::f(x,y) = x^3 y^5, Error = 0.000000
+DEAL::f(x,y) = x^4 y^0, Error = 0.000000
+DEAL::f(x,y) = x^4 y^1, Error = 0.000000
+DEAL::f(x,y) = x^4 y^2, Error = 0.000000
+DEAL::f(x,y) = x^4 y^3, Error = 0.000000
+DEAL::f(x,y) = x^4 y^4, Error = 0.000000
+DEAL::f(x,y) = x^4 y^5, Error = 0.000000
+DEAL::f(x,y) = x^5 y^0, Error = 0.000000
+DEAL::f(x,y) = x^5 y^1, Error = 0.000000
+DEAL::f(x,y) = x^5 y^2, Error = 0.000000
+DEAL::f(x,y) = x^5 y^3, Error = 0.000000
+DEAL::f(x,y) = x^5 y^4, Error = 0.000000
+DEAL::f(x,y) = x^5 y^5, Error = 0.000000