include/algebramix/crt_polynomial.hpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : crt_polynomial.hpp
* DESCRIPTION: Crt support for polynomials
* 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_CRT_POLYNOMIAL_HPP
#define __MMX_CRT_POLYNOMIAL_HPP
#include <algebramix/modular_polynomial.hpp>
#include <algebramix/crt_blocks.hpp>
namespace mmx {

/******************************************************************************
* Default moduli sequence
******************************************************************************/

struct coprime_moduli_sequence_polynomial {
  template<typename C, typename V, typename W, typename Z> static bool
  extend (vector<modulus<polynomial<C,V>, W>, Z>& v, nat k) {
    static polynomial<C,V> x (C(1), 1);
    for (nat i= N(v); i < k; i++) {
        if (i != 0 && C (i) == 0) return false;
        v << modulus<polynomial<C,V>, W> (x - C (i));
    }
    return true; }
};

template<typename C, typename V, typename W>
struct coprime_moduli_helper<modulus<polynomial<C,V>, W> > {
  typedef coprime_moduli_sequence_polynomial sequence;
};

/******************************************************************************
* Default transformer
******************************************************************************/

#define Crt_polynomial_transformer(P)\
  crt_blocks_transformer<\
                         crt_naive_transformer<P>,\
                         crt_dicho_transformer<P> >

template<typename C, typename V>
struct threshold_helper<polynomial<C,V>,crt_blocks_threshold> {
  typedef fixed_value<nat,32> impl;
};

} // namespace mmx
#endif // __MMX_CRT_POLYNOMIAL_HPP

---