]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
added Hierarchical test
authorbrian <brian@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 12 Feb 2003 17:00:44 +0000 (17:00 +0000)
committerbrian <brian@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 12 Feb 2003 17:00:44 +0000 (17:00 +0000)
git-svn-id: https://svn.dealii.org/trunk@7081 0785d39b-7218-0410-832d-ea1e28bc413d

tests/base/polynomial1d.cc
tests/base/polynomial1d.checked

index 975d30db7daabfd528fd2e3ce1c343973b3f3894..e7a94529411a24f43597cf123dfb584240456a5c 100644 (file)
@@ -15,6 +15,7 @@
 // 1d-polynomials, to make sure that all changes we make to these classes
 // do not change the results of these classes.
 
+// added Hierarchical test on 2/11/03 - Brian Carnes
 
 #include <iostream>
 #include <fstream>
@@ -114,5 +115,41 @@ int main ()
                    << std::endl;
        }
     }
+
+  deallog << "Hierarchical" << std::endl;
+  
+  p.clear ();
+  for (unsigned int i=0;i<12;++i)
+    {
+      p.push_back (Hierarchical<double>(i));
+    }
+
+  for (unsigned int i=0;i<p.size();++i)
+    {
+      deallog << "N0(P" << i << ")"
+             << " =" << 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;
+    }
+
+  std::vector<double> values (p.size());
+
+  for (unsigned int j=2; j<p.size(); ++j)
+    {
+      double factor = 1.;
+      for (unsigned int k=0; k<j;++k)
+       factor /= 2. * (j-k);
+
+      for (unsigned int i=0;i<p.size();++i)
+       {
+         p[i].value (.5, values);
+         deallog << "N" << j << "(P" << i << ")"
+                 << " =" << values[j] * factor << std::endl;
+       }
+    }
 }
 
index 2f5d26cd6efd1a1a83b829a1eec7c00487b6a27f..b7295538b0a9f3d252c6a18b623e1fcdc6fd8ada 100644 (file)
@@ -2,14 +2,14 @@
 DEAL::Legendre
 DEAL::P0 * P0 =1.00
 DEAL::P1 * P0 =0.00
-DEAL::P1 * P1 =1.00
+DEAL::P1 * P1 =1.000
 DEAL::P2 * P0 =0.00
 DEAL::P2 * P1 =0.00
 DEAL::P2 * P2 =1.00
 DEAL::P3 * P0 =0.00
 DEAL::P3 * P1 =0.00
 DEAL::P3 * P2 =0.00
-DEAL::P3 * P3 =1.00
+DEAL::P3 * P3 =1.000
 DEAL::P4 * P0 =0.00
 DEAL::P4 * P1 =0.00
 DEAL::P4 * P2 =0.00
@@ -27,7 +27,7 @@ DEAL::P6 * P2 =0.00
 DEAL::P6 * P3 =0.00
 DEAL::P6 * P4 =0.00
 DEAL::P6 * P5 =0.00
-DEAL::P6 * P6 =1.00
+DEAL::P6 * P6 =1.000
 DEAL::P7 * P0 =0.00
 DEAL::P7 * P1 =0.00
 DEAL::P7 * P2 =0.00
@@ -44,7 +44,7 @@ DEAL::P8 * P4 =0.00
 DEAL::P8 * P5 =0.00
 DEAL::P8 * P6 =0.00
 DEAL::P8 * P7 =0.00
-DEAL::P8 * P8 =1.00
+DEAL::P8 * P8 =1.000
 DEAL::P9 * P0 =0.00
 DEAL::P9 * P1 =0.00
 DEAL::P9 * P2 =0.00
@@ -65,7 +65,7 @@ DEAL::P10 * P6 =0.00
 DEAL::P10 * P7 =0.00
 DEAL::P10 * P8 =0.00
 DEAL::P10 * P9 =0.00
-DEAL::P10 * P10 =1.00
+DEAL::P10 * P10 =1.000
 DEAL::P11 * P0 =0.00
 DEAL::P11 * P1 =0.00
 DEAL::P11 * P2 =0.00
