X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Fnormal.h;h=5a735e346a8f485bfeaf1811ba329951b90c39c9;hp=304d51c2b50bbb492e45dc03c4bc4fe313e411f4;hb=c277943ace48d2f22f7016ba9439ff7cb949a231;hpb=dbd9c306a74f1cb258c0d15a346b973b39deaad2 diff --git a/ginac/normal.h b/ginac/normal.h index 304d51c2..5a735e34 100644 --- a/ginac/normal.h +++ b/ginac/normal.h @@ -6,7 +6,7 @@ * computation, square-free factorization and rational function normalization. */ /* - * GiNaC Copyright (C) 1999-2003 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2008 Johannes Gutenberg University Mainz, Germany * * 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 @@ -20,7 +20,7 @@ * * 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., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ #ifndef __GINAC_NORMAL_H__ @@ -30,29 +30,52 @@ namespace GiNaC { +/** + * Flags to control the behaviour of gcd() and friends + */ +struct gcd_options +{ + enum { + /** + * Usually GiNaC tries heuristic GCD algorithm before PRS. + * Some people don't like this, so here's a flag to disable it. + */ + no_heur_gcd = 2, + /** + * GiNaC tries to avoid expanding expressions when computing + * GCDs. This is a good idea, but some people dislike it. + * Hence the flag to disable special handling of partially + * factored polynomials. DON'T SET THIS unless you *really* + * know what are you doing! + */ + no_part_factored = 4 + }; +}; + class ex; class symbol; // Quotient q(x) of polynomials a(x) and b(x) in Q[x], so that a(x)=b(x)*q(x)+r(x) -extern ex quo(const ex &a, const ex &b, const symbol &x, bool check_args = true); +extern ex quo(const ex &a, const ex &b, const ex &x, bool check_args = true); // Remainder r(x) of polynomials a(x) and b(x) in Q[x], so that a(x)=b(x)*q(x)+r(x) -extern ex rem(const ex &a, const ex &b, const symbol &x, bool check_args = true); +extern ex rem(const ex &a, const ex &b, const ex &x, bool check_args = true); // Decompose rational function a(x)=N(x)/D(x) into Q(x)+R(x)/D(x) with degree(R, x) < degree(D, x) -extern ex decomp_rational(const ex &a, const symbol &x); +extern ex decomp_rational(const ex &a, const ex &x); // Pseudo-remainder of polynomials a(x) and b(x) in Q[x] -extern ex prem(const ex &a, const ex &b, const symbol &x, bool check_args = true); +extern ex prem(const ex &a, const ex &b, const ex &x, bool check_args = true); // Pseudo-remainder of polynomials a(x) and b(x) in Q[x] -extern ex sprem(const ex &a, const ex &b, const symbol &x, bool check_args = true); +extern ex sprem(const ex &a, const ex &b, const ex &x, bool check_args = true); // Exact polynomial division of a(X) by b(X) in Q[X] (quotient returned in q), returns false when exact division fails extern bool divide(const ex &a, const ex &b, ex &q, bool check_args = true); // Polynomial GCD in Z[X], cofactors are returned in ca and cb, if desired -extern ex gcd(const ex &a, const ex &b, ex *ca = NULL, ex *cb = NULL, bool check_args = true); +extern ex gcd(const ex &a, const ex &b, ex *ca = NULL, ex *cb = NULL, + bool check_args = true, unsigned options = 0); // Polynomial LCM in Z[X] extern ex lcm(const ex &a, const ex &b, bool check_args = true); @@ -66,6 +89,9 @@ extern ex sqrfree_parfrac(const ex & a, const symbol & x); // Collect common factors in sums. extern ex collect_common_factors(const ex & e); +// Resultant of two polynomials e1,e2 with respect to symbol s. +extern ex resultant(const ex & e1, const ex & e2, const ex & s); + } // namespace GiNaC #endif // ndef __GINAC_NORMAL_H__