X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Fnormal.h;h=ff418bbea617fc39520fc957278c90efe04bf170;hp=b461dcf0723407930f2b9ca9e9d1ade6f50e9bc5;hb=e5e2ac399c225a8c7248e1879e192eba00067538;hpb=4e3a4ac2bcb0837611ea31bc8fc05d84a20c33ac diff --git a/ginac/normal.h b/ginac/normal.h index b461dcf0..ff418bbe 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-2001 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2004 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 @@ -26,19 +26,27 @@ #ifndef __GINAC_NORMAL_H__ #define __GINAC_NORMAL_H__ +#include "lst.h" + namespace GiNaC { 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 ex &x); + +// Pseudo-remainder of polynomials a(x) and b(x) in Q[x] +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 Z[x] -extern ex prem(const ex &a, const ex &b, const symbol &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 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); @@ -50,7 +58,16 @@ extern ex gcd(const ex &a, const ex &b, ex *ca = NULL, ex *cb = NULL, bool check extern ex lcm(const ex &a, const ex &b, bool check_args = true); // Square-free factorization of a polynomial a(x) -extern ex sqrfree(const ex &a, const symbol &x); +extern ex sqrfree(const ex &a, const lst &l = lst()); + +// Square-free partial fraction decomposition of a rational function a(x) +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