include/algebramix/vector_simd.hpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : vector_simd.hpp
* DESCRIPTION: vectors using SIMD instructions
* COPYRIGHT  : (C) 2008  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__VECTOR_SIMD__HPP
#define __MMX__VECTOR_SIMD__HPP
#include <numerix/simd.hpp>
#include <algebramix/vector_unrolled.hpp>
#include <algebramix/vector_aligned.hpp>
#include <algebramix/vector_sse.hpp>

namespace mmx {

/******************************************************************************
* Variant for SIMD support
******************************************************************************/

template<nat len1, nat len2, typename V= vector_naive>
struct vector_simd: public vector_aligned<
                             vector_unrolled<len1,V>,
                             vector_unrolled<len2,V> >
{
  typedef vector_simd<len1,len2,typename V::Naive> Naive;
  typedef vector_simd<len1,len2,typename V::Aligned> Aligned;
  typedef vector_unrolled<len1,typename V::No_simd> No_simd;
  typedef vector_simd<len1,len2,typename V::No_thread> No_thread;
};

template<typename F, typename V, nat len1, nat len2, typename W>
struct implementation<F,V,vector_simd<len1,len2,W> >:
  public implementation<F,V,vector_aligned<
                              vector_unrolled<len1,W>,
                              vector_unrolled<len2,W> > > {};

/******************************************************************************
* Suggested specilizations for an Intel plateform supporting SSE2
******************************************************************************/

#define Vector_simd_variant(C) vector_simd_variant_helper<C>::VV

template<typename C>
struct vector_simd_variant_helper {
  typedef vector_simd<4,4> VV; };

STMPL
struct vector_simd_variant_helper<double> {
  typedef vector_simd<8,4> VV; };

STMPL
struct vector_simd_variant_helper<char> {
  typedef vector_simd<32,4> VV; };

STMPL
struct vector_simd_variant_helper<signed char> {
  typedef vector_simd<32,4> VV; };

STMPL
struct vector_simd_variant_helper<unsigned char> {
  typedef vector_simd<32,4> VV; };

STMPL
struct vector_simd_variant_helper<short int> {
  typedef vector_simd<32,4 > VV; };

STMPL
struct vector_simd_variant_helper<short unsigned int> {
  typedef vector_simd<32, 4> VV; };

STMPL
struct vector_simd_variant_helper<int> {
  typedef vector_simd<16, 4> VV; };

STMPL
struct vector_simd_variant_helper<unsigned int> {
  typedef vector_simd<16, 4> VV; };

STMPL
struct vector_simd_variant_helper<long int> {
  typedef vector_simd<16, 4> VV; };

STMPL
struct vector_simd_variant_helper<long unsigned int> {
  typedef vector_simd<16, 4> VV; };

STMPL
struct vector_simd_variant_helper<long long int> {
  typedef vector_simd<8,4> VV; };

STMPL
struct vector_simd_variant_helper<long long unsigned int> {
  typedef vector_simd<8,4> VV;
};

} // namespace mmx

#endif // __MMX__VECTOR_SIMD__HPP

---