include/algebramix/matrix_simd.hpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : matrix_simd.hpp
* DESCRIPTION: matrices using SIMD instructions
* COPYRIGHT  : (C) 2007  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_SIMD__HPP
#define __MMX__MATRIX_SIMD__HPP
#include <algebramix/vector_simd.hpp>
#include <algebramix/matrix_unrolled.hpp>
#include <algebramix/matrix_aligned.hpp>
#include <algebramix/matrix_sse.hpp>

namespace mmx {

/******************************************************************************
* SIMD matrix variant
******************************************************************************/

template<nat len1, nat len2, typename V=matrix_naive>
struct matrix_simd: public matrix_aligned<
                             matrix_unrolled<len1, V>,
                             matrix_unrolled<len2, V> >
{
  typedef vector_simd<len1,len2,typename V::Vec> Vec;
  typedef matrix_simd<len1,len2,typename V::Naive> Naive;
  typedef matrix_simd<len1,len2,typename V::Positive> Positive;
  typedef matrix_unrolled<len1,typename V::No_simd> No_simd;
  typedef matrix_simd<len1,len2,typename V::No_thread> No_thread;
  typedef matrix_simd<len1,len2,typename V::No_scaled> No_scaled;
};

template<typename F, typename V, nat len1, nat len2, typename W>
struct implementation<F,V,matrix_simd<len1,len2,W> >:
  public implementation<F,V,matrix_aligned<
                              matrix_unrolled<len1, W>,
                              matrix_unrolled<len2, W> > > {};

/******************************************************************************
* Suggested specializations
******************************************************************************/

#define Matrix_simd_variant(C) matrix_simd_variant_helper<C>::MV

template<typename C>
struct matrix_simd_variant_helper {
  typedef matrix_simd<4,4> MV; };

STMPL
struct matrix_simd_variant_helper<double> {
  typedef matrix_simd<8,4> MV; };

STMPL
struct matrix_simd_variant_helper<char> {
  typedef matrix_simd<16,8> MV; };

STMPL
struct matrix_simd_variant_helper<signed char> {
  typedef matrix_simd<16,8> MV; };

STMPL
struct matrix_simd_variant_helper<unsigned char> {
  typedef matrix_simd<16,4> MV; };

STMPL
struct matrix_simd_variant_helper<short int> {
  typedef matrix_simd<16,4 > MV; };

STMPL
struct matrix_simd_variant_helper<short unsigned int> {
  typedef matrix_simd<32, 4> MV; };

STMPL
struct matrix_simd_variant_helper<int> {
  typedef matrix_simd<16, 4> MV; };

STMPL
struct matrix_simd_variant_helper<unsigned int> {
  typedef matrix_simd<16, 4> MV; };

STMPL
struct matrix_simd_variant_helper<long int> {
  typedef matrix_simd<16, 4> MV; };

STMPL
struct matrix_simd_variant_helper<long unsigned int> {
  typedef matrix_simd<16, 4> MV; };

STMPL
struct matrix_simd_variant_helper<long long int> {
  typedef matrix_simd<8,4> MV; };

STMPL
struct matrix_simd_variant_helper<long long unsigned int> {
  typedef matrix_simd<8,4> MV;
};

} // namespace mmx
#endif // __MMX__MATRIX_SIMD__HPP

---