add a constructor for a zero polynomial of degree n

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

diff --git a/deal.II/base/include/base/polynomial.h b/deal.II/base/include/base/polynomial.h
index 0dfaf086c6b9c5c841d2ac68b9b266b22431a8a1..330a6c411fcc4010a4f00ea002d9172392e43814 100644 (file)
--- a/deal.II/base/include/base/polynomial.h
+++ b/deal.II/base/include/base/polynomial.h
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006 by the deal.II authors
+//    Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2009 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -43,7 +43,7 @@ namespace Polynomials
  * happens through the Horner scheme which provides both numerical
  * stability and a minimal number of numerical operations.
  *
- * @author Ralf Hartmann, Guido Kanschat, 2000, 2006
+ * @author Ralf Hartmann, Guido Kanschat, 2000, 2006, 2009
  */
   template <typename number>
   class Polynomial : public Subscriptor
@@ -68,6 +68,12 @@ namespace Polynomials
                                         */
       Polynomial (const std::vector<number> &coefficients);
 
+                                      /**
+                                       * Constructor creating a zero
+                                       * polynomial of degree @p n.
+                                       */
+      Polynomial (const unsigned int n);
+      
                                       /**
                                        * Default constructor creating
                                        * an illegal object.

diff --git a/deal.II/base/source/polynomial.cc b/deal.II/base/source/polynomial.cc
index a460dc75be1b13a2488bce062c81a76eec66ca68..df38a18697f194c651238dcf6e6db64d4b83f7c2 100644 (file)
--- a/deal.II/base/source/polynomial.cc
+++ b/deal.II/base/source/polynomial.cc
@@ -2,7 +2,7 @@
 //      $Id$   
 //    Version: $Name$
 //
-//    Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006 by the deal.II authors
+//    Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2009 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -52,6 +52,14 @@ namespace Polynomials
 
 
 
+  template <typename number>
+  Polynomial<number>::Polynomial (const unsigned int n)
+                  :
+                  coefficients(n+1, 0.)
+  {}
+
+
+
   template <typename number>
   void
   Polynomial<number>::value (const number         x,