Fix a bug in SIMDVector under scalar mode that was causing compilation failure for when actual primitive types were mixed in the arithmetics with SIMDVector types

Author: romeric

simd_vector/simd_vector_base.h

@@ -216,18 +216,18 @@ FASTOR_INLINE SIMDVector<T,ABI> operator+(const SIMDVector<T,ABI> &a, const SIMD
         out.value[i] = a.value[i] + b.value[i];
     return out;
 }
-template<typename U, typename T, int ABI>
-FASTOR_INLINE SIMDVector<T,ABI> operator+(const SIMDVector<T,ABI> &a, U b) {
+template<typename T, int ABI>
+FASTOR_INLINE SIMDVector<T,ABI> operator+(const SIMDVector<T,ABI> &a, T b) {
     SIMDVector<T,ABI> out;
     for (FASTOR_INDEX i=0; i<SIMDVector<T,ABI>::Size; ++i)
-        out.value[i] = a.value[i] + static_cast<T>(b);
+        out.value[i] = a.value[i] + b;
     return out;
 }
-template<typename U, typename T, int ABI>
-FASTOR_INLINE SIMDVector<T,ABI> operator+(U a, const SIMDVector<T,ABI> &b) {
+template<typename T, int ABI>
+FASTOR_INLINE SIMDVector<T,ABI> operator+(T a, const SIMDVector<T,ABI> &b) {
     SIMDVector<T,ABI> out;
     for (FASTOR_INDEX i=0; i<SIMDVector<T,ABI>::Size; ++i)
-        out.value[i] = static_cast<T>(a) + b.value[i];
+        out.value[i] = a + b.value[i];
     return out;
 }
 template<typename T, int ABI>
@@ -242,18 +242,18 @@ FASTOR_INLINE SIMDVector<T,ABI> operator-(const SIMDVector<T,ABI> &a, const SIMD
         out.value[i] = a.value[i] - b.value[i];
     return out;
 }
-template<typename U, typename T, int ABI>
-FASTOR_INLINE SIMDVector<T,ABI> operator-(const SIMDVector<T,ABI> &a, U b) {
+template<typename T, int ABI>
+FASTOR_INLINE SIMDVector<T,ABI> operator-(const SIMDVector<T,ABI> &a, T b) {
     SIMDVector<T,ABI> out;
     for (FASTOR_INDEX i=0; i<SIMDVector<T,ABI>::Size; ++i)
-        out.value[i] = a.value[i] - static_cast<T>(b);
+        out.value[i] = a.value[i] - b;
     return out;
 }
-template<typename U, typename T, int ABI>
-FASTOR_INLINE SIMDVector<T,ABI> operator-(U a, const SIMDVector<T,ABI> &b) {
+template<typename T, int ABI>
+FASTOR_INLINE SIMDVector<T,ABI> operator-(T a, const SIMDVector<T,ABI> &b) {
     SIMDVector<T,ABI> out;
     for (FASTOR_INDEX i=0; i<SIMDVector<T,ABI>::Size; ++i)
-        out.value[i] = static_cast<T>(a) - b.value[i];
+        out.value[i] = a - b.value[i];
     return out;
 }
 template<typename T, int ABI>
@@ -271,18 +271,18 @@ FASTOR_INLINE SIMDVector<T,ABI> operator*(const SIMDVector<T,ABI> &a, const SIMD
         out.value[i] = a.value[i] * b.value[i];
     return out;
 }
-template<typename U, typename T, int ABI>
-FASTOR_INLINE SIMDVector<T,ABI> operator*(const SIMDVector<T,ABI> &a, U b) {
+template<typename T, int ABI>
+FASTOR_INLINE SIMDVector<T,ABI> operator*(const SIMDVector<T,ABI> &a, T b) {
     SIMDVector<T,ABI> out;
     for (FASTOR_INDEX i=0; i<SIMDVector<T,ABI>::Size; ++i)
-        out.value[i] = a.value[i] * static_cast<T>(b);
+        out.value[i] = a.value[i] * b;
     return out;
 }
-template<typename U, typename T, int ABI>
-FASTOR_INLINE SIMDVector<T,ABI> operator*(U a, const SIMDVector<T,ABI> &b) {
+template<typename T, int ABI>
+FASTOR_INLINE SIMDVector<T,ABI> operator*(T a, const SIMDVector<T,ABI> &b) {
     SIMDVector<T,ABI> out;
     for (FASTOR_INDEX i=0; i<SIMDVector<T,ABI>::Size; ++i)
-        out.value[i] = static_cast<T>(a) * b.value[i];
+        out.value[i] = a * b.value[i];
     return out;
 }
 
@@ -293,18 +293,18 @@ FASTOR_INLINE SIMDVector<T,ABI> operator/(const SIMDVector<T,ABI> &a, const SIMD
         out.value[i] = a.value[i] / b.value[i];
     return out;
 }
-template<typename U, typename T, int ABI>
-FASTOR_INLINE SIMDVector<T,ABI> operator/(const SIMDVector<T,ABI> &a, U b) {
+template<typename T, int ABI>
+FASTOR_INLINE SIMDVector<T,ABI> operator/(const SIMDVector<T,ABI> &a, T b) {
     SIMDVector<T,ABI> out;
     for (FASTOR_INDEX i=0; i<SIMDVector<T,ABI>::Size; ++i)
-        out.value[i] = a.value[i] / static_cast<T>(b);
+        out.value[i] = a.value[i] / b;
     return out;
 }
-template<typename U, typename T, int ABI>
-FASTOR_INLINE SIMDVector<T,ABI> operator/(U a, const SIMDVector<T,ABI> &b) {
+template<typename T, int ABI>
+FASTOR_INLINE SIMDVector<T,ABI> operator/(T a, const SIMDVector<T,ABI> &b) {
     SIMDVector<T,ABI> out;
     for (FASTOR_INDEX i=0; i<SIMDVector<T,ABI>::Size; ++i)
-        out.value[i] = static_cast<T>(a) / b.value[i];
+        out.value[i] = a / b.value[i];
     return out;
 }