Run grand replace for 2015.

[lilypond.git] / flower / include / polynomial.hh

diff --git a/flower/include/polynomial.hh b/flower/include/polynomial.hh
index df79274bdc96d45e9408953b4655327cde2d77f3..f80b53c2eec1d62355bbb984bf617ed3859d3e0f 100644 (file)
--- a/flower/include/polynomial.hh
+++ b/flower/include/polynomial.hh
@@ -1,84 +1,99 @@
+/*
+  This file is part of LilyPond, the GNU music typesetter.
+
+  Copyright (C) 1993--2015 Han-Wen Nienhuys <hanwen@xs4all.nl>
+
+  LilyPond is free software: you can redistribute it and/or modify
+  it under the terms of the GNU General Public License as published by
+  the Free Software Foundation, either version 3 of the License, or
+  (at your option) any later version.
+
+  LilyPond is distributed in the hope that it will be useful,
+  but WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+  GNU General Public License for more details.
+
+  You should have received a copy of the GNU General Public License
+  along with LilyPond.  If not, see <http://www.gnu.org/licenses/>.
+*/
 
 /*
- * poly.h -- routines for manipulation of polynomials in one var
- *
- * (c) 1993, 1996,1999 Han-Wen Nienhuys
+ * polynomial.hh -- routines for manipulation of polynomials in one var
  */
 
 #ifndef POLY_H
 #define POLY_H
 
-#include "array.hh"
+#include <sys/types.h>
+#include "std-vector.hh"
 #include "arithmetic-operator.hh"
 #include "real.hh"
 
-/// structure for a polynomial in one var. 
+/// structure for a polynomial in one var.
 struct Polynomial
 {
   /// degree of polynomial
-  int degree ()const;
+  ssize_t degree ()const;
 
-  /// coefficients 
-  Array<Real>     coefs_;
+  /// coefficients
+  vector<Real> coefs_;
 
   // leading coef
-  Real &lc();
+  Real &lc ();
 
   // leading coef
-  Real lc() const;
-  void print() const ; 
-  Real eval(Real) const ;
-  void print_sols(Array<Real>) const ;
-  void check_sols(Array<Real>) const ;
-  void check_sol(Real x) const;
-  static Polynomial multiply(const Polynomial & p1, const Polynomial & p2);
-  static Polynomial power(int exponent, const Polynomial & src);
+  Real lc () const;
+  void print () const;
+  Real eval (Real) const;
+  Real minmax (Real, Real, bool) const;
+  void print_sols (vector<Real>) const;
+  void check_sols (vector<Real>) const;
+  void check_sol (Real x) const;
+  static Polynomial multiply (const Polynomial &p1, const Polynomial &p2);
+  static Polynomial power (int exponent, const Polynomial &src);
 
   /// chop low coefficients
-  void clean();
+  void clean ();
 
   /// eliminate #x#  close to  zero
-  void real_clean();
-  static Polynomial add(const Polynomial & p1, const Polynomial & p2);
-  void scalarmultiply(Real fact);
-  void operator *= (Real f) { scalarmultiply(f); }
-  void operator /= (Real f) { scalarmultiply(1/f); }
+  void real_clean ();
+  void scalarmultiply (Real fact);
+  void operator *= (Real f) { scalarmultiply (f); }
+  void operator /= (Real f) { scalarmultiply (1 / f); }
   void operator += (Polynomial const &p2);
   void operator *= (Polynomial const &p2);
   void operator -= (Polynomial const &p2);
-  Polynomial (Real a, Real b =0.0);
-  Polynomial (){}
-  static Polynomial subtract(const Polynomial & p1, const Polynomial & p2);
-  void set_negate(const Polynomial & src);
-    
+  Polynomial (Real a, Real b = 0.0);
+  Polynomial () {}
+  void set_negate (const Polynomial &src);
+
   /// take the derivative
-  void differentiate();
-  int set_mod(const Polynomial &u, const Polynomial &v);
+  void differentiate ();
+  int set_mod (const Polynomial &u, const Polynomial &v);
 
-  void debug_clean();
+  void debug_clean ();
 
-  Array<Real> solve_quadric()const;
-  Array<Real> solve_cubic()const;
-  Array<Real> solve_linear()const;
+  vector<Real> solve_quadric ()const;
+  vector<Real> solve_cubic ()const;
+  vector<Real> solve_linear ()const;
 
-  Array<Real> solve () const;
+  vector<Real> solve () const;
 };
 
-
-IMPLEMENT_ARITHMETIC_OPERATOR(Polynomial, - );
-IMPLEMENT_ARITHMETIC_OPERATOR(Polynomial, + );
-IMPLEMENT_ARITHMETIC_OPERATOR(Polynomial, * );
+IMPLEMENT_ARITHMETIC_OPERATOR (Polynomial, -);
+IMPLEMENT_ARITHMETIC_OPERATOR (Polynomial, +);
+IMPLEMENT_ARITHMETIC_OPERATOR (Polynomial, *);
 
 inline Polynomial
 operator * (Polynomial p, Real a)
 {
-  p *=a;
+  p *= a;
   return p;
 }
 inline Polynomial
-operator * (Real a,Polynomial p)
+operator * (Real a, Polynomial p)
 {
-  p *=a;
+  p *= a;
   return p;
 }
 #endif