include/algebramix/polynomial_kronecker.hpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : polynomial_kronecker.hpp
* DESCRIPTION: Multiplication of polynomials via Kronecker substitution
* COPYRIGHT  : (C) 2009  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_KRONECKER_HPP
#define __MMX_POLYNOMIAL_KRONECKER_HPP
#include <algebramix/polynomial_naive.hpp>

namespace mmx {
#define TMPL template<typename C> 

/******************************************************************************
* Kronecker variant (for the product)
******************************************************************************/

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

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

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

template<typename V,typename W>
struct implementation<polynomial_multiply,V,polynomial_kronecker<W> >:
  public implementation<polynomial_linear,V>
{
  typedef polynomial_multiply_threshold<polynomial_kronecker<W> > Th;
  typedef implementation<polynomial_multiply,W> Fallback;

template<typename C,typename K> static inline void
mul (C* dest, const C* s1, const K* s2, nat n1, nat n2) {
  Fallback::mul (dest, s1, s2, n1, n2); }

TMPL static inline void
tmul (C* dest, const C* s1, const C* s2, nat n1, nat n2) {
  Fallback::tmul (dest, s1, s2, n1, n2); }

TMPL static inline void
mul (C* dest, const C* s1, const C* s2, nat n1, nat n2) {
  if (min (n1, n2) < Threshold(C,Th))
    Fallback::mul (dest, s1, s2, n1, n2);
  else mul_kronecker (dest, s1, n1, s2, n2); }

TMPL static inline void
square (C* dest, const C* s, nat n) {
  if (n < Threshold(C,Th)) Fallback::square (dest, s, n);
  else square_kronecker (dest, s, n); }

}; // implementation<polynomial_multiply,V,polynomial_kronecker<W> >

#undef TMPL
} // namespace mmx
#endif // __MMX_POLYNOMIAL_KRONECKER_HPP

---