/usr/include/polymake/tropical/polynomial_tools.h is in libpolymake-dev-common 3.2r2-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 | /*
This program 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 2
of the License, or (at your option) any later version.
This program 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 this program; if not, write to the Free Software
Foundation, Inc., 51 Franklin Street, Fifth Floor,
Boston, MA 02110-1301, USA.
---
Copyright (C) 2011 - 2015, Simon Hampe <simon.hampe@googlemail.com>
FIXME Most (or all) of these should at some time be implemented in
Polynomial.
*/
#ifndef POLYMAKE_ATINT_POLYNOMIAL_TOOLS_H
#define POLYMAKE_ATINT_POLYNOMIAL_TOOLS_H
#include "polymake/client.h"
#include "polymake/Matrix.h"
#include "polymake/Vector.h"
#include "polymake/Rational.h"
#include "polymake/Polynomial.h"
#include "polymake/TropicalNumber.h"
namespace polymake { namespace tropical {
template <typename Addition>
Rational evaluate_polynomial(const Polynomial<TropicalNumber<Addition>> &p, const Vector<Rational> &v){
Matrix<Rational> monoms(p.monomials_as_matrix());
Vector<TropicalNumber<Addition>> coefs(p.coefficients_as_vector());
TropicalNumber<Addition> result = TropicalNumber<Addition>::zero();
for(int m = 0; m < monoms.rows(); m++) {
result += (coefs[m] * TropicalNumber<Addition>(monoms.row(m)*v));
}
return Rational(result);
}//END evaluate_polynomial
template <typename Coefficient>
Vector<int> degree_vector(const Polynomial<Coefficient> &p){
return accumulate( cols(p.monomials_as_matrix()), operations::add());
}
template <typename Coefficient>
int polynomial_degree(const Polynomial<Coefficient> &p) {
if(p.monomials_as_matrix().rows() == 0) return -1;
return accumulate( degree_vector(p), operations::max());
}
template <typename Coefficient>
bool is_homogeneous(const Polynomial<Coefficient> &p) {
if(p.monomials_as_matrix().rows() == 0) return true;
Vector<int> dv = degree_vector(p);
return dv == dv[0] * ones_vector<int>(dv.dim());
}
} }
#endif
|