2 This file is part of LilyPond, the GNU music typesetter.
4 Copyright (C) 1993--2014 Han-Wen Nienhuys <hanwen@xs4all.nl>
6 LilyPond is free software: you can redistribute it and/or modify
7 it under the terms of the GNU General Public License as published by
8 the Free Software Foundation, either version 3 of the License, or
9 (at your option) any later version.
11 LilyPond is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 GNU General Public License for more details.
16 You should have received a copy of the GNU General Public License
17 along with LilyPond. If not, see <http://www.gnu.org/licenses/>.
21 * polynomial.hh -- routines for manipulation of polynomials in one var
27 #include <sys/types.h>
28 #include "std-vector.hh"
29 #include "arithmetic-operator.hh"
32 /// structure for a polynomial in one var.
35 /// degree of polynomial
36 ssize_t degree ()const;
47 Real eval (Real) const;
48 Real minmax (Real, Real, bool) const;
49 void print_sols (vector<Real>) const;
50 void check_sols (vector<Real>) const;
51 void check_sol (Real x) const;
52 static Polynomial multiply (const Polynomial &p1, const Polynomial &p2);
53 static Polynomial power (int exponent, const Polynomial &src);
55 /// chop low coefficients
58 /// eliminate #x# close to zero
60 void scalarmultiply (Real fact);
61 void operator *= (Real f) { scalarmultiply (f); }
62 void operator /= (Real f) { scalarmultiply (1 / f); }
63 void operator += (Polynomial const &p2);
64 void operator *= (Polynomial const &p2);
65 void operator -= (Polynomial const &p2);
66 Polynomial (Real a, Real b = 0.0);
68 void set_negate (const Polynomial &src);
70 /// take the derivative
71 void differentiate ();
72 int set_mod (const Polynomial &u, const Polynomial &v);
76 vector<Real> solve_quadric ()const;
77 vector<Real> solve_cubic ()const;
78 vector<Real> solve_linear ()const;
80 vector<Real> solve () const;
83 IMPLEMENT_ARITHMETIC_OPERATOR (Polynomial, -);
84 IMPLEMENT_ARITHMETIC_OPERATOR (Polynomial, +);
85 IMPLEMENT_ARITHMETIC_OPERATOR (Polynomial, *);
88 operator * (Polynomial p, Real a)
94 operator * (Real a, Polynomial p)