+/*
+ This file is part of LilyPond, the GNU music typesetter.
+
+ Copyright (C) 1993--2011 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 "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;
- /// coefficients
- Array<Real> coefs_;
+ /// coefficients
+ vector<Real> coefs_;
// leading coef
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 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);
+ static Polynomial multiply (const Polynomial &p1, const Polynomial &p2);
+ static Polynomial power (int exponent, const Polynomial &src);
/// chop low coefficients
void clean ();
void real_clean ();
void scalarmultiply (Real fact);
void operator *= (Real f) { scalarmultiply (f); }
- void operator /= (Real f) { scalarmultiply (1/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 (){}
- 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 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