include/algebramix/matrix_int.hpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : matrix_int.hpp
* DESCRIPTION: matrices over hardware integers
* 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_MATRIX_INT_HPP
#define __MMX_MATRIX_INT_HPP
#include <algebramix/vector_unrolled.hpp>
#include <algebramix/vector_simd.hpp>
#include <algebramix/matrix.hpp>
#include <algebramix/matrix_ring_naive.hpp>
#include <algebramix/matrix_strassen.hpp>
#include <algebramix/matrix_threads.hpp>
#include <algebramix/matrix_simd.hpp>
namespace mmx {
#define TMPL template<typename V>
#define Matrix matrix<int,V>
/******************************************************************************
* Special case for matrix<int>
******************************************************************************/

TMPL inline Matrix operator + (const Matrix& x, const int& y) {
  return x + promote (y, x); }
TMPL inline Matrix operator + (const int& x, const Matrix& y) {
  return promote (x, y) + y; }
TMPL inline Matrix operator - (const Matrix& x, const int& y) {
  return x - promote (y, x); }
TMPL inline Matrix operator - (const int& x, const Matrix& y) {
  return promote (x, y) - y; }
TMPL inline Matrix operator * (const Matrix& x, const int& y) {
  return x * promote (y, x); }
TMPL inline Matrix operator * (const int& x, const Matrix& y) {
  return promote (x, y) * y; }
TMPL inline Matrix operator / (const Matrix& x, const int& y) {
  return x / promote (y, x); }
TMPL inline Matrix operator / (const int& x, const Matrix& y) {
  return promote (x, y) / y; }

#undef TMPL
#undef Matrix

/******************************************************************************
* Roughly 'optimized' variants for an '__SSE2__' Intel processor
******************************************************************************/

DEFINE_VARIANT_1 (typename I, I, matrix_int_simd,
                    matrix_strassen<
                        typename Matrix_simd_variant(I) >)

DEFINE_VARIANT (matrix_int,
                matrix_ring_naive<
                  matrix_strassen<
                    matrix_unrolled<16> > >)

STMPL
struct matrix_variant_helper<signed char> {
  typedef matrix_int_simd<signed char> MV; };
STMPL
struct matrix_variant_helper<unsigned char> {
  typedef matrix_int_simd<unsigned char> MV; };
STMPL
struct matrix_variant_helper<short int> {
  typedef matrix_int_simd<short int> MV; };
STMPL
struct matrix_variant_helper<unsigned short int> {
  typedef matrix_int_simd<unsigned short int> MV; };
STMPL
struct matrix_variant_helper<int> {
  typedef matrix_int MV; };
STMPL
struct matrix_variant_helper<unsigned int> {
  typedef matrix_int MV; };
STMPL
struct matrix_variant_helper<long int> {
  typedef matrix_int MV; };
STMPL
struct matrix_variant_helper<unsigned long int> {
  typedef matrix_int MV; };
STMPL
struct matrix_variant_helper<long long int> {
  typedef matrix_int MV; };
STMPL
struct matrix_variant_helper<unsigned long long int> {
  typedef matrix_int MV; };
}
#endif // __MMX_MATRIX_INT_HPP

---