From 32e7b17e69e5ca8f67036b6cc2efa45e766a854d Mon Sep 17 00:00:00 2001 From: guido Date: Tue, 16 Dec 2003 11:29:22 +0000 Subject: [PATCH] check new arithmetic functions git-svn-id: https://svn.dealii.org/trunk@8268 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/base/polynomial1d.cc | 45 +++++++++++++++ .../base/polynomial1d.output | 57 +++++++++++++++++++ 2 files changed, 102 insertions(+) diff --git a/tests/base/polynomial1d.cc b/tests/base/polynomial1d.cc index 3e48238535..303775f7df 100644 --- a/tests/base/polynomial1d.cc +++ b/tests/base/polynomial1d.cc @@ -44,13 +44,58 @@ double scalar_product (const Polynomial& p1, return sum; } + +void polynomial_arithmetic () +{ + std::vector c1(4); + c1[0] = 2; + c1[1] = 3; + c1[2] = 5; + c1[3] = 4; + Polynomial p1(c1); + + std::vector c2(3); + c2[0] = .4; + c2[1] = .7; + c2[2] = -1.3; + Polynomial p2(c2); + Monomial p3(5); + std::cerr << "P1" << std::endl; + p1.print(std::cerr); + std::cerr << "P2" << std::endl; + p2.print(std::cerr); + p1 += p2; + std::cerr << "P1+P2" << std::endl; + p1.print(std::cerr); + p2 += p1; + std::cerr << "P1+2P2" << std::endl; + p2.print(std::cerr); + std::cerr << "P1+P2+x^5" << std::endl; + p1 += p3; + p1.print(std::cerr); + p1 *= 2.; + std::cerr << "*2" << std::endl; + p1.print(std::cerr); + + for (unsigned int i=0;i<7;++i) + { + std::cerr << "derive" << std::endl; + p1 = p1.derivative(); + p1.print(std::cerr); + } +} + + int main () { std::ofstream logfile("polynomial1d.output"); logfile.precision(3); deallog.attach(logfile); deallog.depth_console(0); + deallog.log_cerr(); + polynomial_arithmetic(); + std::vector > p; std::vector > q; diff --git a/tests/results/i686-pc-linux-gnu+gcc3.2/base/polynomial1d.output b/tests/results/i686-pc-linux-gnu+gcc3.2/base/polynomial1d.output index 9920cfa1f6..4c20674288 100644 --- a/tests/results/i686-pc-linux-gnu+gcc3.2/base/polynomial1d.output +++ b/tests/results/i686-pc-linux-gnu+gcc3.2/base/polynomial1d.output @@ -1,4 +1,61 @@ +P1 +4.00000 x^3 +5.00000 x^2 +3.00000 x^1 +2.00000 x^0 +P2 +-1.30000 x^2 +0.700000 x^1 +0.400000 x^0 +P1+P2 +4.00000 x^3 +3.70000 x^2 +3.70000 x^1 +2.40000 x^0 +P1+2P2 +4.00000 x^3 +2.40000 x^2 +4.40000 x^1 +2.80000 x^0 +P1+P2+x^5 +1.00000 x^5 +0.00000 x^4 +4.00000 x^3 +3.70000 x^2 +3.70000 x^1 +2.40000 x^0 +*2 +2.00000 x^5 +0.00000 x^4 +8.00000 x^3 +7.40000 x^2 +7.40000 x^1 +4.80000 x^0 +derive +10.0000 x^4 +0.00000 x^3 +24.0000 x^2 +14.8000 x^1 +7.40000 x^0 +derive +40.0000 x^3 +0.00000 x^2 +48.0000 x^1 +14.8000 x^0 +derive +120.000 x^2 +0.00000 x^1 +48.0000 x^0 +derive +240.000 x^1 +0.00000 x^0 +derive +240.000 x^0 +derive +0.00000 x^0 +derive +0.00000 x^0 DEAL::Legendre DEAL::P0 * P0 =1.00 DEAL::P1 * P0 =0.00 -- 2.39.5