uintpoly_flint.h

Go to the documentation of this file.

 
#ifndef SYMENGINE_UINTPOLY_FLINT_H
#define SYMENGINE_UINTPOLY_FLINT_H
 
#include <symengine/polys/upolybase.h>
 
#ifdef HAVE_SYMENGINE_FLINT
#include <symengine/flint_wrapper.h>
using fzp_t = SymEngine::fmpz_poly_wrapper;
using fqp_t = SymEngine::fmpq_poly_wrapper;
 
namespace SymEngine
{
 
template <typename Container, template <typename X, typename Y> class BaseType,
          typename Poly>
class UFlintPoly : public BaseType<Container, Poly>
{
public:
    using Cf = typename BaseType<Container, Poly>::coef_type;
 
    UFlintPoly(const RCP<const Basic> &var, Container &&dict)
        : BaseType<Container, Poly>(var, std::move(dict))
    {
    }
 
    int compare(const Basic &o) const
    {
        SYMENGINE_ASSERT(is_a<Poly>(o))
        const Poly &s = down_cast<const Poly &>(o);
 
        if (this->get_poly().degree() != s.get_poly().degree())
            return (this->get_poly().degree() < s.get_poly().degree()) ? -1 : 1;
 
        int cmp = this->get_var()->compare(*s.get_var());
        if (cmp != 0)
            return cmp;
 
        for (unsigned int i = 0; i < this->get_poly().length(); ++i) {
            if (this->get_poly().get_coeff(i) != s.get_poly().get_coeff(i))
                return (this->get_poly().get_coeff(i)
                        < s.get_poly().get_coeff(i))
                           ? -1
                           : 1;
        }
        return 0;
    }
 
    static Container container_from_dict(const RCP<const Basic> &var,
                                         std::map<unsigned, Cf> &&d)
    {
        Container f;
        for (auto const &p : d) {
            if (p.second != 0) {
                typename Container::internal_coef_type r(get_mp_t(p.second));
                f.set_coeff(p.first, r);
            }
        }
        return f;
    }
 
    static RCP<const Poly> from_vec(const RCP<const Basic> &var,
                                    const std::vector<Cf> &v)
    {
        // TODO improve this (we already know the degree)
        Container f;
        for (unsigned int i = 0; i < v.size(); i++) {
            if (v[i] != 0) {
                typename Container::internal_coef_type r(get_mp_t(v[i]));
                f.set_coeff(i, r);
            }
        }
        return make_rcp<const Poly>(var, std::move(f));
    }
 
    template <typename FromPoly>
    static enable_if_t<is_a_UPoly<FromPoly>::value, RCP<const Poly>>
    from_poly(const FromPoly &p)
    {
        Container f;
        for (auto it = p.begin(); it != p.end(); ++it) {
            typename Container::internal_coef_type r(get_mp_t(it->second));
            f.set_coeff(it->first, r);
        }
        return Poly::from_container(p.get_var(), std::move(f));
    }
 
    Cf eval(const Cf &x) const
    {
        typename Container::internal_coef_type r(get_mp_t(x));
        return to_mp_class(this->get_poly().eval(r));
    }
 
    Cf get_coeff(unsigned int x) const
    {
        return to_mp_class(this->get_poly().get_coeff(x));
    }
 
    // can't return by reference
    Cf get_coeff_ref(unsigned int x) const
    {
        return to_mp_class(this->get_poly().get_coeff(x));
    }
 
    typedef ContainerForIter<Poly, Cf> iterator;
    typedef ContainerRevIter<Poly, Cf> r_iterator;
    iterator begin() const
    {
        return iterator(this->template rcp_from_this_cast<Poly>(), 0);
    }
    iterator end() const
    {
        return iterator(this->template rcp_from_this_cast<Poly>(), size());
    }
    r_iterator obegin() const
    {
        return r_iterator(this->template rcp_from_this_cast<Poly>(),
                          (long)size() - 1);
    }
    r_iterator oend() const
    {
        return r_iterator(this->template rcp_from_this_cast<Poly>(), -1);
    }
 
