littrack/parallel/vectors/PSEUDO_F32vec1.h

Go to the documentation of this file.

#ifndef LIT_F32VEC1_H
#define LIT_F32VEC1_H
 
#include <cmath>
#include <iostream>
 
/**********************************
 *
 *   Vector of one single float
 *
 **********************************/
 
//const union {
//   int i[1];
//   float m;
//}
//__f32vec1_true_cheat     = {0xFFFFFFFF},
//__f32vec1_false_cheat    = {0x00000000};
 
//#define _f32vec1_true     ((F32vec1)__f32vec1_true_cheat.m)
//#define _f32vec1_false    ((F32vec1)__f32vec1_false_cheat.m)
 
class F32vec1 {
public:
  float v;
 
  float& operator[](int i) { return v; }
  float operator[](int i) const { return v; }
 
  F32vec1() {}
  F32vec1(const float& v0) { v = v0; }
 
  /* Conversion function */
  operator float() const { return v; } /* Convert to __m128 */
 
  /* Arithmetic Operators */
 
  /* Functions */
  //friend F32vec1 min( const F32vec1 &a, const F32vec1 &b ){ return a<b ?a :b; }
  //friend F32vec1 max( const F32vec1 &a, const F32vec1 &b ){ return a>b ?a :b; }
 
  /* Square Root */
 
  /* Reciprocal( inverse) Square Root */
  //friend F32vec1 rsqrt( const F32vec1 &a ){ return 1./sqrt(a); }
 
  /* Reciprocal (inversion) */
  friend F32vec1 rcp(const F32vec1& a) { return 1. / a; }
 
  /* Absolute value */
  //friend F32vec1 fabs(const F32vec1 &a){ return fabs(a); }
 
  /* Sign */
  //friend F32vec1 sgn(const F32vec1 &a){ return a<0 ?-1 :(a>0 ?1 :0); }
 
  /* Logical */
  /*
   friend F32vec1 operator&( const F32vec1 &a, const F32vec1 &b ){ // mask returned
     F32vec1 tmp;
     int *x = (int*)&tmp;
     int *y = (int*)&a;
     int *z = (int*)&b;
     x[0] = y[0] & z[0];
     x[1] = y[1] & z[1];
     return tmp;
   }
   */
  /* Non intrinsic functions */
 
  /* Define all operators for consistensy */
 
  friend void operator+=(F32vec1& a, const F32vec1& b) { a = a + b; }
  friend void operator-=(F32vec1& a, const F32vec1& b) { a = a - b; }
  friend void operator*=(F32vec1& a, const F32vec1& b) { a = a * b; }
  friend void operator/=(F32vec1& a, const F32vec1& b) { a = a / b; }
 
  friend std::ostream& operator<<(std::ostream& strm, const F32vec1& a) {
    strm << a[0];
    return strm;
  }
 
  friend std::istream& operator>>(std::istream& strm, F32vec1& a) {
    float tmp;
    strm >> tmp;
    a = tmp;
    return strm;
  }
 
};  //__attribute__ ((aligned(16)));;
 
typedef F32vec1 fvec;
const int fvecLen = 1;
typedef float fscal;
//#define fvec_true  _f32vec1_true
//#define fvec_false _f32vec1_false
#define _fvecalignment
 
#endif