/Users/mourrain/Devel/mmx/algebramix/include/algebramix/polynomial_polynomial.hpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : polynomial_polynomial.hpp
* DESCRIPTION: Kronecker multiplication of bivariate polynomials
* COPYRIGHT  : (C) 2004  Joris van der Hoeven and Gregoire Lecerf
*******************************************************************************
* This software falls under the GNU general public license and comes WITHOUT
* ANY WARRANTY WHATSOEVER. See the file $TEXMACS_PATH/LICENSE for more details.
* If you don't have this file, write to the Free Software Foundation, Inc.,
* 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
******************************************************************************/

#ifndef __MMX_POLYNOMIAL_POLYNOMIAL_HPP
#define __MMX_POLYNOMIAL_POLYNOMIAL_HPP
#include <algebramix/polynomial.hpp>
#include <algebramix/polynomial_ring_dicho.hpp>
#include <algebramix/kronecker_polynomial.hpp>
#include <algebramix/polynomial_kronecker.hpp>

namespace mmx {

/******************************************************************************
* Variant
******************************************************************************/

template<typename V>
struct polynomial_polynomial: public V {
  typedef typename V::Vec Vec;
  typedef typename V::Naive Naive;
  typedef polynomial_polynomial<typename V::Positive> Positive;
  typedef polynomial_polynomial<typename V::No_simd> No_simd;
  typedef polynomial_polynomial<typename V::No_thread> No_thread;
  typedef polynomial_polynomial<typename V::No_scaled> No_scaled;
};

template<typename F, typename V, typename W>
struct implementation<F,V,polynomial_polynomial<W> >:
  public implementation<F,V,W> {};

/******************************************************************************
* Variable names, default is C[x][y][z]
******************************************************************************/

template<typename V,typename W>
struct implementation<polynomial_defaults,V,polynomial_polynomial<W> > {

  template<typename P>
  class global_variables {
    static inline generic& dyn_name () {
      static generic name = "x";
      return name; }
  public:
    static inline void set_variable_name (const generic& x) { dyn_name () = x; }
    static inline generic get_variable_name () { return dyn_name (); }
  };

  template<typename C, typename V1, typename V2>
  class global_variables<polynomial<polynomial<C,V1>, V2> > {
    static inline generic& dyn_name () {
      static generic name = "y";
      return name; }
  public:
    static inline void set_variable_name (const generic& x) {
      dyn_name () = x; }
    static inline generic get_variable_name () {
      return dyn_name (); }
  };

  template<typename C, typename V1, typename V2, typename V3>
  class global_variables<polynomial<polynomial<
                                           polynomial<C,V1>, V2>, V3 > > {
    static inline generic& dyn_name () {
      static generic name = "z";
      return name; }
  public:
    static inline void set_variable_name (const generic& x) {
      dyn_name () = x; }
    static inline generic get_variable_name () {
      return dyn_name (); }
  };

};

/******************************************************************************
* Kronecker exact division
******************************************************************************/

template<typename V,typename W>
struct implementation<polynomial_exact_divide,V,
                      polynomial_polynomial<W> >:
  public implementation<polynomial_linear,V>
{
  typedef polynomial_exact_divide_threshold <polynomial_polynomial<W> > Th;
  typedef implementation<polynomial_exact_divide,W> Fallback;

template<typename C> static inline void
div (C* dest, const C* s1, const C* s2, nat n1, nat n2) {
  if (min (n1, n2) < Threshold(C,Th))
    Fallback::div (dest, s1, s2, n1, n2);
  else div_kronecker (dest, s1, n1, s2, n2); }

};

/******************************************************************************
* Kronecker multiplication
******************************************************************************/

template<typename C, typename V>
struct polynomial_variant_helper<polynomial<C,V> > {
  typedef polynomial_gcd_ring_dicho<
            polynomial_kronecker<
              polynomial_polynomial<
                polynomial_dicho<
                  polynomial_naive> > > > PV;
};

template<typename C, typename V1, typename V2>
struct polynomial_variant_helper<polynomial<polynomial<C,V1>,V2> > {
  typedef V2 PV;
};

} // namespace mmx
#endif // __MMX_POLYNOMIAL_POLYNOMIAL_HPP