    int size() const
    {
        return numeric_cast<int>(this->get_poly().length());
    }
};
 
class UIntPolyFlint : public UFlintPoly<fzp_t, UIntPolyBase, UIntPolyFlint>
{
public:
    IMPLEMENT_TYPEID(SYMENGINE_UINTPOLYFLINT)
    UIntPolyFlint(const RCP<const Basic> &var, fzp_t &&dict);
    hash_t __hash__() const override;
 
}; // UIntPolyFlint
 
class URatPolyFlint : public UFlintPoly<fqp_t, URatPolyBase, URatPolyFlint>
{
public:
    IMPLEMENT_TYPEID(SYMENGINE_URATPOLYFLINT)
    URatPolyFlint(const RCP<const Basic> &var, fqp_t &&dict);
    hash_t __hash__() const override;
}; // URatPolyFlint
 
template <typename Container, template <typename X, typename Y> class BaseType,
          typename Poly>
std::vector<std::pair<RCP<const Poly>, long>>
factors(const UFlintPoly<Container, BaseType, Poly> &a)
{
    auto fac_wrapper = a.get_poly().factors();
    auto &fac = fac_wrapper.get_fmpz_poly_factor_t();
    std::vector<std::pair<RCP<const Poly>, long>> S;
    if (fac->c != 1_z)
        S.push_back(std::make_pair(
            make_rcp<const Poly>(a.get_var(), numeric_cast<int>(fac->c)), 1));
    SYMENGINE_ASSERT(fac->p != NULL and fac->exp != NULL)
    for (long i = 0; i < fac->num; i++) {
        fzp_t z;
        z.swap_fmpz_poly_t(fac->p[i]);
        S.push_back(std::make_pair(
            make_rcp<const Poly>(a.get_var(), std::move(z)), fac->exp[i]));
    }
    return S;
}
 
template <typename Container, template <typename X, typename Y> class BaseType,
          typename Poly>
RCP<const Poly> gcd_upoly(const UFlintPoly<Container, BaseType, Poly> &a,
                          const Poly &b)
{
    if (!(a.get_var()->__eq__(*b.get_var())))
        throw SymEngineException("Error: variables must agree.");
    return make_rcp<const Poly>(a.get_var(), a.get_poly().gcd(b.get_poly()));
}
 
template <typename Container, template <typename X, typename Y> class BaseType,
          typename Poly>
RCP<const Poly> lcm_upoly(const UFlintPoly<Container, BaseType, Poly> &a,
                          const Poly &b)
{
    if (!(a.get_var()->__eq__(*b.get_var())))
        throw SymEngineException("Error: variables must agree.");
    return make_rcp<const Poly>(a.get_var(), a.get_poly().lcm(b.get_poly()));
}
 
template <typename Container, template <typename X, typename Y> class BaseType,
          typename Poly>
RCP<const Poly> pow_upoly(const UFlintPoly<Container, BaseType, Poly> &a,
                          unsigned int p)
{
    return make_rcp<const Poly>(a.get_var(), a.get_poly().pow(p));
}
 
template <typename Container, template <typename X, typename Y> class BaseType,
          typename Poly>
bool divides_upoly(const UFlintPoly<Container, BaseType, Poly> &a,
                   const Poly &b, const Ptr<RCP<const Poly>> &res)
{
    if (!(a.get_var()->__eq__(*b.get_var())))
        throw SymEngineException("Error: variables must agree.");
 
    typename Poly::container_type q, r;
 
    b.get_poly().divrem(q, r, a.get_poly());
    if (r == 0) {
        *res = make_rcp<Poly>(a.get_var(), std::move(q));
        return true;
    } else {
        return false;
    }
}
} // namespace SymEngine
 
#endif // HAVE_SYMENGINE_FLINT
 
#endif // SYMENGINE_UINTPOLY_FLINT_H