From 68e3d0f7e4916dc3b8b5cf12901c438c779d63a7 Mon Sep 17 00:00:00 2001 From: Guido Kanschat 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 ++++++++++++++++++++++++++++++++++++++ 1 file changed, 45 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; -- 2.39.5