@@ -115,3 +115,148 @@ 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.00
+DEAL::N0(P2) =0.00
+DEAL::N0(P3) =0.00
+DEAL::N0(P4) =0.00
+DEAL::N0(P5) =0.00
+DEAL::N0(P6) =0.00
+DEAL::N0(P7) =0.00
+DEAL::N0(P8) =0.00
+DEAL::N0(P9) =0.00
+DEAL::N0(P10) =0.00
+DEAL::N0(P11) =0.00
+DEAL::N1(P0) =0.00
+DEAL::N1(P1) =1.00
+DEAL::N1(P2) =0.00
+DEAL::N1(P3) =0.00
+DEAL::N1(P4) =0.00
+DEAL::N1(P5) =0.00
+DEAL::N1(P6) =0.00
+DEAL::N1(P7) =0.00
+DEAL::N1(P8) =0.00
+DEAL::N1(P9) =0.00
+DEAL::N1(P10) =0.00
+DEAL::N1(P11) =0.00
+DEAL::N2(P0) =0.00
+DEAL::N2(P1) =0.00
+DEAL::N2(P2) =1.00
+DEAL::N2(P3) =0.00
+DEAL::N2(P4) =0.00
+DEAL::N2(P5) =0.00
+DEAL::N2(P6) =0.00
+DEAL::N2(P7) =0.00
+DEAL::N2(P8) =0.00
+DEAL::N2(P9) =0.00
+DEAL::N2(P10) =0.00
+DEAL::N2(P11) =0.00
+DEAL::N3(P0) =0.00
+DEAL::N3(P1) =0.00
+DEAL::N3(P2) =0.00
+DEAL::N3(P3) =1.00
+DEAL::N3(P4) =0.00
+DEAL::N3(P5) =0.00
+DEAL::N3(P6) =0.00
+DEAL::N3(P7) =0.00
+DEAL::N3(P8) =0.00
+DEAL::N3(P9) =0.00
+DEAL::N3(P10) =0.00
+DEAL::N3(P11) =0.00
+DEAL::N4(P0) =0.00
+DEAL::N4(P1) =0.00
+DEAL::N4(P2) =0.00
+DEAL::N4(P3) =0.00
+DEAL::N4(P4) =1.00
+DEAL::N4(P5) =0.00
+DEAL::N4(P6) =0.00
+DEAL::N4(P7) =0.00
+DEAL::N4(P8) =0.00
+DEAL::N4(P9) =0.00
+DEAL::N4(P10) =0.00
+DEAL::N4(P11) =0.00
+DEAL::N5(P0) =0.00
+DEAL::N5(P1) =0.00
+DEAL::N5(P2) =0.00
+DEAL::N5(P3) =0.00
+DEAL::N5(P4) =0.00
+DEAL::N5(P5) =1.00
+DEAL::N5(P6) =0.00
+DEAL::N5(P7) =0.00
+DEAL::N5(P8) =0.00
+DEAL::N5(P9) =0.00
+DEAL::N5(P10) =0.00
+DEAL::N5(P11) =0.00
+DEAL::N6(P0) =0.00
+DEAL::N6(P1) =0.00
+DEAL::N6(P2) =0.00
+DEAL::N6(P3) =0.00
+DEAL::N6(P4) =0.00
+DEAL::N6(P5) =0.00
+DEAL::N6(P6) =1.00
+DEAL::N6(P7) =0.00
+DEAL::N6(P8) =0.00
+DEAL::N6(P9) =0.00
+DEAL::N6(P10) =0.00
+DEAL::N6(P11) =0.00
+DEAL::N7(P0) =0.00
+DEAL::N7(P1) =0.00
+DEAL::N7(P2) =0.00
+DEAL::N7(P3) =0.00
+DEAL::N7(P4) =0.00
+DEAL::N7(P5) =0.00
+DEAL::N7(P6) =0.00
+DEAL::N7(P7) =1.000
+DEAL::N7(P8) =0.00
+DEAL::N7(P9) =0.00
+DEAL::N7(P10) =0.00
+DEAL::N7(P11) =0.00
+DEAL::N8(P0) =0.00
+DEAL::N8(P1) =0.00
+DEAL::N8(P2) =0.00
+DEAL::N8(P3) =0.00
+DEAL::N8(P4) =0.00
+DEAL::N8(P5) =0.00
+DEAL::N8(P6) =0.00
+DEAL::N8(P7) =0.00
+DEAL::N8(P8) =1.000
+DEAL::N8(P9) =0.00
+DEAL::N8(P10) =0.00
+DEAL::N8(P11) =0.00
+DEAL::N9(P0) =0.00
+DEAL::N9(P1) =0.00
+DEAL::N9(P2) =0.00
+DEAL::N9(P3) =0.00
+DEAL::N9(P4) =0.00
+DEAL::N9(P5) =0.00
+DEAL::N9(P6) =0.00
+DEAL::N9(P7) =0.00
+DEAL::N9(P8) =0.00
+DEAL::N9(P9) =1.000
+DEAL::N9(P10) =0.00
+DEAL::N9(P11) =0.00
+DEAL::N10(P0) =0.00
+DEAL::N10(P1) =0.00
+DEAL::N10(P2) =0.00
+DEAL::N10(P3) =0.00
+DEAL::N10(P4) =0.00
+DEAL::N10(P5) =0.00
+DEAL::N10(P6) =0.00
+DEAL::N10(P7) =0.00
+DEAL::N10(P8) =0.00
+DEAL::N10(P9) =0.00
+DEAL::N10(P10) =1.00
+DEAL::N10(P11) =0.00
+DEAL::N11(P0) =0.00
+DEAL::N11(P1) =0.00
+DEAL::N11(P2) =0.00
+DEAL::N11(P3) =0.00
+DEAL::N11(P4) =0.00
+DEAL::N11(P5) =0.00
+DEAL::N11(P6) =0.00
+DEAL::N11(P7) =0.00
+DEAL::N11(P8) =0.00
+DEAL::N11(P9) =0.00
+DEAL::N11(P10) =0.00
+DEAL::N11(P11) =1.000

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.