include/algebramix/matrix_integer.hpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : matrix_integer.hpp
* DESCRIPTION: Multi-modular multiplication of integer matrices
* COPYRIGHT  : (C) 2009  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_MATRIX_INTEGER_HPP
#define __MMX_MATRIX_INTEGER_HPP
#include <numerix/integer.hpp>
#include <algebramix/crt_integer.hpp>
#include <algebramix/matrix.hpp>
#include <algebramix/matrix_modular_int.hpp>
#include <algebramix/matrix_ring_naive.hpp>
#include <algebramix/matrix_crt.hpp>
#include <algebramix/matrix_balanced.hpp>
#include <algebramix/matrix_bareiss.hpp>

namespace mmx {

/******************************************************************************
* Use multi-modular multiplication by default for integer matrices
******************************************************************************/

DEFINE_VARIANT (matrix_integer,
                matrix_balanced<
                  matrix_bareiss<
                    matrix_crt<
                      matrix_ring_naive<matrix_naive> > > >)

STMPL
struct matrix_variant_helper<integer> {
  typedef matrix_integer MV;
};

/******************************************************************************
* Size helper and default moduli
******************************************************************************/

xnat
matrix_product_bit_size (const integer* s1, nat s1_rs, nat s1_cs,
                         const integer* s2, nat s2_rs, nat s2_cs,
                         nat r, nat l, nat c);

STMPL
struct matrix_crt_multiply_helper<integer> {
private:
  struct mat_crt_naive {
    typedef integer base;
    typedef int modulus_base;
    typedef modulus_int_preinverse<20> modulus_base_variant;
    typedef Modulus_variant(integer) modulus_variant;
  };
  struct mat_crt_dicho {
    typedef integer base;
    typedef integer modulus_base;
    typedef modulus_integer_naive modulus_base_variant;
    typedef Modulus_variant(integer) modulus_variant;
  };
  typedef modulus<mat_crt_naive::modulus_base,
                  mat_crt_naive::modulus_base_variant> M;
public:
  static const nat dimension_threshold= 7;
  static const nat ratio_threshold= 100; // this is a percentage

  typedef crt_naive_transformer<integer,mat_crt_naive> crt_naive;
  typedef crt_dicho_transformer<integer,mat_crt_dicho> crt_dicho;
  typedef crt_blocks_transformer<crt_naive,crt_dicho,160> crt_transformer;

  typedef moduli_helper<integer,M,fft_prime_sequence_int<20> > moduli_sequence;

  static inline xnat size (const integer* s1, nat s1_rs, nat s1_cs,
                           const integer* s2, nat s2_rs, nat s2_cs,
                           nat r, nat l, nat c) {
    return matrix_product_bit_size (s1, s1_rs, s1_cs,
                                    s2, s2_rs, s2_cs,
                                    r, l, c); }
};

template<typename W>
struct threshold_helper<integer,matrix_multiply_threshold<W> > {
  typedef fixed_value<nat,2> impl;
};

} // namespace mmx
#endif // __MMX_MATRIX_INTEGER_HPP

---