for (unsigned int i=0;i<p.size();++i)
for (unsigned int j=0;j<=i;++j)
deallog << 'P' << i << " * P" << j
- << " =" << scalar_product(p[i], p[j]) << std::endl;
+ << " = " << scalar_product(p[i], p[j]) << std::endl;
deallog << "LagrangeEquidistant" << std::endl;
for (unsigned int i=0;i<p.size();++i)
for (unsigned int j=0;j<p.size();++j)
deallog << 'P' << i << "(x" << j
- << ") =" << p[i].value((double) j/p.size())+1.0001 << std::endl;
+ << ") = " << p[i].value((double) j/p.size())+1.0001 << std::endl;
for (unsigned int i=0;i<p.size();++i)
{
for (unsigned int i=0;i<p.size();++i)
{
deallog << "N0(P" << i << ")"
- << " =" << p[i].value(0.) << std::endl;
+ << " = " << p[i].value(0.) << std::endl;
}
for (unsigned int i=0;i<p.size();++i)
{
deallog << "N1(P" << i << ")"
- << " =" << p[i].value(1.) << std::endl;
+ << " = " << p[i].value(1.) << std::endl;
}
std::vector<double> values (p.size());
{
p[i].value (.5, values);
deallog << "N" << j << "(P" << i << ")"
- << " =" << values[j] * factor << std::endl;
+ << " = " << values[j] * factor << std::endl;
}
}
}
derive
0.00000 x^0
DEAL::Legendre
-DEAL::P0 * P0 =1.00
-DEAL::P1 * P0 =0
-DEAL::P1 * P1 =1.000
-DEAL::P2 * P0 =0
-DEAL::P2 * P1 =0
-DEAL::P2 * P2 =1.00
-DEAL::P3 * P0 =0
-DEAL::P3 * P1 =0
-DEAL::P3 * P2 =0
-DEAL::P3 * P3 =1.000
-DEAL::P4 * P0 =0
-DEAL::P4 * P1 =0
-DEAL::P4 * P2 =0
-DEAL::P4 * P3 =0
-DEAL::P4 * P4 =1.00
-DEAL::P5 * P0 =0
-DEAL::P5 * P1 =0
-DEAL::P5 * P2 =0
-DEAL::P5 * P3 =0
-DEAL::P5 * P4 =0
-DEAL::P5 * P5 =1.00
-DEAL::P6 * P0 =0
-DEAL::P6 * P1 =0
-DEAL::P6 * P2 =0
-DEAL::P6 * P3 =0
-DEAL::P6 * P4 =0
-DEAL::P6 * P5 =0
-DEAL::P6 * P6 =1.000
-DEAL::P7 * P0 =0
-DEAL::P7 * P1 =0
-DEAL::P7 * P2 =0
-DEAL::P7 * P3 =0
-DEAL::P7 * P4 =0
-DEAL::P7 * P5 =0
-DEAL::P7 * P6 =0
-DEAL::P7 * P7 =1.00
-DEAL::P8 * P0 =0
-DEAL::P8 * P1 =0
-DEAL::P8 * P2 =0
-DEAL::P8 * P3 =0
-DEAL::P8 * P4 =0
-DEAL::P8 * P5 =0
-DEAL::P8 * P6 =0
-DEAL::P8 * P7 =0
-DEAL::P8 * P8 =1.000
-DEAL::P9 * P0 =0
-DEAL::P9 * P1 =0
-DEAL::P9 * P2 =0
-DEAL::P9 * P3 =0
-DEAL::P9 * P4 =0
-DEAL::P9 * P5 =0
-DEAL::P9 * P6 =0
-DEAL::P9 * P7 =0
-DEAL::P9 * P8 =0
-DEAL::P9 * P9 =1.00
-DEAL::P10 * P0 =1.46e-10
-DEAL::P10 * P1 =2.12e-10
-DEAL::P10 * P2 =1.96e-10
-DEAL::P10 * P3 =1.47e-10
-DEAL::P10 * P4 =1.26e-10
-DEAL::P10 * P5 =0
-DEAL::P10 * P6 =0
-DEAL::P10 * P7 =0
-DEAL::P10 * P8 =0
-DEAL::P10 * P9 =0
-DEAL::P10 * P10 =1.000
-DEAL::P11 * P0 =-6.91e-10
-DEAL::P11 * P1 =-5.76e-10
-DEAL::P11 * P2 =-4.43e-10
-DEAL::P11 * P3 =-5.78e-10
-DEAL::P11 * P4 =-1.55e-10
-DEAL::P11 * P5 =0
-DEAL::P11 * P6 =1.68e-10
-DEAL::P11 * P7 =-3.62e-10
-DEAL::P11 * P8 =-3.44e-10
-DEAL::P11 * P9 =0
-DEAL::P11 * P10 =0
-DEAL::P11 * P11 =1.00
+DEAL::P0 * P0 = 1.00
+DEAL::P1 * P0 = 0
+DEAL::P1 * P1 = 1.000
+DEAL::P2 * P0 = 0
+DEAL::P2 * P1 = 0
+DEAL::P2 * P2 = 1.00
+DEAL::P3 * P0 = 0
+DEAL::P3 * P1 = 0
+DEAL::P3 * P2 = 0
+DEAL::P3 * P3 = 1.000
+DEAL::P4 * P0 = 0
+DEAL::P4 * P1 = 0
+DEAL::P4 * P2 = 0
+DEAL::P4 * P3 = 0
+DEAL::P4 * P4 = 1.00
+DEAL::P5 * P0 = 0
+DEAL::P5 * P1 = 0
+DEAL::P5 * P2 = 0
+DEAL::P5 * P3 = 0
+DEAL::P5 * P4 = 0
+DEAL::P5 * P5 = 1.00
+DEAL::P6 * P0 = 0
+DEAL::P6 * P1 = 0
+DEAL::P6 * P2 = 0
+DEAL::P6 * P3 = 0
+DEAL::P6 * P4 = 0
+DEAL::P6 * P5 = 0
+DEAL::P6 * P6 = 1.000
+DEAL::P7 * P0 = 0
+DEAL::P7 * P1 = 0
+DEAL::P7 * P2 = 0
+DEAL::P7 * P3 = 0
+DEAL::P7 * P4 = 0
+DEAL::P7 * P5 = 0
+DEAL::P7 * P6 = 0
+DEAL::P7 * P7 = 1.00
+DEAL::P8 * P0 = 0
+DEAL::P8 * P1 = 0
+DEAL::P8 * P2 = 0
+DEAL::P8 * P3 = 0
+DEAL::P8 * P4 = 0
+DEAL::P8 * P5 = 0
+DEAL::P8 * P6 = 0
+DEAL::P8 * P7 = 0
+DEAL::P8 * P8 = 1.000
+DEAL::P9 * P0 = 0
+DEAL::P9 * P1 = 0
+DEAL::P9 * P2 = 0
+DEAL::P9 * P3 = 0
+DEAL::P9 * P4 = 0
+DEAL::P9 * P5 = 0
+DEAL::P9 * P6 = 0
+DEAL::P9 * P7 = 0
+DEAL::P9 * P8 = 0
+DEAL::P9 * P9 = 1.00
+DEAL::P10 * P0 = 1.46e-10
+DEAL::P10 * P1 = 2.12e-10
+DEAL::P10 * P2 = 1.96e-10
+DEAL::P10 * P3 = 1.47e-10
+DEAL::P10 * P4 = 1.26e-10
+DEAL::P10 * P5 = 0
+DEAL::P10 * P6 = 0
+DEAL::P10 * P7 = 0
+DEAL::P10 * P8 = 0
+DEAL::P10 * P9 = 0
+DEAL::P10 * P10 = 1.000
+DEAL::P11 * P0 = -6.91e-10
+DEAL::P11 * P1 = -5.76e-10
+DEAL::P11 * P2 = -4.43e-10
+DEAL::P11 * P3 = -5.78e-10
+DEAL::P11 * P4 = -1.55e-10
+DEAL::P11 * P5 = 0
+DEAL::P11 * P6 = 1.68e-10
+DEAL::P11 * P7 = -3.62e-10
+DEAL::P11 * P8 = -3.44e-10
+DEAL::P11 * P9 = 0
+DEAL::P11 * P10 = 0
+DEAL::P11 * P11 = 1.00
DEAL::LagrangeEquidistant
-DEAL::P0(x0) =2.00
-DEAL::P0(x1) =1.00
-DEAL::P0(x2) =1.00
-DEAL::P0(x3) =1.00
-DEAL::P0(x4) =1.00
-DEAL::P0(x5) =1.00
-DEAL::P1(x0) =1.00
-DEAL::P1(x1) =2.00
-DEAL::P1(x2) =1.00
-DEAL::P1(x3) =1.00
-DEAL::P1(x4) =1.00
-DEAL::P1(x5) =1.00
-DEAL::P2(x0) =1.00
-DEAL::P2(x1) =1.00
-DEAL::P2(x2) =2.00
-DEAL::P2(x3) =1.00
-DEAL::P2(x4) =1.00
-DEAL::P2(x5) =1.00
-DEAL::P3(x0) =1.00
-DEAL::P3(x1) =1.00
-DEAL::P3(x2) =1.00
-DEAL::P3(x3) =2.00
-DEAL::P3(x4) =1.00
-DEAL::P3(x5) =1.00
-DEAL::P4(x0) =1.00
-DEAL::P4(x1) =1.00
-DEAL::P4(x2) =1.00
-DEAL::P4(x3) =1.00
-DEAL::P4(x4) =2.00
-DEAL::P4(x5) =1.00
-DEAL::P5(x0) =1.00
-DEAL::P5(x1) =1.00
-DEAL::P5(x2) =1.00
-DEAL::P5(x3) =1.00
-DEAL::P5(x4) =1.00
-DEAL::P5(x5) =2.00
+DEAL::P0(x0) = 2.00
+DEAL::P0(x1) = 1.00
+DEAL::P0(x2) = 1.00
+DEAL::P0(x3) = 1.00
+DEAL::P0(x4) = 1.00
+DEAL::P0(x5) = 1.00
+DEAL::P1(x0) = 1.00
+DEAL::P1(x1) = 2.00
+DEAL::P1(x2) = 1.00
+DEAL::P1(x3) = 1.00
+DEAL::P1(x4) = 1.00
+DEAL::P1(x5) = 1.00
+DEAL::P2(x0) = 1.00
+DEAL::P2(x1) = 1.00
+DEAL::P2(x2) = 2.00
+DEAL::P2(x3) = 1.00
+DEAL::P2(x4) = 1.00
+DEAL::P2(x5) = 1.00
+DEAL::P3(x0) = 1.00
+DEAL::P3(x1) = 1.00
+DEAL::P3(x2) = 1.00
+DEAL::P3(x3) = 2.00
+DEAL::P3(x4) = 1.00
+DEAL::P3(x5) = 1.00
+DEAL::P4(x0) = 1.00
+DEAL::P4(x1) = 1.00
+DEAL::P4(x2) = 1.00
+DEAL::P4(x3) = 1.00
+DEAL::P4(x4) = 2.00
+DEAL::P4(x5) = 1.00
+DEAL::P5(x0) = 1.00
+DEAL::P5(x1) = 1.00
+DEAL::P5(x2) = 1.00
+DEAL::P5(x3) = 1.00
+DEAL::P5(x4) = 1.00
+DEAL::P5(x5) = 2.00
DEAL::Hierarchical
-DEAL::N0(P0) =1.00
-DEAL::N0(P1) =0
-DEAL::N0(P2) =0
-DEAL::N0(P3) =0
-DEAL::N0(P4) =0
-DEAL::N0(P5) =0
-DEAL::N0(P6) =0
-DEAL::N0(P7) =0
-DEAL::N0(P8) =0
-DEAL::N0(P9) =0
-DEAL::N0(P10) =0
-DEAL::N0(P11) =0
-DEAL::N1(P0) =0
-DEAL::N1(P1) =1.00
-DEAL::N1(P2) =0
-DEAL::N1(P3) =0
-DEAL::N1(P4) =0
-DEAL::N1(P5) =0
-DEAL::N1(P6) =0
-DEAL::N1(P7) =0
-DEAL::N1(P8) =0
-DEAL::N1(P9) =0
-DEAL::N1(P10) =0
-DEAL::N1(P11) =0
-DEAL::N2(P0) =0
-DEAL::N2(P1) =0
-DEAL::N2(P2) =1.00
-DEAL::N2(P3) =0
-DEAL::N2(P4) =0
-DEAL::N2(P5) =0
-DEAL::N2(P6) =0
-DEAL::N2(P7) =0
-DEAL::N2(P8) =0
-DEAL::N2(P9) =0
-DEAL::N2(P10) =0
-DEAL::N2(P11) =0
-DEAL::N3(P0) =0
-DEAL::N3(P1) =0
-DEAL::N3(P2) =0
-DEAL::N3(P3) =1.00
-DEAL::N3(P4) =0
-DEAL::N3(P5) =0
-DEAL::N3(P6) =0
-DEAL::N3(P7) =0
-DEAL::N3(P8) =0
-DEAL::N3(P9) =0
-DEAL::N3(P10) =0
-DEAL::N3(P11) =0
-DEAL::N4(P0) =0
-DEAL::N4(P1) =0
-DEAL::N4(P2) =0
-DEAL::N4(P3) =0
-DEAL::N4(P4) =1.00
-DEAL::N4(P5) =0
-DEAL::N4(P6) =0
-DEAL::N4(P7) =0
-DEAL::N4(P8) =0
-DEAL::N4(P9) =0
-DEAL::N4(P10) =0
-DEAL::N4(P11) =0
-DEAL::N5(P0) =0
-DEAL::N5(P1) =0
-DEAL::N5(P2) =0
-DEAL::N5(P3) =0
-DEAL::N5(P4) =0
-DEAL::N5(P5) =1.00
-DEAL::N5(P6) =0
-DEAL::N5(P7) =0
-DEAL::N5(P8) =0
-DEAL::N5(P9) =0
-DEAL::N5(P10) =0
-DEAL::N5(P11) =0
-DEAL::N6(P0) =0
-DEAL::N6(P1) =0
-DEAL::N6(P2) =0
-DEAL::N6(P3) =0
-DEAL::N6(P4) =0
-DEAL::N6(P5) =0
-DEAL::N6(P6) =1.00
-DEAL::N6(P7) =0
-DEAL::N6(P8) =0
-DEAL::N6(P9) =0
-DEAL::N6(P10) =0
-DEAL::N6(P11) =0
-DEAL::N7(P0) =0
-DEAL::N7(P1) =0
-DEAL::N7(P2) =0
-DEAL::N7(P3) =0
-DEAL::N7(P4) =0
-DEAL::N7(P5) =0
-DEAL::N7(P6) =0
-DEAL::N7(P7) =1.000
-DEAL::N7(P8) =0
-DEAL::N7(P9) =0
-DEAL::N7(P10) =0
-DEAL::N7(P11) =0
-DEAL::N8(P0) =0
-DEAL::N8(P1) =0
-DEAL::N8(P2) =0
-DEAL::N8(P3) =0
-DEAL::N8(P4) =0
-DEAL::N8(P5) =0
-DEAL::N8(P6) =0
-DEAL::N8(P7) =0
-DEAL::N8(P8) =1.000
-DEAL::N8(P9) =0
-DEAL::N8(P10) =0
-DEAL::N8(P11) =0
-DEAL::N9(P0) =0
-DEAL::N9(P1) =0
-DEAL::N9(P2) =0
-DEAL::N9(P3) =0
-DEAL::N9(P4) =0
-DEAL::N9(P5) =0
-DEAL::N9(P6) =0
-DEAL::N9(P7) =0
-DEAL::N9(P8) =0
-DEAL::N9(P9) =1.000
-DEAL::N9(P10) =0
-DEAL::N9(P11) =0
-DEAL::N10(P0) =0
-DEAL::N10(P1) =0
-DEAL::N10(P2) =0
-DEAL::N10(P3) =0
-DEAL::N10(P4) =0
-DEAL::N10(P5) =0
-DEAL::N10(P6) =0
-DEAL::N10(P7) =0
-DEAL::N10(P8) =0
-DEAL::N10(P9) =0
-DEAL::N10(P10) =1.00
-DEAL::N10(P11) =0
-DEAL::N11(P0) =0
-DEAL::N11(P1) =0
-DEAL::N11(P2) =0
-DEAL::N11(P3) =0
-DEAL::N11(P4) =0
-DEAL::N11(P5) =0
-DEAL::N11(P6) =0
-DEAL::N11(P7) =0
-DEAL::N11(P8) =0
-DEAL::N11(P9) =0
-DEAL::N11(P10) =0
-DEAL::N11(P11) =1.000
+DEAL::N0(P0) = 1.00
+DEAL::N0(P1) = 0
+DEAL::N0(P2) = 0
+DEAL::N0(P3) = 0
+DEAL::N0(P4) = 0
+DEAL::N0(P5) = 0
+DEAL::N0(P6) = 0
+DEAL::N0(P7) = 0
+DEAL::N0(P8) = 0
+DEAL::N0(P9) = 0
+DEAL::N0(P10) = 0
+DEAL::N0(P11) = 0
+DEAL::N1(P0) = 0
+DEAL::N1(P1) = 1.00
+DEAL::N1(P2) = 0
+DEAL::N1(P3) = 0
+DEAL::N1(P4) = 0
+DEAL::N1(P5) = 0
+DEAL::N1(P6) = 0
+DEAL::N1(P7) = 0
+DEAL::N1(P8) = 0
+DEAL::N1(P9) = 0
+DEAL::N1(P10) = 0
+DEAL::N1(P11) = 0
+DEAL::N2(P0) = 0
+DEAL::N2(P1) = 0
+DEAL::N2(P2) = 1.00
+DEAL::N2(P3) = 0
+DEAL::N2(P4) = 0
+DEAL::N2(P5) = 0
+DEAL::N2(P6) = 0
+DEAL::N2(P7) = 0
+DEAL::N2(P8) = 0
+DEAL::N2(P9) = 0
+DEAL::N2(P10) = 0
+DEAL::N2(P11) = 0
+DEAL::N3(P0) = 0
+DEAL::N3(P1) = 0
+DEAL::N3(P2) = 0
+DEAL::N3(P3) = 1.00
+DEAL::N3(P4) = 0
+DEAL::N3(P5) = 0
+DEAL::N3(P6) = 0
+DEAL::N3(P7) = 0
+DEAL::N3(P8) = 0
+DEAL::N3(P9) = 0
+DEAL::N3(P10) = 0
+DEAL::N3(P11) = 0
+DEAL::N4(P0) = 0
+DEAL::N4(P1) = 0
+DEAL::N4(P2) = 0
+DEAL::N4(P3) = 0
+DEAL::N4(P4) = 1.00
+DEAL::N4(P5) = 0
+DEAL::N4(P6) = 0
+DEAL::N4(P7) = 0
+DEAL::N4(P8) = 0
+DEAL::N4(P9) = 0
+DEAL::N4(P10) = 0
+DEAL::N4(P11) = 0
+DEAL::N5(P0) = 0
+DEAL::N5(P1) = 0
+DEAL::N5(P2) = 0
+DEAL::N5(P3) = 0
+DEAL::N5(P4) = 0
+DEAL::N5(P5) = 1.00
+DEAL::N5(P6) = 0
+DEAL::N5(P7) = 0
+DEAL::N5(P8) = 0
+DEAL::N5(P9) = 0
+DEAL::N5(P10) = 0
+DEAL::N5(P11) = 0
+DEAL::N6(P0) = 0
+DEAL::N6(P1) = 0
+DEAL::N6(P2) = 0
+DEAL::N6(P3) = 0
+DEAL::N6(P4) = 0
+DEAL::N6(P5) = 0
+DEAL::N6(P6) = 1.00
+DEAL::N6(P7) = 0
+DEAL::N6(P8) = 0
+DEAL::N6(P9) = 0
+DEAL::N6(P10) = 0
+DEAL::N6(P11) = 0
+DEAL::N7(P0) = 0
+DEAL::N7(P1) = 0
+DEAL::N7(P2) = 0
+DEAL::N7(P3) = 0
+DEAL::N7(P4) = 0
+DEAL::N7(P5) = 0
+DEAL::N7(P6) = 0
+DEAL::N7(P7) = 1.000
+DEAL::N7(P8) = 0
+DEAL::N7(P9) = 0
+DEAL::N7(P10) = 0
+DEAL::N7(P11) = 0
+DEAL::N8(P0) = 0
+DEAL::N8(P1) = 0
+DEAL::N8(P2) = 0
+DEAL::N8(P3) = 0
+DEAL::N8(P4) = 0
+DEAL::N8(P5) = 0
+DEAL::N8(P6) = 0
+DEAL::N8(P7) = 0
+DEAL::N8(P8) = 1.000
+DEAL::N8(P9) = 0
+DEAL::N8(P10) = 0
+DEAL::N8(P11) = 0
+DEAL::N9(P0) = 0
+DEAL::N9(P1) = 0
+DEAL::N9(P2) = 0
+DEAL::N9(P3) = 0
+DEAL::N9(P4) = 0
+DEAL::N9(P5) = 0
+DEAL::N9(P6) = 0
+DEAL::N9(P7) = 0
+DEAL::N9(P8) = 0
+DEAL::N9(P9) = 1.000
+DEAL::N9(P10) = 0
+DEAL::N9(P11) = 0
+DEAL::N10(P0) = 0
+DEAL::N10(P1) = 0
+DEAL::N10(P2) = 0
+DEAL::N10(P3) = 0
+DEAL::N10(P4) = 0
+DEAL::N10(P5) = 0
+DEAL::N10(P6) = 0
+DEAL::N10(P7) = 0
+DEAL::N10(P8) = 0
+DEAL::N10(P9) = 0
+DEAL::N10(P10) = 1.00
+DEAL::N10(P11) = 0
+DEAL::N11(P0) = 0
+DEAL::N11(P1) = 0
+DEAL::N11(P2) = 0
+DEAL::N11(P3) = 0
+DEAL::N11(P4) = 0
+DEAL::N11(P5) = 0
+DEAL::N11(P6) = 0
+DEAL::N11(P7) = 0
+DEAL::N11(P8) = 0
+DEAL::N11(P9) = 0
+DEAL::N11(P10) = 0
+DEAL::N11(P11) = 1.000